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NOMENCLATURE 


A, B, C 
C 

U 

D 

d 

E 

G 

K 

k 

L 

9 . 

p, a, r 
q 

Re 

->■ 

v 

V = (u, V, w) 

x, r, 9 

x, y, z 
Z 
a 
e 


calibration constants 
k-e model constant 
test section diameter 
inlet nozzle diameter 
hot-wire voltage 
pitch factor 
yaw factor 

kinetic energy of turbulence 
contraction nozzle downstream distance 
turbulence dissipation length scale 
selected hot-wire probe positions 
local velocity vector 
Reynolds number 

time-mean vector velocity magnitude 
time-mf»3n velocity (in x- , r-, 6-directions) in 
facility coordinates 

axial, radial, and azimuthal cylindrical polar coordi 
nates 

Cartesian coordinates 

effective cooling velocity acting on a wire 
side-wall expansion angle 
turbulence energy dissipation rate 
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* 

M 

P 


correlation coefficient (estimated) between cooling 
velocities of adjacent wire orientations 
swirl vane angle with respect to facility axis 
turbulent viscosity 
time-mean density 

Subscripts 


1,2,3, 4,5, 6 
P, 0., R 
i , J 
o 

rms 


refer to the six probe measuring positions 
refer tc the three selected cooling velocities 
summation indices 
value at inlet to flowfield 
root-mean-squared quantity 


( ) 

( )' 
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Superscri pts 
time-mean average 
fluctuating quantity 
relative to probe coordinates 
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CHAPTER I 


INTRODUCTION 

1.1 Combustor Design Phenomena 

The function of the gas turbine combustion chamber, illustrated in 
Figure I, is to heat the air from the compressor up to the design maxi- 
mum temperature. That is, the temperature that the combustion chamber 
and turbine will stand without undue creep or distortion for a reason- 
able period of time. The combustor must also give reasonably uniform 
velocity and temperature distributions onto the leading edge of the tur- 
bine biades. Complete combustion to minimize the formation of carbon 
particles and pollutants end ignition reliability at high altitudes are 
both very important criteria for the efficient performance of a gas tur- 
bine combustor. 

To improve c !ame stabilization and promote mixing between the incom- 
ing fuel and the hot products of combustion, swirl is often imparted to 
the incoming combustion air prior to entry into the primary zone. The 
purpose of swirl is to assist in the creation of a toroidal flow reversal 
to obtain high rates of mixing. This type of recirculation provides bet- 
ter mixing than might be obtained by other means, such as bluff bodies 
[!]. 

In design situations, the engineer has to analyze the swirling re- 
circulating turbulent flow within combustors whose complexity is increas- 
ed manyfold by the processes of combustion and heat transfer. At present 
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the combustor designer is faced with making only slight modifications to 
present combustors because of the economic limitations of full-scale test- 
ing of combustion chambers. Before a route can be provided which leads to 
the accomplishments of design objectives more quickly and less expensive- 
ly than current practice permits greater knowledge is required concerning 
the flowfields inside gas turbine combustion chambers. 

1.2 Tne Present Research 

In gas turbine and ramjet combustion chamber development, designers 
are aided by both experimental and theoretical studies [2], Both methods 
of approach are complementary, theoretical modeling being aided by care- 
fully chosen experiments [ 3 ]. A major area of research need, with which 
the present research work is concerned, is the measurement of time-mean 
and turbulence properties in nonreacting swirling flow in various combus- 
tor geometries, so as to aid in the understanding and development of sui t- 
able turbulence models [4]. 

Some recent researchers [5, 6, 7] have measured the time-mean and 
turbulence properties of swirling flowfields in the presence of chemical 
reaction. However, a more fundamental knowledge of the fluctuating velo- 
cities and their cross-correlations for a variety of swirl strengths 
under isothermal conditions is needed before the complex interactions of 
chemical reaction and turbulent mixing are fully understood. The present 
research is restricted to isothermal flows. 

The specific problem being investigated experimentally is concerned 
with the steady turbulent flow of air in axi symmetr i c geometries under 
low speed conditions. Specifically, flow exits from a round-sectioned 
contoured nozzle with known velocity and turbulence distributions and 
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enters a round pipe with an expansion ratio D/d *= 2 or I, as illustrated 
in Figure 2. The incoming air may possess a swiri component of velocity 
via passage through an annular vane swirler with blades at angle 4 1 to 
the undisturbed flow direction. If swirl is present, the resulting flow- 
field domain may possess a central toroidal recirculation zone (CTRZ) . 

In the sudden expansion case a corner recirculation zone (CRZ) is provok- 
ed because of the abrupt change in geometry of the boundary. It is well 
known that both of these recirculation zones contain large turbulent 
eddies which promote turbulent mixing. To simulate ramjet combustion 
chambers, the test section may be equipped with a strong contraction noz- 
zle of area ratio 

Of vital concern is the characterization of these turbulent, swirl- 
ing, confined jet flows in terms of the effects of swirl strength, expan 
sion ratio, and downstream blockage on the time-mean velocities and Rey- 
nolds normal and shear stresses. Realization that these flows also play 
an important role in many other engineering applications has inspired t 
definite desire for investigation in this area. 

1.3 Experimental Investigation 

Measurements of the time-mean and turbulence characteristics of con 
fined swirling jet flows are being used to aid in the evolution of turbu 
lence models for these complex flow situations. These measurements are 
part of a larger ongoing project at Oklahoma State University into the 
fundamental concepts of swirling flows in combustor geometries. 

The present research has included flow visualization techniques us- 
ing kerosene smoke and neutral ly-buoyant helium-filled soap bubbles re- 
corded on still [8] and movie [9] photographs. Time-mean velocities in 
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swirling and nonswirling confined flowfields have been found by use of a 
five-hole pitot probe [8, 10]. Turbulence properties using the single 
hot-wire mul t i -or i entat ion technique [11] and a crossed hot-wire probe 
[12] have been measured In nonswirling flowfields. The performance of 
the annular vane swirler used in this report has also been studied using 
a five-hole pitot probe [13]. Predictions of the flowfields correspond- 
ing to the ones studied have been made using the measured inlet condi- 
tions [14, 15, 16]. The predictions are achieved by using the STARPIC 
computer code which solves the appropriate partial differential equa- 
tions and simulates turbulence by means of the two-equation k-e turbu- 
1 ence mode I [4] . 

All this previous work has been performed using the low speed wind 
tunnel, as seen in Figure 3, in the Mechanical Engineering laboratory at 
Oklahoma State University. The facility consists of an axial-flow fan 
whose speed can be varied. Numerous fine screens and straws produce flow 
in the settling chamber of relatively low turbulence intensity. The con- 
toured nozzle leading to the test section has been designed to produce a 
minimum adverse pressure gradient on the boundary layer and thus avoid 
any regions of separated flow within the nozzle. The test section con- 
sists of either a 15 cm or 30 cm diameter plexiglass tube depending on 
which expansion ratio is desired. Plexiglass was chosen so as to allow 
accurate probe positioning and permit flow visualization. The substan- 
tial size of the test section provides excellent probe resolution. The 
swirler, placed upstream of the test section, consists of ten variable 
angle flat blades with a pi tch- to-chord ratio of 0.68. Each individual 
blade can be set to any desired angle between 0 and 80 degrees. Three 
swirl flow types are investigated in the present study: nonswirling 
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flow ($ = 0° swirler removed), moderately swirling flow ($ = 38 and 45°) 
and strongly swirling Flow ($ = 60 and 70°). 

1.4 Scope and Objectives 

The focus of the present research is the measurement of time-mean 
and turbulence properties of a confined swirling jet. The measurements 
provide a data base for the development of turbulence models for these 
flow situations. Measurements have been made for a full range of swirl 
strengths and geometric configurations over the entire flowfieldof inter- 
est. Specifically, these objectives include: 

1. Measurements of time-mean velocities and complete Reynolds 
stress tensor to find the effect of swirl on a suddenly expanding con- 
fined jet with an expansion ratio of 2. Swirl vane angles of 0 (swirler 
removed), 38 , 45, 60, and 70 degrees are considered. 

2. The effect of a strong contraction nozzle with an area reduc- 
tion ratio of 4 located at two large chamber diameters downstream of the 
jet exit is investigated. Time-mean velocities and turbulent stresses 
are measured for three cases of swirl, <j) => 0, 45, and 70 degrees. 

3- A directional sensitivity analysis is performed on the measur- 
ing technique to determine the relative accuracy of the technique to ap- 
proach time- mean velocity. 

4. Measurements of time-mean velocities and Reynolds stress tensor 
in a nonexpanding confined jet flowfield (D/d = 1) to find the effect of 
the sudden expansion on three cases of swirl, $ = 0, 45, and 70 degrees. 

5. Determination of the length scale of energy containing eddies 
in suddenly expanding and nonexpanding confined jets . The swirl strengths 
of interest are nonswirling, medium, and strong swirl. 


CHAPTER II 


REVIEW OF PREVIOUS EXPERIMENTAL STUDIES 
IN RECIRCULATING FLOWS 

Previous experimental work in axisymmetric recirculating flows pro- 
vide an important background for this investigation. Information regard- 
ing nonswirling and swirling turbulent flows have been extensively re- 
viewed. Techniques, in the measurement of turbulence are constantly be- 
ing updated; these require reviewing for their applicability to the 
present study. Noteworthy results are included in the summary where the 
emphasis is on nonreacting turbulent swirling flows. 

2.1 Confined Nonswirling Flows 

One of the first researchers to make a thorough investigation into 
confined nonswirling jet flow was Chaturvedi [17), who measured the mean 
velocity and pressure and characteristics of the turbulence in a confin- 
ed jet with an expansion ratio D/d = 2. The mean velocity and pressure 
measurements were made by pitot probe and stagnation tube, respectively. 
Measurements of velocity head were made using a pitot tube of very small 
diameter with a 15 degree taper machined internally on the leading edge. 
This was due to regions of high turbulence intensity or where the direc- 
tion of the velocity vector was unknown. The turbulence characteristics 
were measured first with a single hot-wire to determine u' and again 
with cross-wires to determine the radial intensity v'; the measurements 
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also produced the turbulent shear u'v', It was reported trial the corner 
recirculation zone (CRZ) extended to a length of two chamber diameters 
downstream of the jet exit for the sudden expansion case with high turbu- 
lence intensities on the edge of this zone. 

Many researchers have studied the recirculation zone behind a pipe 
step [18-22], in particular, the reattachment length. Krall and Sparrow 
[18] utilized an orifice in an electrically heated tube to create flow 
separation with water as the modeling fluid. It was presumed that the 
peak Nusselt number occurred at the reattachment point which they found 
existed at 1.25 to 2,50 pipe diameters downstream from the onset of sepa- 
ration. Their results agreed with those of Phaneuf and Netzer [19] in 
that the reattachment length was unaffected by the Reynolds number in the 
turbulent range. Also, they found that the peak Nusselt number became 
spread out into a broad zone for strong separations caused by large step 
heights, This transition from a point to a broad zone occurred between 
h/D = 0.l6 and h/D = 0.25, where D is the pipe diameter and h is the step 
height. It was also found that increasing the step height moved the 
peak point slightly downstream. 

Back and Roshke [20] and Roshke and Back [21] also investigated f.he 
reattachment length by visual observation in water flows. The diameter 
expansion ratio used was 2.6 and a conical contraction section just up- 
stream reduced the boundary layer thickness. Small 0.76 -nm diameter 
holes were located on the larger diameter tube one step height apart. 

Reattachment locations were determined by slowly metering dye through 
the holes and observing whether it moved upstream or downstream. jj 

Phaneuf and Netzer [19] measured wall static pressure as well as 
mean axial velocities behind abrupt steps where the expansion ratios 
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were 0.303 and 0.195* They found that the reattachment zone spreads out 
with both Increasing Reynolds number (for a given step height) and in- 
creasing step height providing the Reynolds number is in the turbulence 
transition range. It was further concluded that the point of maximum 
heat transfer does not coincide with the zone of reattachment and further- 
more it is dependent upon the inlet mass flux. 

Three different step heights of 0 . 6 , 1 . 6 , and 2.0 were studied by 
Ha Minh and Chassaing [ 22 ]. They measured the time-mean axial veloci- 
ties and the turbulence characteristics of the flow. Agreement between 
hot-wire and pitot probe measurements were found to exist for the time- 
mean velocities. The Reynolds stresses were successfully measured using 
the rotating, inclined, single hot-wire technique, and their results com- 
pared well with previous studies [17]. 

Laser doppler anemometry was used by Mocn and Rudinger [23] to study 
the behavior of recirculating flows. Even though an L.D.A. system was 
used, only mean axial velocities were measured, particularly those found 
in regions of recirculation. The L.D.A. system with the aid of a Bragg 
cell could detect the velocity reversals in the recirculation region. 

Johnson and Bennett [24] developed a measuring technique to deter- 
mine the mean and fluctuating velocity and concentration distributions 
along with their cross correlations by using a laser doppler veiocimeter. 
This method was developed using coaxial jets in sudden expansions with 
D/d = 2.07. A laser- induced flourescence technique was used to measure 
the concentration of a flourescent trace material. 

Janjua [II] measured the time-mean velocities and complete Reynolds 
stress tensor for an abrupt and gradual expanding confined jet. The six- 
orientation single hot-wire technique was employed in a confined jet of 
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expansion ratio 2. The results obtained using this technique were found 
to be in good agreement with those found by Chaturvedi [17]. 

2.2 Free Swirling Jets 

Swirl can be imparted upon a flow by many different means. Rose 
[25] used a very simple method of producing swirl by rotating a pipe at 
9500 rpm so as to provide an approximately fully-developed turbulent 
flow in solid body rotation. By this means it was only possible to ob- 
tain a relatively weak degree of swirl. Because of the low swirl number 
involved, Rose was capable of using a single hot-wire and an X-meter in 
the flowfield to determine the turbulence intensity components and the 
direction of the flow velocities, respectively. He found a free swirl- 
ing jet, in contrast to a nonswirling jet, that the jet spreads at a 
larger angle, entrains fluid more rapidly, and consequently displays a 
more rapid reduction of mean velocity and turbulence intensity. 

Chigier and Beer [26] generated swirl by introducing varying propor- 
tions of air through tangential ports while the remainder of the air was 
introduced axially. General velocity and static pressure distributions, 
form and shape of the central toroidal recirculation zone (CTRZ) were ob- 
tained for differing degrees of weak swirl and varying discharge nozzle 
conf 1 gurat i ons . 

Using a three-dimensional probe, Kerr and Fraser [27] found that the 
jet mass rate, angle of spread, and the reciprocal of the axial velocity 
all vary linearly with the swirl number. They proposed that the swirl 
number S be defined as 
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t 




where T is the nozzle torque, G is the jet thrust, and D g is the effec- 
tive nozzle u.ameter. 

Cbigier and Chervinsky [28] carried out an experimental investiga- 
tion on a series of swirling air jets. The swirl generator was composed 
of axial and tangential air inlets. Measurements of time-mean values of 
velocity and static pressure were made with a five-hole spherical impact 
probe. With this instrument they found that there was a progressive de- 
crease in the length of the potential core for an increasing value of 
degree of swirl. In addition, they noticed that for weak and moderate 
swirl the velocity profiles are effectively similar from an axial dis- 
tance of four nozzle diameters downstream of the jet exit. For stronger 
swirling jets similarity is not achieved until a distance of ten jet dia- 
meters has been reached. 

Using a rotating pipe similar to Rose [25], Pratte and Keffer [29] 
measured the time-mean properties, turbulence, and distribution of kine- 
tic energy via means of fine pitot probes and single hot-wire anemometry. 
Again this was possible because of the weak component of swirl velocity. 

The most common and effective method of producing swirling flows is 
via use of vane swirlers. Mathur and MacCallum [30] conducted an evalua- 
tion of swirler efficiency for both annular (with hub) and hubless swirl- 
ers. Extensive presentation of velocity profiles and recirculation zone 
geometry was given as a function of the degree of swirl. They found that 
for increasing values of swirl a central toroidal recirculation (CTRZ) is 
formed. This is caused when the adverse pressure gradient along the jet 
axis cannot be further overcome by the kinetic energy of the fluid flow- 
ing in the axial direction. They also concluded that tan $ was a good 
measure of swirl for jets issuing from vane swirlers. 


Syred et al. [31] concentrated on turbulence measurements in strong- 
ly swirling flows via a single hot-wire six-orientation technique inside 
the central recirculation zone. Local turbulence intensities were found 
to be extremely high in and near the CTRZ. Measurements of al 1 six tur- 
bulent stress components showed strong variations of absolute kinetic 
energy levels and strong nonisotropy of the stresses and associated tur- 
bulent viscosity. 

Claypole and Syred [32] studied the effect of combustion on the 
CTRZ using a dual beam laser. They found that combustion considerably 
lowers the level of recirculation even at high levels of swirl. A broad- 
er coverage of the measuring techniques and results for free swirling 
flows can be obtained from Reference [33]- 

2.3 Confined Swirling Flows 

Very little in the open literature has been written concerning tur- 
bulence measurements in confined swirling jet flows. The measurements 
that have been reported have been made with highly expensive nonintru- 
sive instrumentation. However, several researchers have reported on 
time-mean veloci'y measurements within confined flows with swirl. 

Mathur and MacCallum [34] studied the characteristics of a swirling 
jet issuing into a model of square cross-section. They used an isother- 
mal air model to obtain time-mean velocity profiles with a five-hole 
pitot probe and a water model to obtain qualitative data. Mathur and 
MacCallum found that the initial rates of spread of enclosed jets is 
more rapid than for free swirling. Also, that the central toroidal re- 
circulation zone is much stronger in enclosed jets than free jets. 
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Rhode et al . [8] measured the time-mean velocity properties again 
using a five-hole pitot probe inside a swirling confined jet. By using 
interchangeable expansion blocks with varying degrees of expansion angle 
a with increasing swirl strength, they found that the influence of a is 
negligible for the CRZ. But the effect of the swirl strength was to 
shorten the length of the CRZ and generate the existence of the CTRZ. 

The length of the CTRZ was then found to increase for increasing swirl 
strength. They also noted that on the jet centerline downstream of the 
CTRZ a precessing vortex core (PVC) existed. This PVC then continued to 
the exit of the test section. However, this study only included weak 
and moderate swirl strengths. This work was extended by Yoon and Li 1 ley 
[10] to include higher swirl strengths and to study the effect of differ 
ent contraction nozzles placed at various axial locations in the test 
chamber. They found that a strong contraction nozzle of area ratio k 
greatly affected the size and shape of the CTRZ and PVC. 

Beltagui and MacCallum [35] completed a wide range of experiments 
on annular and hubless swirlers with and without combustion. The expan- 
sion ratios they used were D/d = 2.5 and 5-0 with varying degrees of 
swirl. Mean flow and static pressure measurements were made with a water 
cooled five-hole spherical pitot probe. Their investigation concluded 
that the maximum diameter of the CTRZ is primarily a function of the 
chamber diameter (approximately 0.65 chamber diameters) and only slight- 
ly alters for varying degrees of swirl, by combustion or by D/d ratio. 
They also noticed that the presence of a central hub in the swirler had 
little effect on the flow patterns. 

Using an expansion ratio of D/d = 2 and a five-hole pitot probe, 
time-mean velocity measurements were made in a swirling isothermal and 


12 


reacting confined jet by Syred and Dahman [36]. They studied the effect 
of introducing bluff bodies into the jet on the size and shape of the 
central recirculation zone. Also, they found that by increasing the fur- 
nace exit diameter, only a small effect was noticed upon the aerodyna- 
mics of the system. They also reported that for furnaces requiring low 
calorific value gas, only medium swirl produced the desired recircula- 
tion zone. 

Baker et a!, [7] and Hutchinson et al. [6] used laser doppler anemo- 
metry to measure the values of mean velocity and turbulence intensity in 
the recirculation region of a swirling reacting flame. Methane was ig- 
nited at the burner and the flow then seeded with TiO^ particles to en- 
hance the light scattering. As flow reversals are found in swirling 
flames, a frequency-shifting technique was used to solve the directional 
ambiguity of the doppler signal. Both groups of researchers found high 
levels of turbulence near the exit of the jet for even a moderate swirl 
number. 

Laser velocimeter measurements of a confined turbulent diffusion 
flame burner was made by Owen [5]. Mean and rms turbulent velocity lev- 
els were measured for axial and tangential components of. the reacting 
flowfieid. He measured four diffe-rent swirl strengths near the exit of 
the burner and found there was no significant variation on the time-mean 
profiles for increasing swirl. He also found that the turbulence is 
high with significant deviation from isotropy over the initial mixing re- 
gion of the jet. In this area large scale fluctuations decrease for in- 
crease in swirl strength. A bragg cell was used for directional sensi- 
tivity of the flow. 
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Habib and Whitelaw [37) measured the time-mean velocities and direc- 
tional turbulence intensities in a swirling confined coaxial jet. Two 
measuring techniques were used to measure the turbulence, hot-wire anemo- 
metry, and L.D.A. The two methods were in close agreement even in the 
corner recirculation zone. These researchers only used low values of 
swirl strength; therefore, no central recirculation zone was formed. 

Using a five-hole pitot probe and a single normal hot-wire, Vu and 
Gouldin [ 38 ] measured the time-mean and turbulent properties of a confin- 
ed isothermal coaxial jet under co-swirling and counter-swirling condi- 
tions. They found that levels of high turbulent fluctuations and large 
dissipation rates characterize the central flow region for both co- and 
counter-swirl. They suggested that with chemical reaction the internal 
structure of the recirculation zone is of much less importance to the 
combustion process than the turbulent transport properties in the vicin- 
ity of the CTRZ and the interjet shear layer. However, they reported 
that no central recirculation zone was formed for the co-swirling case. 

Gouldin et ai. [ 39] continued this work with special emphasis on 
•co-swirling flows with combustion. They found that with the presence of 
chemical reaction a CTRZ was formed. This seems surprising as mass con- 
servation requires axial flow acceleration with combustion and therefore 
a favorable pressure gradient must exist over a portion of the test sec- 
tion. Since the formation of the CTRZ requires an adverse pressure gra- 
dient on the centerline, one would expect combustion to inhibit recircu- 
lation zone formation. However, the researchers showed by application 
of Bernoulli's equation that the pressure on the centerline increases 
wi th combustion. 
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No PVC was found by Brum and Samuel son [40] in dilute swirl combus- 
tors with or without combustion. They suggest that the PVC is an arti- 
fact of the modeling process and is not present in full scale combustors 
as dilution jets are present. The researchers found that if only 70 per- 
cent of the inlet air was swirled and the remainder used as dilution air, 
no PVC existed. 

The six-orientation single hot-wire technique was used by Janjua 
etal. [41] in an isothermal nonswi rl ing and swirl ing confined jet wi th 
an expansion ratio of 2. Measurements in the nonswirling jet compared 
favorably with those of previous researchers [8, 17]. The swirling jet 
results presented were quite limited but did show that around regions of 
recirculation large values of turbulence intensity and shear existed, 
even for moderate swirl. The researchers also performed an uncertainty 
analysis on the hot-wire technique. This involved changing certain in- 
put parameters to the data reduction and determining the percentage 
change in output parameters. They found that the most inaccurate output 
term was u 'w 1 . 

Sommer [42] measured time-mean and fluctuating properties of con- 
fined coaxial jets and found that a CTRZ did exist for the co-swirling 
case. He found that the recirculation zone consisted of a one-cel led 
toroidal vortex characterized by low internal swirl velocity, high tur- 
bulence intensity, and large turbulent dissipation rates. 

The previous experimental work presented here gives a good indica- 
tion of the levels of the properties that are quantified in this study. 
However, before any direct comparison can be made between these reported 
results and the present measured data, the flows must be similar. This 
must involve not only geometric similarity but also, and more importantly, 
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the inlet flow conditions must be dynamically similar. This often in- 
volves the use of a similar inlet contoured contraction nozzle and the 
same number and shape of swirl vane blades. 
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CHAPTER ill 


MEASUREMENT TECHNIQUES AND ANALYSIS 

Turbulence measurements in a complex flowfield have always been a 
complicated problem encountered by engineers. In the past turbulence 
phenomena have been discussed by various authors in detail and various 
methods of turbulence measurements have been suggested [ 3 3 » A3, AA] • One 
of the most widely-used instruments to obtain turbulence quantities is 
the hot-wire anemometer, the most common of which is the single hot-wire. 
When used at a single orientation and in a two-dimensional flow with a 
dominant flow direction, a single hot-wire can measure the streamwise 
components of the time-mean velocity and the root-mean-square velocity 
fluctuations at a particular location in the flowfield. A two-wire 
probe can be used to determine the time-mean velocities, streamwise and 
cross stream turbulence intensities, and the cross correlation between 
the two components of the velocity fluctuations [A3, AS]. 

Hot-wire measurements in a complex three-dimensional flowfield are 
considerably more difficult than in one- or two-dimensional flowfields 
in which the mean flow is predominantly in one direction. To measure 
the three velocities and their corresponding fluctuating components in a 
three-dimensional flowfield such as encountered in combustor simulations, 
a triple hot-wire technique is typically used. 

The three-wire probe technique permits the necessary simultaneous 
measurements from which three instantaneous velocity components can be 


17 


determined. The appropriate signal processing can produce estimates of 
mean velocity components and the Reynolds stress tensor. However, the 
three-wire probe technique is very complex, requiring a multi-dimens ional 
probe drive to align the probe in the mean flow direction. Also, sophis- 
ticated signal processing electronics is required to handle the three in- 
stantaneous hot-wire voltages. Finally, because of the fact that three 
hot-wires are involved, the dimensions of the instrument tend to be large, 
causing poor spatial resolution. 

In recent years the use of single-wire multiple-orientation tech- 
niques is rapidly growing in popularity. This is because of their abil- 
ity to measure all the time-mean velocities and all the components of 
normal and shear turbulent stresses in complex flowfields. The method 
is very cost effective, requiring only standard hot-wire electronics. 
Assumptions regarding the nature of the turbulence are that it is sta- 
tionary and that it follows a normal probability distribution, it is 
for these reasons that makes the six-orientation, single hot-wire tech- 
nique suitable for swirling confined jet flows. 

3.1 Single-Wire Six-Orientation Technique 

Multi-orientation of a single hot-wire is a novel way to measure 
the three components of a velocity vector and their fluctuating compo- 
nents, A method devised by Dvorak and Syred [46] uses a single normal 
hot-w i re or iented at three different positions such that the center one 
is separated by 45 degrees from the other two. The velocity vector at a 
location is related to the three orthogonal components using pitch and 
yaw factors as defined by Jorgensen [47 j . The data are obtained in the 
form of mean and root-mean-square voltages at each orientation. However, 
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the measurements done with a single wire do not supply all the informa- 
tion needed to obtain the turbulence quantities. Therefore, in addition 
to a single wire, Dvorak and Syred used a cross-wire probe to obtain the 
covariances between the voltages obtained at adjacent hot-wire orienta- 
t i ons . 

King t^8] modified the technique developed by Dvorak and Syred. His 
method calls for a normal hot-wire to be oriented through six different 
positions, each orientation separated by 30 degrees from the adjacent 
one. Mean and root-mean-square voltages are measured at each orienta- 
tion. The data reduction is performed using some assumptions regarding 
the statistical nature of turbulence, making it possible to solve for 
the three time-mean velocities, the three normal turbulent stresses, and 
the three turbulent shear stresses. 

The six-orientation, hot-wire technique requires a single, straight, 
hot-wire to be calibrated for three different flow directions in order 
to determine the directional sensitivity of such a probe. The three 
directions and three typical calibration curves are shown in Figures k 
and 5- Each of the calibration curves follows a second order, least 
square fit of the form: 


-2 . A _ - 1/2 . r ~ 
E. = A. + B.u. + C.u. 

1 1 11 11 


( 3.0 


which is an extension of the commonly used King's law. In this equation, 
A., B., and C. are calibration constants, and u. can take on a value of 
u, v, and w for the three calibration curves, respectively. 

When the wire is placed in a three-dimensional flowfield, the effec- 
tive cooling velocity experienced by the hot-wire is: 

2 = v +Gu +Kw \3*2J 
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where G and K are the pitch and yaw factors defined by Jorgensen [^7] 
to be: 


and 


v(u and w = 0) 
u(v and w ■> 0) 


„ v(u and w = 0) 

\\ a ii ■ ■ 

w(u and v * 0) 


(3.3a) 


(3.3b) 


which are evaluated from the three calibration curves (Figure 5) for a 
constant value of E . Equation (3-3) shows that the pitch and yaw fac- 
tors are calculated with the v component i = 2 in Equation (3-0 of the 
effective cooling velocity as the reference. Therefore, the calibration 
constants used in Equation (3.1) are the coefficients in the E versus v 
calibration of Figure 5i i.e., in a general flowfield: 

2 l /2 

E = A 2 + B 2 Z W£ + C 2 Z 

with Z as given in Equation (3-2) above. 

Figure 6 shows the pitch and yaw factors as a function of hot-wire 
voltage determined from the calibration curve of Figure 5. Both factors 
vary with hot-wire voltage, but the yaw factor is far more sensitive. 

The sensitivity analysis discussed in the next section demonstrates that 
uncertainties associated with the varying pitch and yaw factors do not 
seriously affect the accuracy of the estimated flow quantities. 

To carry out measurements in the confined jet flowfield, the wire 
is aligned in the flow in such a way that in the first orientation the 
wire is normal to the flow in the axial direction and the probe coordi- 
nates coincide with the coordinates of the experimental facility. Thus, 
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the six equations for the instantaneous cooling velocities at the six 
orientations, as given by King [A8], are: 


,2 Z x J 2 x „2 2 . 

Zj ° v +Gu +Kw (3.4a) 

Z 2 = v 2 + G 2 (u cos 30° + w sin 30°) 2 

+ K 2 (w cos 30° - u sin 30°) 2 (3.Ab) 

Z 2 » v 2 + G 2 (u cos 60° + w sin 60°) 2 

+ K 2 (w cos 60° - u sin 60°) 2 (3.4c) 

Z^=v +Gw + K u (3.4d) 

Z 2 « v 2 + G 2 (w sin 120° + u cos I20 e ) 2 

+ K 2 (u sin 120° - w cos l 20°) 2 (3.4e) 

Z^ = v 2 + G 2 (w sin 150° + u cos 1 50 ° ) 2 

+ K 2 (u sin 150° - w cos 150°) 2 (3.i*f) 


Solving simultaneously any three adjacent equations provide expressions 
for the instantaneous values of the three velocity components, u, w, and 
v in terms of the equivalent cooling velocities (Zj, Z^, and Zy for ex- 
ample, when the first three equations are chosen). It is then possible 
to obtain the three time-mean velocity components and the three normal 
turbulent stresses of the Reynolds Stress tensor, in a manner described 
in Reference [1*1]. 

To determine the cross-correlations of the turbulent fluctuations, 
a correlation coefficient has to be assumed. King [48] argued that if 
two wires are separated by an angle of 30 degrees, the fluctuating sig- 
nals from the wires at the two locations would be such that their con- 
tribution to the cooling of the wire would be related by the cosine of 
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the angle between the wires. This assumption leads to the following two 
values of the correlation coefficients: 

Y ZpZQ ° cos ° 0,887 (3.5a) 

YZqZr = cos 30 = 0.867 (3-5b) 

To relate Y^p^Q with YzqZr an8 ^ZrZr’ King [48] Introduced the following 
rel at i onsh ip: 

TZpZR ~ nY ZpZQ " Yzqz r (3*6) 

where n is given a value of 0.8. A sensitivity analysis presented in 
Chapter VI demonstrates that there is no significant error magnification 
In the data reduction due to the correlation terms. 

To obtain measurements of the correlation coefficients in a swirl- 
ing flowfield, a triple hot-wire and a time correlator would have to be 
used. Dvorak and Syred C^6 ] used a double-wire in a two-dimensional 
flow and found that the above assumptions are quite good. 


3.2 Eddy Dissipation Rate Measurements 

Studies of the en**>*gy processes in shear flows have all suffered 
from uncertainty in determining the rate of turbulent dissipation. Hinze 
[1*9] shows that thi s process can be expressed in cartesian coordinates by 



(3-7) 


(where u! is the component of the fluctuating velocity in the x. direc- 
tion). Assuming incompressible, homogeneous, small scale turbulence, 
Hinze reduces this equation to the form: 


22 



Lawn [50] suggests that complete Isotropy is not required for Equa- 
tion (3-8), but the condition that the transfer of energy from one velo- 
city comnonent to another be relatively small suggests that the normal 
stresses due to these components be approximately equal. This condition 
is defined as ''second class local isotropy" by Bradshaw [1*3]. The assump- 
tion does not neglect the existence of shear, and in fact includes an in- 
herent nonisotropy by giving special emphasis or; the mean flow direction. 
This reasoning states that all the turbulent energy is initially produc- 
ed in the component of the flow direction and is then distributed among 
other components. This turbulent energy could be distributed by means 
of the pressure fluctuations [ 50 ]. 

if Bradshaw's definition applies, and assuming that all the turbu- 
lent dissipation occurs in the mean flow direction, then Equation (3-8) 
reduces to: 



(3-9) 


where s is i n the direction of the local time-mean velocity. 

Invokinq Taylor's hypothesis (S/8t 3/9s) , that the fluctuations 

at any point in the flowfield are caused by the whole turbulent flowfieid 
passing that point with a constant velocity V, it follows that 


E 




(3-10) 
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Heskestad [51] has suggested that the use of Taylor's hypothesis can 
possibly overestimate the dissipation rate and therefore introduced a 
modification to the convective velocity: 


S 2 


(— ) „ 

V 3t 1 


*■8 x ’ 


(1 + 


,2 


+2 

V 


+ 2 


+ v/ ' 
? 2 


(3.H) 


However, Antonia et al. [52] did not apply this modification to mea- 
surements in the fully developed region of axisymmetric free jet at axial 
locations greater than x/D = 40. They found that there was a linear de- 
crease in dissipation rate for increasing axial location. 

Laufer [53] used Equation (3-7) to determine the eddy dissipation 
rate. He measured five of the nine contributions to the sum and assumed 
that isotropic relations may be used to derive the remainder. His mea- 
surements were performed in a fully-developed turbulent pipe flow, a re- 
latively easy flowfield. However, his results contained several incon- 
sistencies mainly because of the errors incurred while measuring five 
different derivatives. Lawn [50] also measured the energy dissipation 
in turbulent pipe flow. He concluded from his own and Laufer c s results 

[53] that Equation (3.3), which assumes complete isotropy, will not be 

4 

grossly in error if the pipe Reynolds number is greater than 9 x 10 . Any 
slight errors that do exist could be attributed to poor high frequency 
response of hot-wires and to finite wire length. Isotropy is approached 

as the Reynolds number increases and is a good approximation at Re = 9 x 

4 - L 

10 . Lawn [50] used Equation (3-8) with a Reynolds number of 4x10 

with a single hot-wire and then Equation (3-7) with a double wire; he 

found that the results did not show a great deal of scatter between the 
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two methods, indicating that a single hot-wire could be used at the lower 
Reynolds number. 

Habib and Whit-elaw [54] measured the dissipation rate in a confined 
coaxial jet without swirl. They used the isotropic relationship of Equa- 
tion (3.8) and compared their results with those obtained with the TEACH 
computer code [55l- The researchers found that the calculated results 
were a factor of two higher than the measured results. However, the mea- 
surements were only of a preliminary nature and the authors concluded 
that neither the precision nor the detail of the measurements were suffi- 
cient for satisfactory appraisal. 

In confined coaxial swirling flows, Vu and Gouldin [38] used Equa- 
tion (3-10) to determine the dissipation rate. From this they were able 
tc calculate the length scale of dissipation from the well-known k-e re- 
lationship. In addition, estimates of the Taylor microscale and the 
Kolmogorov length scales were obtained using isotropic relationships and 
the eddy dissipation rates. The fact that Vu and Gouldin [ 38 ] have no 
comparisons to make wi th their results in swirling recirculating flows 
attests to the fact that accurate measurements in such a flow are quite 
difficult. 

To avoid errors associated with taking many turbulent derivatives, 
the isotropic relationship of Equation (3.10) is used in this study to 
determine the eddy dissipation rate. The high Reynolds numbers employed 
in this project also gives another point in favor for us i ng the i sot rop i c 
assumpt ions . 

3.3 Contribution to Turbulence Modeling 

The practical need of the turbulence measurements in swirling flows 
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can be seen from studying the Reynolds equations. In solving the Rey- 
nolds equations for time-mean velocities, specification is required of 
the turbulent stresses, usually via the turbulent viscosity p^_. While 
the molecular viscosity is a real property of the fluid, existing whe- 
ther the fluid is in motion or at rest, the turbulent viscosity requires 
some flow of fluid to become effective and is thus not a property of the 
fluid but a parameter of the fluid motion. This parameter describes the 
behavior of the turbulent stresses in terms of the time-mean velocity 
gradients, implying the assumption that the turbulence transport is of 
the gradient type. Although simple models are availau-e for the 

specification of p , especially in boundary layer flows, for example, 
they are not useful for recirculating flows. 

Increasingly complex flowfields are governed by increasingly compli- 
cated equations for the time-mean properties, which contain increasingly 
more components of the turbulent (Reynolds) stress tensor. For example, 
the increasing number of components required are: 

1. Nonswirling boundary layer f low--rx-component only; 

2. Swirling boundary layer flow--rx and r0 components; 

3- .Nonswirling axisymmetric recirculating flow — rx, xx, rr, and 90 
componen ts ; 

k . Swirling axisymmetric recirculating flow — rx, r9, x0, xx, rr, 
and 00 components; 

with more complex flows requiring, again, all six components. 

In nonswirling axisymmetric recirculating flow, the turbulent (Rey- 
nolds) equations of conservation of axial and radial momentum may be 
solved (along with the conservation of mass [continuity] equation) so as 
to obtain t'me-ntean values of p, u, and v (time-mean pressure, axial and 
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radial velocities) throughout the flowfield. The equations for u and v 
contain only the following components of the turbulent (Reynolds) stress 
tensor jl: 


,2 

T = *pU’ 

XX 

= 2p tl7-f Pk 

(3.12a) 

~2 

T = -PV* 

rr 

- 2p t If - 1 pk 

(3.12b) 

~/2 

T ea = " pw 

= ' 2,J tTr J pk 

(3.12c) 

T = *PU' V 1 

rx 

,3u 3v\ 

(3-1 2d) 


which are usual' 1 y related by way of an isotropic turbulent viscosity m 
to the appropriate components of strain (velocity gradients) in the form 
also shown. Incompressibility has also been assumed and the terms con- 
taining dlv v^ have been omitted. Notice that the normal components also 
contain the term pk (which is often omitted in practice) in order 

that summation of the normal components and application of the continu- 

2 2 2 

i ty equation are consistent with the definition of k = (u 1 + v' +w’ )/2. 

If swirl is also present, coupled with the above is an equation for 
conservation of swirl momentum from which the time-mean swirl velocity w 
may be obtained throughout the flowfield. Additional terms also appear 
in the other momentum equations. From the turbulence simulation point 
of view, additional components of the turbulent (Reynolds) stress tensor 
appear in the swirl equation: 

T„, - -PUSF - „ £ ■ (3.Ua) 

t r0 - -pv'w 1 « ur (~ ) (3.13b) 
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These correlations are usually related by way of an isotropic turbulent 
viscosity p to the appropriate components of strain {velocity gradients) 
in the form shown on the extreme right hand of these equations. Thus 
all six different components appear in the governing equations. Notice 
that axisymmetry has reduced the complexity, not only of the governing 
equations themselves but also of the constitutive relations. 

If isotropic turbulent viscosity is not assumed, different values 
of u may be appropriate to different components of the turbulent (Rey- 
nolds) stress tensor j_. That is, a double suffix is added so as to indi- 
cate the component in question. Then by analogy with Prandtl -Schmi dt 
numbers, viscosity numbers may be defined so as to relate other viscosi- 
ties to the primary component of turbulent viscosity p = p . (This 
particular one is chosen, since it is the only component which remains 
in the governing equations of a nonswirling boundary layer flow, on which 
much previous research has taken place. Swirling boundary layers con- 
tain p in addition to p and only one viscosity ratio appears.) Thus 
one is interested in whether or not the following viscosity numbers are 
unity (isotropic) or not (nonisotropl c) : 


= p/u 

(3.14a) 

XX XX 


rr = ^rr 

(3- 1 4b) 

00 = U/a 00 

(3-1 Ac) 

rB = " /M re 

(3.1 Ad) 

x0 = p/p x0 

(3- I3e) 


There is some evidence that these values are not unity [31, 55], but gen 
era I recommendations are not available [56]. The remaining problem is 
how the main component of turbulent viscosity is to be specified by a 
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turbulence model. Or, rather, the more complete question of how are all 
the components of turbulent (Reynolds) stress to be specified, either 
directly (advanced differential or algebraic stress mode I i ng) or i nd i rect - 
ly (via the turbulent viscosity concept). 

A turbulence model is a set of differential and/or algebraic equa- 
tions for things like k, e, u ! uj » ulm^, m^, etc. which connect statisti- 
cal properties of turbulence (correlations) with each other and with 
terms appearing in the time-averaged equations of conservation of mass, 
momentum, energy, and chemical species. Auxiliary relations are like 


u!ul = -const * 3u,/3x. 

' J i J 

— 1/2 

R £ = -const (m^ ) pek 


(3- ! 5a) 
(3.15b) 


Models are classified according to the turbulent flux hypothesis 
(whether or not turbulent exchange coefficients are introduced) and the 
number of extra differential equations to be solved. If introduced, ex- 
change coefficients have generally been assumed isotropic until recently, 
even in flows with swirl, but recent experimental, inverse and predic- 
tion works have disputed this for swirling flows. Briefly the choice 
available is: 

2 1 /2 

1. Prandtl mixing length p = pi ( ZA : a ) 

■ 2-1 

2. Energy- 1 ength k e ; 

3- Differential stress modeling Dr /Dt = P + D +R + c : 

rx rx rx rx rx’ 

k . A1 geb ra i c st ress mode ling x - f (other t 1 s , k , e. A) . 

rx 

Here A is the time-mean flow rate of strain tensor, k and c are turbu- 
lent kinetic energy and dissipation rate, and P, D, R, and e stand for 
production, diffusion, redistribution, and dissipation In the turbulent 
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stress equation. The first two choices are examples of theories of the 
exchange coefficient type; the second two are of direct stress specifica- 
tion type. 

Currently two-equation energy-length models are to be recommended 
for application in practical engineering situations; in particular the 
k-c model where i .= k ’ /£ and 9 . is the macrolength scale of turbulence. 
Models of this type have been successfully modified and applied to swirl- 
ing flows, even with nonisotropic assumptions [58, 59], Recent relevant 
work is discussed elsewhere [60]. Numerical and analytical inverse 

solution of the turbulent swirl flow boundary layer equations (which 
allow turbulence model development directly from time-mean experimental 
data) have been useful for this purpose. Details of the general turbu- 
lence model development in both nonreacting and reacting flows are avail- 
able via recent conferences [61, 62]. 

The advantage of a turbulent k i neti c energy TKE model over a simple 
algebraic "local 1 ' model is that solution of a PDE for k allows upstream 
"historical"’ influence to assert itself on the subsequent flowfield. 
Application of several such models to flowfields of interest in. sudden 
expansion burners have been reviewed elsewhere [63], For high speed tur- 
bulent reacting flows of specific interest to ramjet dump combustors, 
two models show significant promise: the Harsha one-equation model [64] 

and the Rod! two-equation k-e model [57]. Applications to reacting jet 
flows shows that both appear to be suitable [ 65 ] - More advanced turbu- 
lence models, such, as those based upon the Reynol ds-st ress modeling ap- 
proach, are not yet fully developed to warrant their use in recirculat- 
ing flowfield problems as encountered in gas-turbine combustors. (n 
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addition, such an approach will appreciably increase the computation ef- 
fort. 

Measurements of the complete Reynolds stress tensor and the time- 
mean velocity gradients can lead to knowledge of the turbulent eddy vis- 
cosities via use of Equations (3*12) and (3-13)* Determination of the 
rate of eddy dissipation in conjunction with the turbulent viscosity can 
be used to introduce a numerical value to the constant term in Equa- 
tion (3*15). Typically, is given the constant value of 0.09 in swirl- 
ing flows. However, it has been shown that this value is unrealistic 
and is more likely to be a function of the velocity gradients [66]. The 
measured data from this study will be used to obtain a relationship for 
the values of [67] . 


CHAPTER IV 


DIRECTIONAL SENSITIVITY ANALYSIS 

The analysis is performed at any specific flowfield location by ini- 
tially placing the probe in a free jet such that the coordinate system 
of the probe coincides with the coordinate system of the jet, as shown 
in part (a) of Figure 7- Measurements are then taken by rotating the 
probe in the manner of the technique just described. To simulate the ef- 
fect of the flow shifting its dominant flow direction, the probe is ro- 
tated by 9 degree about its z-axis, as shown in part (b) of the figure. 
This rotation causes a misalignment between the probe coordinate system 
and the facility coordinates. This discrepancy can be accounted for by 
use of the Eulerian matrices described in section ^.1. In this configur- 
ation, the measured time-mean values, normal and shear stresses are in a 
coordinate system oblique to the jet coordinate system. However, they 
can be transformed back to the facility coordinate system [68], Notice 
that, the correct directional sense of the rotation must be followed so 
chat standard coordinate transformations may be used on the probe data 
so as to obtain facility coordinate data. Results shown later in Chap- 
ter VI have been obtained in this manner. 

To examine the directional sensitivity of the wire further, the 
probe was subsequently rotated about its new x-axis, thereby forming a 
compound angle between probe and the dominant flow velocity, as also 
shown in Figure 7, part (c) . Again, the time-mean velocities and 
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Reynolds stress tensor can be deduced in terms of the jet coordinate sys- 
tem by the method shown In section A. I, and the results and their accu- 
racy are also discussed in Chapter VI. 

Specifically, the directional sensitivity of the techn ique i s assess- 
ed at five flowfield situations for self-consistency using the following 


five configurations: 

Case 1 9=0° $ = 0° 
Case 2 0 « -45° <j> = 0° 
Case 3 e = -45° <f> « -4 5° 
Case 4 0 = -90° <p = 0° 
Case 5 0 « -90° <j> - -90° 


The above five probe/flow configurations are used at each of five 
representative situations in a free ax i symmetr i c nonswirling jet at x/d 
= 0 (laminar in potential core region), 3 and 10 (turbulent in shear 
layer region) and in a free ax i symmetr ic swirling jet at a location in a 
region of strong shear just downstream of the exit from a variable-angle 
vane swirler with swirl vane angles of 45 and 70 degrees, representing 
moderate and strong swiri cases. In these swirling jet cases, the probe 
was located just downstream of the swirler exit, outside of any regions 
of recirculation: this part of the flow was chosen as it is in an area 

of rapid acceleration and is unlikely to contain any instantaneous flow 
reversals which might cause erroneous readings. Specifically, the fol- 
lowing five situations are used: 

Situation A x/d = 0 Nonswirling laminar region 

Situation 8 x/d = 3 Nonswirling turbulent region 

Situation C x/d = 10' Nonswirling turbulent region 
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i 


1 

j 



Situation 0 


x/d = 0 


Swirling turbulent region with 
swirl vane angle 45 degrees 
Situation E x/d = 0 Swirling turbulent region with 

swirl vane angle 70 degrees 

These locations are illustrated in Figures 8 and 9. It will be seen later 
that the sensitivity analysis assures users that knowledge of local con- 
figuration of probe versus flow direction is not required a priori, and 
useful results are forthcoming and relatively insensitive to specific 
conf igurat ions. 

4,1 Coordinate Transformation 
4.1.1 Rotational Matrices 


To investigate the shifting of the dominant flow direction, the 
orobe holder versus local time-mean flow vector configuration is varied 
through five different cases of interest. To relate probe coordinate 
data back to the facility coordinate system, use is made of Eulerian ro- 
tational matrices. These make it possible to relate all time-mean velo- 
cities and the full Reynolds stress tensor back to the facility coordi- 
nate system after any axis rotation. 

Any (x,y . z) -cartes i an coordinate axes may be rotated about the x, 
y, or z axis, respectively, by an angle 6, with corresponding coordinate 
transformation matrices: 


x9 


0 0 

cos0 -sinfl 

sine cos0 
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COS0 

0 

s 1 n0 

0 

1 

0 

-sin0 

0 

cos8_ 

cos 9 

-sinQ 

0 ‘ 

si n0 

COS0 

0 

0 

0 

1 



For example, a rotation about the z A axis by angle 0 results in old (x, 
y, z) coordinates of a point being related to its new (X, Y, Z) coordi 
nates via: 

[x y zl T = R zQ [X Y Z] T 

Similarly, velocity components {u, v, w) are related to (U, V, W) via 
[u v w] T = R z0 [U V W] T 

in the notation of Figure 1 with velocity components 


u, v, w 

i n 

x, 

r, 0 

facility coordinates 

u, v, w 

i n 

x, 

y. z 

probe Case 1 coordinates 

U, V, w 

i n 

x, 

|N ^ 

Y, z 

probe Case 2 and 4 coordinates 

U, V, w 

i n 

x, 

y, z 

probe Case 3 and 5 coordinates 


the following relationships prevail. 

4.1.2 Case 1 

The facility and probe Case 1 coordinates are coincident and 
[u v w] T =« [u v w] T 
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^.1.3 Cases 2 and b 


In this case, a rotation i)f 0 degree is applied about the z-axis, 
resulting in 

r iT 

[u v w] = R Q [u v w] 

Z 0 

which leads to 


u 


u cos 0 - v s i n 0 

V 


u s i n 0 + v cos 0 




w. 


w 

U 


From this the directional time-mean velocities in facility coordi- 
nates can easily be inferred. The angle terms in the relationships are, 
in fact, the directional cosines [69], which can be used to determine 
the stress transformations in terms of facility coordinates. The direc- 
tional cosines between any new coordinate axes, for example, x, y, and 
z, and the original coordinate axes x, y, and z can be conveniently tabu- 
1 a ted as fol lows : 


! 

L_X_ 

y 

z 

r! 

*1 

m i 

n l 

Y 

*2 

m 2 

rvi 

C 

z 

? '3 

m3 

n 3 


These directional cosines (£., m. , n. for i = 1, 2, and 3) can be used 
in the general three-dimensional stress equations governing the trans- 
formation of coordinates via [68]. 


2 2 2 

x = + m,T-~ + n.i^- + 2JL.m,T"* + 2m.n.t-- + 2n. H.x-'- 

xx 1 xx 1 yy 1 zz 1 1 xy 1 1 yz 1 1 zx 

x = as above with subscript 1 replaced by 2 
YY 
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t z2 = as above with subscript 1 replaced by 3 

T xy * ‘iVii * m l m 2 T yy * "iV^ * U \ m 2 * 

+ (m,n 2 + Vl )T y£ + ( V2 + VP T zx 

T yz = as above with subscripts (I, 2) replaced by (2, 3) 

t = as above with subscripts (1, 2) replaced by (3, 1) 
Z X 


These relationships lead to evaluation of all the stress components in 
facil ity coordinates, from corresponding data originally reduced i n probe 
coordinates. 


4.1.4 Cases 3 and 5 


To go from Cases 2 and 4 to Cases 3 and 5, a rotation of <j> degrees 
is applied about the ne\, x-axis, resulting in 

[u v w] T = R £ fu v w] T 

and 

fu v w] T = R z0 R^ [u v w] T 

which is 


w 


cos 0 - sinO cosifi + sine s i n4> 
sinO + cosO cos (f> - cos 9 sim}> 
0 s i n cos<l> 


*1 

■ 


U 


A 

V 


s 

L W J 


Again, this relationship permits time-mean velocity components to be con 
verted back into facility coordinates. The coefficient matrix defined 
the directional cosines relating the two coordinate systems; these can 
be used for normal and shear stress transformations in the manner de- 
scribed in section 4.1.3. 
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CHAPTER V 


EXPERIMENTAL FACILITY 
5.1 Calibration Equipment 

A small axisymmetric free jet is used for the static calibration of 

the hot-wire. The calibration jet facility, shown in Fi gure 10 , cons i s ts 

of a seamless contoured nozzle geometrically similar to that of the flow- 

field facility. A settling chamber and turbulence management section 

consisting of packed straws is just upstream of the contoured nozzle. 

The nozzle itself has a throat diameter of 3^ mm and is capable of pro- 

c 

ducing Reynolds numbers up to 6x10 (based on throat diameter). The 
air is delivered by a thermally stabilized compressed air generator and 
produces air at ±0.5°C to that of tbo wind tunnel. The air is control- 
led to the nozzle by means of a diaphragm valve and a Fisher and Porter 
mode! I0A1735A rotameter. This arrangement allows 30 calibration points 
to be read over a wide range of known velocities (found by independent 
pitot probe measurements). 

It is necessary to calibrate the hot-wire in a flowfield of known 
properties. The one-dimensional low turbulence level (essentially lami- 
nar) potential core of the free jet was utilized for this calibration. 

The hot-wire is calibrated for three different flow direction, u, v, and 
w, as shewn, in Figure 5. Each of the three calibration curves are ob- 
tained with zero velocity in the other two directions. The calibration 
curves demonstrate that the hot-wire is most efficiently cooled when the 
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flow Is in the direction of the u component, whereas the wire is most 
inefficiently cooled when the flow is in the direction of the w compo- 
nent. The three orientations of the probe were made possible by use of 
a rotary table and DISA 55H151 and 55H153 mounting tubes. 

The familiar King's law relationship 

E 2 = A + BZ n (5.1) 

which is employed as a calibration curve in most conventional hot-wire 
anemometry applications has considerable limitations. Perry and Morri son 
[69] have shown that this relationship is only satisfactory for mean 
velocities over a limited range. These researchers illustrated by dyna- 
mic cal ibrat ion that this form of relationship is only sui tabl e for flows 
with low turbulence levels. Typically, linearizers are used which apply 
the relationship of Equation (..I) to the hot-wire output signal. The 
values of A, B, and n are variable to suit individual wires. The limit- 
ed range of application of such an expression and its restriction to low 
turbulence levels makes the use of linearizers very undesirable in high- 
ly turbulent flows since it could lead to considerable errors. 

Dvorak and Syred [46] found that the expression 

2 1 /2 

E « A + BZ ' + CZ (5.2) 

gave an excellent fit to obtained calibration data with a very small 
standard deviation being obtained over a large velocity range (1 to 60 
m/s). Such a curve fit is a more realistic approach since the conven- 
tional relationship of Equation (5.1) is based on the steady flow heat 
transfer relationship and thus may well not provide an accurate descrip- 
tion of the behavior of individual wires. A relationship of the form of 
Equation (5.2) is used to fit a curve fit to the whole calibration curve 
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using a lease squares curve fitting technique. This curve thus allows 
the instantaneous velocity to wander over the entire range from zero up 
to the maximum value seen in the flow without distortion. 

A hot-wire does not record a true zero velocity voltage while under 
calibration because of the thermal convection velocity induced by the 
bouyancy arising from the wire's own heat. Thus the zero velocity vol- 
tage obtained during calibration is spurious and is not included in the 
points used to obtain the curve fit. 

A further point in favor of Equation (5-2) is made by Dvorak and 
Syred [46 3 - They showed that a curve fit of this form from static cali- 
bration data fitted very well the dynamic calibration data of Perry and 
Morrison [69] and thus stated that dynamic calibration is probably un- 
necessary. This is to be expected, as the diameter of the hot-wire is 
very small (usually 5 pm); therefore, the flow in the vicinity of the 
wire may be considered laminar. Thus the hot-wire sees a time-varying 
laminar flow in the turbulent flow whereas in the calibration rig it 
sees a time-steady laminar flow. This also explains the good approxima- 
tion achieved by the expression in Equation (5-1) despite its having 
been based on laminar heat transfer analysis. 

5.2 Confined Jet Facility 

The swirling confined jet test facility is shown schematically in 
Figure 3- Air flows through an axial flow fan whose speed can be adjust- 
ed by means ol a U.S. varidrive motor. The range of speed of the fan is 
then 1100 to 3 ^ 00 rpm. The air then enters a flow straightener and tur- 
bulence management section, consisting of mesh screens and packed straws. 

Next is the contoured nozzle leading to the test section. This 
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axisymmetric nozzle was designed by the method of Morel to produce a 
minimum adverse pressure gradient in the boundary layer to avoid the flow 
unsteadiness phenomenon associated wi th local separation regions. The 
throat of this contoured nozzle is 15 cm onto which an annular vane 
swirler may be fitted. The swlrler consists of ten blades which may be 
individually set to any angle of interest. A solid hub of <4 cm diameter 
is located at the center of the swirler with a streamlined nose facing 
upstream. The downstream end is simply a flat face, simulating the geo- 
metric shape of a typical fuel spray nozzle. More detailed information 
concerning the swirler can be found elsewhere [13]. 

The test section itself consists of a 90 degree expansion block of 
diameter 30 cm and width 3-2 cm, giving an expansion ratio of 2. A 125 
cm plexiglass tube is then placed over the expansion block to provide 
the main body of the test section. In addition, another plexiglass tube 
of similar legnth but 15 cm diameter is used to provide an expansion ra- 
tio of unity. The substantial size of the test model provides excellent 
probe resolution for the hot-wire measurements. 

The test section has the option of being equipped with a strong con- 
traction nozzle of area ratio 4. This nozzle has a 45 degree slope fac- 
ing upstream and can be located at any axial location. Alignment of the 
contraction nozzle and test section with the wind tunnel is obtained with 
the use of a laser beam. 

A manual traversing mechanism can be attached to the test section 
as shown in Figure !i. This mechanism is equipped with a linear vernier 
scale which is accurately readable to within ±0.25 mm. This allows a 
probe to be accurately traversed across the chamber radius with the capa- 
bility of manually rotating a probe about its axis to an accuracy of 
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±0.2 degree. A more detailed description of the confined swirling jet 
facility at Oklahoma State University can be found in Reference C 1 5] - 

5-3 Hot-Wire Instrumentation 

As mentioned previously, the six-orientation single-wire technique 
employs standard electronic equipment for the determination of the time- 
mean and turbulent stresses in the flowfield. The sensor used in this 
study is a normal hot-wire probe, DISA type 55POI. This probe has two 
prongs set 3 mm apart which support a 5 um diameter tungsten wire. This 
wire is gold plated near the prongs to reduce end effects and strengthen 
the wire. A OISA 5 m cable connects the probe support to the main bridge 
unit. The anemometer used is a DISA type 55M01, constant temperature 
standard bridge. The mean voltage across the hot-wire is measured using 
a Hickok Digital Systems, model DP100, integrating voltmeter with the 
root-mean-squa re component of the voltage being measured using a Hewlett 
Packard, model A00 HR, AC voltmeter. 

To measure the eddy dissipation rate and hence the dissipation 
length scale, a special electronic analog device was constructed to time- 
different iate the fluctuating electronic signal. This device was assem- 
bled using a 7^1 operational amplifier and various capacitors and resis- 
tors- A schematic of the differentiator is shown in Figurel2a. For the 
device to effectively differentiate at the required frequencies, the 
break frequency of the instrument must be approximately ten times that 
of the differentiating frequency. The plot of gain versus frequency in 
Figurel2b shows that the cutoff frequency i s not reached untit after 20 
kHz. The calculated break point of the system is at a gain of 100; in 
addition, a gain of unity was obtained from the differentiator at a 
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frequency of 1,5 kHz. This is well within the range of frequencies to 
be differentiated. The output voltage from the anemometer was first 
passed through a DISA 55D25 filter before analog differentiation. The 
cutoff frequency setting for the wire is critical, since too high set- 
tings would include undesirable amounts of electronic noise and too low 
settings would suppress high frequency components of the signal. The 
settings were determined at typical locations in the flowfield by using 
a Hewlett Packard 3580A spectrum analyzer to identify regions of high 
noise. A unity gain follower was placed before the differentiator to 
prevent any high impedance loading of the differentiator. 

The mean square of the differentiated signal was measured using 
Thermal Systems Inc. true rms voltmeter, model 1076. This instrument 
has a IMfi input resistance to attenuate the low frequency signals. 


CHAPTER VI 


RESULTS AND DISCUSSION 

Nonswirling and swirling nonreacting flows are investigated in an 
ax.i symmetric test section with expansion ratio D/d - which may be 
equipped with a strong contraction nozzle of area ratio 4 located at L/D 
= 2. This nozzle has a 45 degree upstream facing slope as described ear- 
lier [10]. Measurements have also been made in a nonexpanding test sec- 
tion with expansion ratio D/d = I. The nozzle inlet velocities and Rey- 
nolds numbers employed in this investigation are high enough to ensure 
that the flowfields studied are under conditions independent of Reynolds 
number variation. They correspond approximately to conditions reported 
in associated studies [ 1 0] . In each case, radial profiles of interest- 
ing properties at six axial locations are presented. AM properties 
shown are normalized with respect to the swi Her inlet uniform axial 
velocity U q deduced independently from a measurement upstream of the 
swirler. Finally, an accuracy and directional sensitivity analysis in- 
volving the time-mean and turbulent properties of a round free jet has 
been performed. Percentage changes of some of the important input para- 
meters and their effect on the calculated output quantities have been 
studied. Five probe-jet axis orientations have also been used to simu- 
late the shifting of^he dominant flow direction and consequently the 
effectiveness of the hot-wire technique to measure the properties of a 
strongly swirling flow whose local flow direction is unknown. 
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Deduced results from the measured hot-wire signals are presented in 
Tables I through XXIII. 

6. 1 Effect of Sw i rl 
6.1.1 Nonswirling Flows (D/d = 2) 

Figure 13 shows time-mean and turbulence data for the nonswirling 
flow. A nearly-flat axial velocity profile is seen in the entrance re- 
gion of the test section. As expected, there is no measurable swirl 
velocity and only a small radial velocity. Although a corner recircula- 
tion zone (CRZ) is present, there is no evidence of a central toroidal 
recirculation zone (CTRZ) . Indeed, there is no swirl-induced centrifu- 
gal force to encourage its formation. The hot-wire anemometer cannot 
sense flow direction; however, directional properties of the flow can be 
inferred from earlier five-hole pitot probe data [10] , and flow visual- 
ization photography [ 9 ]. Nevertheless, the author has retained the posi- 
tive values on al 1 figures wherever possible uncertainties might exist. 
The time-mean data show good agreement, within 5 percent, wi th that found 
by Chaturvedi [17] in a similar test section, and wi th that found by Yoon 
and Lilley [10] using a five-hole pitot probe in the same test facility 
with identical flow conditions. 

The maximum values of normal stresses appear on the shear layer with 
the axial fluctuation component dominating. Earlier results [ 1 7l indi- 
cated that the axial turbulence intensity was larger than the other two 
components and that the radial turbulence intensity was approximately 
equal to the tangential turbulence intensity. This is confirmed in the 
present study. 
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The six-orientation technique produces positive values of shear 
stress. However, in certain locations in the vicinity of recirculation 
zones, the radial gradient of the axial velocity is predominately posi- 
tive which is associated with negative values of u'v’. Nevertheless, 
all shear stress values in this document are plotted as positive. Be- 
cause of the absence of velocity gradients in the x- and B-di rect ions , 
only one shear stress (the xr-component) is significant in the nonswirl- 
ing case. This shear stress, plotted in Figure 13, tends to be lower 
than the earlier study [17]. It should be noted, however, that the shear 
stresses are the most difficult turbulent quantities to measure accurate 
ly in a complex flowfield. Uncertainties in measurements of time-mean 
velocities and turbulence intensities are increased in the determination 
of shear stresses. 

6.1.2 Moderate Swirl (D/d = 2) 

in confined swirling jet flows the axial and tangential time-mean 
velocities dominate, as can be seen from Figures 14 and 15* The corner 
recirculation zone can be seen clearly at the expansion plane of the 
test section. The CRZ is not seen at any other axial location, indicat- 
ing that the swirling flow greatly reduces the length of the CRZ. Swirl 
i ng flow produces a central toroidal recirculation zone which can also 
be seen in the figures. The CTRZ appears to have a length of approxi- 
mately 1.5 0- Downstream of this point on the jet axis, indication is 
given of a precessing vortex core ( P VC) extending to the exit plane of 
the test section. The PVC is defined as a region of low axial velocity 
and high almost solid body rotation swirl along the axis, and is found 
to be present in the swirling flows considered. Yoon and Li 1 ley [10] 
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also measured the time-mean swi rl ing flowfield using a five-hole pitot 
probe, and the present data are found to be in good agreement. The PVC 
is most clearly observed in flow visualization studies using still photo- 
graphy [8, JO] and videotape recordings [9] . 

The three normal turbulent stresses appear to be fairly isotropic 
at all locations in the test chamber with maximum values occurring in 
regions of recirculation and regions of high shear. It was found that 
large-scale turbulence with big eddies occurs in recirculation regions 
and that small-scale turbulence with small eddies occurs in regions of 
peak velocities. In the downstream regions of the test chamber, the tur- 
bulence levels are low, with a more uniform radial profile, indicating a 
more developed nature of the flow. 

It can also be seen that all three shear stress components are sig- 
nificant in swirling flows, as expected. The maximum values of shear 
stress occur i n the thin shear layer regions but quickly dissipate in the 

. . 2 

downstream direction as the shear layer broadens. The stresses u'v'/u 

o 

2 

and u'w'/u are found to have large values close to the wail, because of 
o 

the steep axial and swirl velocity gradients. 

6.1.3 Strong Swirl (D/d = 2) 

Time-mean velocity profiles for the strongest swirl cases consider- 
ed are shown in Figures 16 and 1 7- Almost all of the flow leaves the 
swirler near the outer edge, producing steep velocity gradients in this 
vicinity. High velocity gradients can also be seen near the wall, es- 
pecially at x/D = 0.5* The strong centrifugal forces present in the in- 
coming flow produce rapid outflow to the confining boundary. Both cen- 
tral and corner recirculation zones can be seen clearly from the time- 
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mean plots. However, it appears that the CTRZ is shorter for th i s strong 
degree of swirl as compared to the moderate swirl case. In contrast, 
the PVC gets wider as the swirl strength increases, as also found infive- 
hole pitot probe data [10]. 

The normal turbulent stresses have increased in magnitude, consis- 
tent with the increase in swirl strength, and still a good degree of iso- 
tropy is observed throughout the entire flowfield. The highest turbu- 
lence levels again occur In regions of recirculation and on the shear 
layers. It can also be seen that high turbulence levels are found in 
the PVC. 

The most dramatic effect of the increase of swirl is the large in- 
crease crease in all three shear stress values. It can be seen that 
very high values of shear stress occur in the shear layers and near the 
walls. The PVC also contains high values of shear stresses and turbu- 
lence levels. Overall, the values of shear stresses are higher than in 
other swirl strengths considered. 

6.2 Effect of Strong Contraction Nozzle 
6.2.1 Nonswirling Flows 

Time-mean and turbulence characteristics for the nonswirling flow 
with a strong contraction nozzle at L/D = 2 are presented in Figure 18. 
The plots show that results vary only slightly from that of the corre- 
sponding flowfield without a contraction nozzle (see Figure 13). The 
major difference appears to be a slight reduction in the length of the 
CRZ. The measured time-mean flowfield compares favorably with previous 
data [ 10] . 
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6. 2. 2 Moderate Swi r 1 


The effects of the contraction nozzle on the moderately swirling 
flow with ij> = 45 degrees are shown in Figure 19. The presence of the 
contraction nozzle accelerates the flow and produces a strong favorable 
pressure gradient over the entire flowfield of interest. This pressure 
gradient conflicts with the adverse pressure gradient inherent in swirl- 
ing flows (associated with recirculation zones). Regions of positive 
axial velocity occur near the centerline at all axial locations. The 
central recirculation zone is now located in annular region around the 
jet axis and is much smaller than in the corresponding open flowfield of 
Figure 15- A narrow central core region is observed to extend through- 
out the length of the test section with strong solid body rotation. Posi- 
tive axial velocities now occur in this region, as opposed to negative 
ones in the corresponding open-efded flow case; this and other time-mean 
data are in excellent agreement with earlier experimental [10] and pre- 
diction [4, 16] studies. 

The directional turbulence intensities do not show any significant 
increase in magnitude compared to the open-ended swirling flow. However, 
the turbulent shear stresses are found to increase near the jet axis as 
the contraction nozzle exit is approached. This is because of the fairly 
high turbulence levels in this region and the effect of strong velocity 
gradients with which these stresses are associated. 

6.2.3 Strong Swi rl 

For swirl vane angle 4> = 70 degrees, measurements are given in Fig- 
ure 20 with the strong contraction blockage located at L/D = 2. The 
axial velocity near the axis is positive, though less so than in the 45 
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degree case of Figure 1 9 » since the favorable pressure gradient has nts. 
to overcome an even stronger unfavorable pressure gradient. The central 
recirculation region is now very small, extending in an annular region 
to less than x/D = 1.0, considerably less than the no-blockage case of 
Figure 1 7 . At the axial station x/D = 1.0, forward flow occurs across 
the whole test section. Very strong swirl velocity magnitudes and gradi- 
ents are seen, which contrast sharply with the corresponding open-ended 
flow situation. A wide core region is again noticed along the jet axis 
containing strong solid body rotation; movie photography reveals the pre- 
cessing nature of this phenomenon. Again, time-mean data compare very 
well with previous work [10], 

The normal components of the Reynolds stress tensor show an increase 
in turbulence alcng the jet axis as the contraction is approached--more 
so than in the AS degree swirl case of Figure 19, but similar to those 
found in the open-ended case of Figure 1 ?, N ow, larger values are found 
near the axis, associated with the strong vortex core region. Shear 
stress levels at the entrance to the test section tend to be slightly 
lower for the blockage case. But the levels are an order of magnitude 
higher in the core region near the contraction nozzle. Again, this is 
because of the strong time-mean velocity gradients in this area. 

C.3 Effect of Expansion Ratio 
6.3.1 Nonswirling Flow (D/d = 1) 

Figure 2! shows the very flat axial time-mean velocity profiles 
associated with a one-dimensional flow. The flowfieid is seen to be 
homogeneous throughout the whole test section. This is because of pres- 
sure "losses due to abrupt expansions being absent in the nonexpanding 
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confined jet. Radial velocities are found to be present across the 
whole test section. This was not to be expected as there is no expan- 
sion of the flow in this case. These velocities are assumed to be an 
error of the six-orientation technique which is discussed later in this 
chapter. No recirculation zones are evident at any section of the flow- 
field, as expected. 

The successful design of the wind tunnel and contraction nozzle is 
evident by the low levels of turbulence intensity measured. The absence 
of time-mean velocity gradients and flow discontinuities also explains 
the fact that all components of the Reynolds stress tensor are negligi- 
ble. 

6.3.2 Moderate Swirl (D/d = 1) 

Strong axial and tangential time-mean velocities dominate in the 
nonexpanding confined jet. These velocities were found to be approxi- 
mately twice as large as those in the corresponding expanding flowfield 
(D/d = 2). A momentum balance across the test section indicates that 
these increases in magnitude are to be expected [70] . These results can 
be seen in Figure 22. This is because the wall of the test section pro- 
duces a favorable pressure gradient across the entire flowfield which 
conflicts with the adverse pressure gradient within the swirling flow. 
However, a core region around the jet centerline is found to contain 
strong solid body rotation and positive axial velocities. Evidence can 
be seen from the time-mean velocity plots of the central hub of the 
swirler. Negative radial velocities were expected and measured in this 
particular flowfield. 
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Directional turbulence intensities are found to be similar in mag- 
nitude for this case as in the expanding flowfieid (D/d = 2) except in 
the entrance region. This is because of the lack of regions of recircu- 
lation normally associated with high turbulence levels. The turbulence 
in this flowfieid is primarily caused by passage of the air through the 
swlrler and because of the velocity gradients. Because of these velo- 
city gradients, shear stresses are significant throughout the entire 
flowfieid. The magnitude of these stresses were found to be comparable 
to those in the expanding case. 

6 . 3.3 Strong Swirl (D/d = 1) 

It can be seen from Figure 23 that the tangential time-mean velo- 
city completely dominates the flowfieid. This strong swirl velocity pro- 
duces a large core region with low axial velocity on the jet centerline. 
This is because of the centrifugal forces within the flew. Again, no re- 
circulation zones are present within the flowfieid. Radial velocities 
have increased due to the increase in swirl velocity. 

The large gradients associated with the core region produce high 
levels of turbulence intensity along the jet centerline and correspond- 
ingly a large increase in all of the shear stresses in this region. 

Shear stresses in the entrance region of this flowfieid are much lower 
than in the expanding case. Again, this is because of the lack of recir- 
culation zones. 

S.k Turbulence Parameters 

The eduy dissipation rate measurements have been performed in flows 
with two different expansion ratios with three degrees of swirl. The 
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levels tsf confinement considered are with expansion ratios of D/d = 2 
and I with swirl strengths of nonswirling, moderate, and strong swirl. 

The measurements are obtained using the time differentiating technique 
described in Chapter III. Estimates of the different length scales of 
the flows can be obtained from the measurements of the dissipation rates. 
The eddy dissipation length scale (macroscale) is estimated via: 

1,3/2 

l = £ 

e 

the Taylor microscale is deduced from the isotropic relation: 


and finally the Kolmogorov length scale is giv.en by: 



2 3 

The eddy dissipation, e, has dimensions of L /T , the same as velo- 
city cubed divided by a length scale. At any axial station, the station 

2 

mean axial velocity (* U (d/D ) ) and chamber diameter D are used for 

v 

this purpose, Accordingly, plots of 


u 3 dV 

o 


are given. The dissipation iength scale is nond i mens i ona 1 i zed with re- 
spect to the test chamber diameter D. 


6.4.1 Expanding Fiowfield (D/d = 2) 

Measu; -.ments of the turbulent dissipation rate in a confined, ex- 
panding, nonswirling jet are shown in Table XII and Figure 24. The 
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plots show that most of the dissipation occurs on the thin shear layer 
in the entrance region. As the shear layer spreads across the entire 
test section, the corresponding dissipation rate follows. The dissipa- 
tion rate is also seen tc be high in the CRZ. As expected, the dissipa- 
tion rates in the essentially laminar portion of the potential core are 
very low. These low dissipation rates in the laminar core region tend 
to cause high values of length scale, as seen in Table XIII. The length 
scales are largest on the jet centerline, outside the potential core, 
and then decrease as the test section is radially traversed. This is be- 
cause the dissipation rates are low on the centerline, then increase 
radially; but the kinetic energy of the turbulence is somewhat constant 
across the test chamber. 

'(hen a moderate degree of swirl is introduced to the flow, the dis- 
sipation rate is found to increase significantly in the entrance region 
of the test section (Figure 24 and Table XIV). However, these values 
are rapidly reduced in the downstream region of the chamber. This im- 
plies that most of the production and dissipation of the turbul ent energy 
is being performed in the entrance region of the test secti on where areas 
of recirculation and high shear dominate the flow. On the shear layer 
eddies are formed which contribute to the large rates of dissipation. 
Within the recirculation zones large scale eddies cascade energy to its 
smaller scales which eventually leads to its dissipation.. 

The measurements of dissipation rate in the entrance region of the 
test section must only be considered as qualitative as several sources 
of error may be present at these axial locations. In the recirculation 
zones the time-mean velocity direction is not well defined; therefore, 
perfect alignment of the sensor in the mean flow direction is not assured. 


54 


Two other distinct types of error were certainly present. The first was 
discussed by Kovasznay [ 71 ] and is the one arising of loss of response 
of the anemometer equipment to high frequency turbulence components. The 
anemometer system has a cutoff frequency of approximately 30 kHz which 
is below the frequencies at which turbulent dissipation occurs. The mea- 
sured dissipation rate must therefore be extrapolated to its full rate 
by assuming a theoretical turbulence spectrum. Pao [72] suggested such 
a spectrum but it has yet to be conclusive in free or confined shear 
flows. Estimates of the Kolmogorov length scale indicate that it is much 
smaller than the hot-wire length. This leads to the second source of 
error of having to average over the hot-wire length. Wyngaard [73] has 
shown that the underestimation of the dissipation rate is in the order 
of 6 percent for the present flow conditions. As this is considered 
within the scatter of the data, no hot-wire length corrections to the 
measured data have been made. 

The distribution of length scale for <j> = 45 degrees can be seen in 
Table XV. The table shows that low values of length scale can be expect- 
ed for moderately swirling flows, especially in regions of high accelera- 
tion. further downstream the length scales are close to unity, indicat- 
ing that the flow is probably of a more uniform nature. The low dissipa- 
tion rates at the exit of the test section also lend to this conclusion. 

The effect of increasing the swirler angle to 70 degrees can be 
seen in Figure 2^ and Table XVI. A note of interjection: the wind tun- 
nel test facility cannot provide equal Reynolds number flows for all of 
the swirl strengths considered because of the blockage caused by the 
swirler vanes. Therefore, the nond imens ional i zi ng inlet velocity, U o> 
varies for each cf the swirl strengths. As this value is raised to the 
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third power, this number becomes very significant and hence the change 
in scales of Figures 24 and 25- Nevertheless, the dissipation rates in 
the entrance region of the test section with 70 degrees of swirl are 
still higher than the corresponding cases with 45 degrees and zero swirl. 
Peak regions of dissipation still occur in the same regions as in the 45 
degrees of swirl flowfield, but the reduction of dissipation rates across 
the test section is more rapid. 

The dissipation length scales, shown in Table XVII, are smaller 
than the other flow cases considered, mainly because of the existence of 
larger recirculation zones and greater shear. 

6.4.2 Nonexpanding Flowfield (S/d - 1) 

Eddy dissipation measurements in nonexpanding confined jets are pre- 
sented in Tables XVII! through XX! 11 and Figure 25- The same scales are 
kept on the axis of the plots to correspond to the expanding flowfield. 

As can be seen from Figure 25, the eddy dissipation rate is greatly de- 
creased for the nonexpanding flow. This can be explained by the fact 
that there are no areas of recirculation present for any of the swirl 
strengths considered. However, on the centerline of the test section 
the dissipation rates are independent of expansion ratio for the swirl- 
ing flows. This is because velocity gradients are present due to the 
tangential velocity. 

The dissipation length scales in the nonexpanding flowfield tend to 
be generally lower than In the expandi ng case but more uniform across 
the entire test section. Table XXI II shows that these length scales are 
smaller than the test chamber diameter, indicating that turbulent dissi- 
pation i s occurring. 
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6.4.3 Important Turbulence Properties 


Other properties in turbulence modeling are the kinetic energy of 
turbulence, k, and the turbulent dissipation length model, £. These are 
defined by 


k = { u * 2 + v’ 2 + w ' 2 )/2 


% = k 3/2 /e 


and can be obtained from data presented earlier, but are conveniently 
deduced and plotted now. 

Plots of the kinetic energy of turbulence for two varieties of ex- 
panding flowfields for a range of swirl strengths are shown in Figures 
26 and 27- These figures show that most of the kinetic energy of the 
turbulence in the nonswirling expanding flowfield is occurring in the 
downstream regions of the test section away from the laminar-type flows 
encountered in the entrance region. In the swirling flowfield high 
levels of kinetic energy are observed ?n two distinct regions. One is 
the entrance region vfh-are there are two thin shear layers. The other is 
in areas of recirculation including the PVC. 

In the nonexpanding flowfield the levels of kinetic energy are very 
low in the nonswirling flow. This is because of the absence of any time- 
mean velocity gradients except very close to the confining wall boundary. 
In the swirling flowfields energy levels are again fairly large in the 
entrance regions and in the core region around the jet centerline. It 
can be seen from both Figures 26 and 27 that outside areas of recircula- 
tion and shear the radial profiles of the kinetic energy are fairly uni- 
form, indicating a developed nature of the flow. 
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Figures 28 and 23 , and Tables Xil, XV, XVII, XIX, XXI, and XXII 
show plots of the calculated dissipation length scale, in the expanded 
and nonexpanded geometries, respectively. The accuracy of these plots 
is dependent upon the accuracy of the measuring techniques used for de- 
termining the dissipation rate and the directional turbulence intensi- 
ties. In deducing the kinetic energy of turbulence the root-mean-square 
of the three turbulence intensities is measured. Therefore, any error 
in these measured quantities is compounded in the determination of the 
kinetic energy. In calculating the dissipation length scale the kinetic 
energy is raised to the power 3/2, this again magnifying any error in the 
measurements. Any errors likely to be found in the eddy dissipation 
rate measurements have been discussed previously in section 6.4.1. It 
is expected that the entrance region of the test section would be an area 
of the flow to give the least accurate measurements of the dissipation 
rate. It can be concluded that areas of high shear and recirculation 
(and hence high turbulence levels) will produce the worst estimates of 
dissipation length scales. 

The plots presented here show that the largest dissipation length 
scales can be expected in the nonswirling expanding flowfield, especial- 
ly in the downstream regions where the kinetic energy is high and the 
dissipation rate is fairly low. In the entrance region of the swirling 
flowfield the dissipation rates are very high with most of the turbulent 
dissipation quickly spreading over the entire test section. This leads 
to small dissipation length scales in entrance regions of the flowfield. 

In the non5wirling, nonexpanding flowfield the dissipation length 
scales are very low across the entire test section. This is because of 
the very low turbulence levels throughout the flowfield. For the 
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swirling flows the length scales are small in the entrance region but 
Increase in size inside the core region as the exit to the test section 
is approached. This is because of the large values of kinetic energy 
found In the core regions. 

Overall, it is felt that these plots of the dissipation length 
scale are indicative of the levels of turbulence transport found i n swirl- 
ing confined flows with the downstream regions of the test section giv- 
ing the more accurate quantitative data due to the inherent problems of 
hot-wire anemometry in the entrance region of the test chamber. 

6.5 Uncertai nty Analysi s 

The uncertainty analysis includes a determination of the sensitiv- 
ity of the six-orientation hot-wire data reduction to various input para- 
meters which have major contributions in the response equations. Pitch 
and yaw factors (G and K) are used in the response equations described 
in Chapter lit in order to account for the directional sensitivity of 
the single hot-wire probe. Figure 6 shows the pitch and yaw factors 
plotted against the hot-wire mean effective voltage. Both the pitch and 
yaw factors are functions of the hot-wire mean effective vol tage, but the 
yaw factor is far more sensitive. A one percent i ncrease In the hot-wire 
voltage reduces the pitch factor by 1 . 3 percent and the yaw factor by 56 
percent. For the present study the values of these factors are chosen 
at an average hot-wire voltage experienced in the flowfield. This was 
appropriate since the output quantities (u, u' , u'v 1 , etc.) are only 
weakly dependent on the value of K. This can be seen in the data of 
Table XXIV which summarizes an analysis performed on the data reduction 
program at a representative position in the flowfield. 
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Table XXIV demonstrates the percent change in the output quantities 
for a 1 percent change in most of the important input quantities. For 
the data presented in this table, only quantities calculated from the 
probe orientation combination 5, 6, and 1 are used, for simplicity. The 
situation is that of a moderately swirling confined flowfield from a 
swirl generator with vane angle of 38 degrees. In this swirling flow 
orientation, 6 was the minimum of the six mean effective cooling veloci- 
ties. King [^8] has argued that the probe orientation combination ap- 
proximately centered around the minimum effective cooling velocity pro- 
duces more accurate estimates of calculated turbulence quantities than 

do the other orientation combinations. However, all previously reported 
data have been obtained by averaging all the six possible combinations. 

It is not unusual in hot-wire anemometry to have the mean velocity 
components and turbulence quantities that are measured be quite sensi- 
tive to Jiges in mean hot-wire voltage. For interpretive purposes, 
the mean hot-wire voltage variations can be thought of as being either 
errors in measuring the mean voltage, or shifts in the individual wire 
calibrations due to contamination or strain "aging" of the wire. The 
data of Table XXIV demonstrate that the most serious inaccuracies in the 
measurement and data reduction technique will be in the estimates of tur- 
bulent shear stresses , the most inaccurate output term being u 'w' . 

As already discussed in Chapter III, an ad hoc assumption is made 
regarding the numerical values of the correlation coefficients used in 
the deduction of time-mean and turbulence quantities. The results of 
the uncertainty analysis (Table XXIV) show the time-mean and turbulence 
quantities to be relatively insensitive to variations in the correlation 
coefficients. Therefore, the major ad hoc assumption made in the 
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technique does not seem to have a great effect on the output quantities 
compared to the effect of other input quantities. 

The accuracy of the hot-wire technique is dependent upon two types 
of input parameters. One set is the measured input data and the other 
is the assumed data of the reduction procedure. The accuracy of the re- 
sulting deduced values are from the d.c. voltages (Ep E^» E^) , the rms 

of the fluctuating voltages (Ej , E‘ , E£ ), the pitch and yaw 

rms 3 rms rms 

factors (G and K) , and the correlation coefficients (y ). Each of 

2 i Z j 

these quantities is not equally likely to have a I percent error. For 
example, an error of 1 percent in the mean voltage is very large and is 
more realistically in the order of 0.25 percent. However, it is useful 
to appreciate the likely error in each of the input parameters , an assess- 
ment of which is now given. 

The quality of the measured input data of mean voltages and rms 
fluctuating voltages depends upon the instrumentation used and the com- 
plexity of the flowfield at the measuring location. In laminar stable 
flows, these values are likely to be very accurate (less than 0 . 25 %), 
whereas In turbulent shear regions the mean voltages are measurable to 
0.5 percent and the rms voltages to 1.0 percent. Pitch and yaw factors 
are obtained from calibration curves through a sequence of measurements 
made in a very stable flowfield in the potential core of the calibration 
jet. They are not likely to be in error by more than j. 2 percent. 

The assumed input parameters of the data reduction, the correlation 
coefficients, could be in error by as much as 5 percent. This is be- 
cause in most turbulent flows the correlation coefficients are usually 
not less than 0.8 and no greater than 0.9. As mentioned in Chapter III, 
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the correlation coefficient vlaues of 0.867 are obtained from measure- 
ments in a round free Jet, a flow similar to the type being studied. 

Overall, the worst possible scenario would be the cumulative effect 
of reinforcement of errors in these parameters as follows. In laminar 
regions the maximum errors are: 9 percent for time-mean velocities, 15 

percent for normal stress, 22 percent for shear stresses, and 55 percent 
for the x6-component of the shear-stress tensor. In turbulent shear re- 
gions the maximum errors are: 18 percent for time-mean values t 2** per- 

cent for normal stresses, and 29 percent for shear stresses, except 98 
percent for the x9-component. 

It may also be noted that absolute accuracy of the technique depends 
also on the physical size of the probe interfering with the flow, the 
position of the wire at each of its rise orientations, and the effect of 
the dominant flow direction on the measured cooling velocities. These 
are assessed in the following section with section 6.9 concluding the 
discussion. On this basis, it is felt that maximum errors of 5, 10, and 
15 percent are realistic deviations from the true time-mean, directional 
turbulence intensity, and turbulent shear-stress values over most flow- 
field locations. 

6.6 Directional Sensitivity Analysis 
in Laminar Jet 

The directional sensi tivi ty of the technique is assessed at the five 
locations A through E (see Figures G and 9) corresponding to five differ 
ent flow situations. The first of these is in the laminar potentional 
core region, at x/d = 0 and r/d = 0. Table XX, V give- the results with 
the probe coordinates aligned with facility coordinates, as Case 1 of 


Figure 7 illustrates. This is referred to as Situation A Case l with 
analogous statements later. The time-mean velocities, nondimens ional i z- 
ed with the jet exit velocity deduced from an independent measurement, 
are shown. In this one-dimensional flowfield the axial velocity is ex- 
pected to be unity with the other two components of the velocity vector 
to be equal to zero. Results using each of the six possible combina- 
tions of three adjacent wire-orientations are presented, together with 
the mean of the values. The standard deviation and its ratio with the 
mean are also presented to show the amount of scatter in the readings. 

As can be seen, the error for the axial and swirl velocities are very 
low for each combination. The radial velocity error tends to be larger, 
possibly because of slight probe misalignment with the normal to the jet 
axis. The mean of these quantities brings the data to well within accep- 
table limits. 

Results of the probe being rotated by 45 degrees about the z-axis 
(as in Case 2 with 0 = -45°) are shown in Table XXVI. This is Situation 
A Case 2. The probe coordinate system is now different from the jet co- 
ordinate system but the measured velocities can be related to the facil- 
ity coordinates by use of the rotational matrices given in Chapter IV. 

The values given in Table XXVI and indeed all the tables are presented 
in terms of the facility coordinate system. The results show that this 
misalignment of the probe with the dominant flow direction still gives 
excellent values of velocities in the laboratory coordinate system with 
the use of any of the six possible wire combinations. Consequently, 
averaging of the data also gives good results. 

Rotation about the z-axis by -90 degrees to obtain Case 4 results 
in no data being generated by the technique. This is because in a 
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steady one-dimensional flowfield, all the instantaneous cooling veloci- 
ties acting on the hot-wire ai"e equal for all six rotations. As the 
data analysis requires that their cooling velocities be subtracted from 
pi?ch other (see Reference [41]), all the output terms are deduced as 
zero. 

To further investigate the shifting of the dominant flow direction, 
the probe was rotated twice (-45° about its z-axis, followed by -45° 
about its new x-axis) so as to conform to Case 3- The results of these 
axis rotations can be seen in Table XXVII. Again the laboratory coordi- 
nate deduced values and deviations from expected values are relatively 
low, although not quite as good as in the previous cas^. The advantages 
of averaging can be seen in Table XXVI i, where the under- and/or over- 
estimation of the velocities for the individual positions are smoothed 
out after averaging. 

Rotations of 0 = -90 degrees and <£ => -90 degrees were also carried 
out at the same flow location, thus obtaining Case 5. Table XXVII l 
shows the results of these rotations. Good axial and tangential velo- 
city values can be seen but with a decrease in the accuracy of the radi- 
al velocity. 


6.7 Directional Sensitivity Analysis in 
Turbulent Nonswirling Jet 

Similar rotations of the probe about its axes have been performed 
on the shear layer of a free nonswirling axi symmetric jet at two axial 
locations, x/d » 3 and 10, referred to as Situations B and C. At these 
points in the flow the axial velocity dominates with a small contribu- 
tion from the radial velocity. In axi symmetric jets it is well known 
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that the axial directional turbulence intensity is larger than its other 
two components [74] . The only significant shear stress inthis fiowfield 
is the rx-shear. Tables XXIX and XXX, obtained with Case I probe con- 
figuration in Situations B and C, confirm this. Incidentally, all the 
data presented in this study are a consequence of a typical set of read- 
ings obtained from the hot-wire technique. If a large amount of scatter 
is found in deduced results at a particular location, the problem is fur- 
ther investigated and/or remeasured before being accepted as valid. Data 
reported elsewhere [4l] have been thoroughly ana I y; ad and checked for re- 
peatabi 1 i ty. 

Tables XXIX and XXX are used as a standard for nonswirling flow- 
field values at locations B and C, so as to be able to compare the re- 
sults obtained from other probe configuration cases. However , the radial 
time-mean velocity appears to be very large for this particular flow- 
field and could possibly be in error. The coefficient of variation 
(ct/x) is seen to be acceptable for most of the flow properties except 
for the shear stresses. These large variations are caused by the shear 
stresses being two orders of magnitudes lower than the time-mean veloci- 
ties. Sometimes the data reduction will not resolve a particular para- 
meter. This is usually the consequence of subtracting two almost equal 
effective cooling velocities as described earlier. If a large propor- 
tion of the data is not resolved from the different combinations of cool- 
ing velocities used, the parameter deduced is taken to be zero. 

Results for probe configurations of Cases 2 and k for location B 
are given in Tables XXXI and XXXII; corresponding data for ’ocation C 
are given in Tables XXXIII and XXXIV. It can be seen, from Tables XXXI I 
and XXXIV, that the Case 4 results give poor estimates of the time-mean 
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radial and swirl velocities and underestimates of the axial components 
of the velocity vector. This is expected because the time-mean velocity 
vector is almost parallel to the hot-wire support axis, and rotating the 
probe through its six orientations yields little effect on the sensed 
data. This also decreases values of the axial turbulence intensity but 
significantly increases the other two normal stress components. Accep- 
table values of the axial normal stress are obtained via 8 = -45 degrees 
(Case 2 probe configuration) but again the tangential and radial turbu- 
lence intensities are increased; however, not as significantly as in the 
previous Case 4. Shear stress values for both Cases 2 and 4 are found 
to be very poor with large coefficients of variation. Also, many oi the 
turbulence quantities are not resolved. 

It would appear from these results that the six-orientation hot- 
wire technique is a poor tool to use if the flow is dominately in the 
direction of the probe support, for. reasons just described. If this 
occurs, simply reconfiguring the probe holder versus flow direction can 
overcome the problem. 

Tables XXXV through XXXVIII correspond to the last four tables, but 
now probe configurations of Cases 3 and 5 are used at locations B (Tables 
XXXV and XXXVI) and C (Tables XXXVII and XXXVIII). Now the time-mean 
velocity components are seen to be in excellent agreement with the val- 
ues determined from the standard configuration of Case 1. The axial nor- 
mal stress tends to be underestimated at x/d = 3 and overestimated at 
x/d - 10, relative to the Case 1 calculations. The radial turbulence 
intensity is consistently overestimated for both compound configurations 
dhd both axial locations. This could infer a failing of the hot-wire 
technique. The tangential turbulence i ntens i ty measuremen ts are found 
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to provide acceptable results. The dominant shear stress (the rx-compo- 
nent) in this flow is found to measure very well in the configurations 
of Cases 3 and 5, relative to Case 1 values. The coefficient of varia- 
tion is pot too large considering the magnitude of the numbers involved. 
The f v ria 1 components of the Reynolds stress tensor, although appearing 
co be measurable, exhibit a great dec’ f.-' scatter, perhaps indicating 
that these values are close to zero. 

6.8 Directional Sensitivity Analysis in 
Turbulent Swirling Jet 

Two free swirling jets were considered for further assessment of 
the hot-wire technique. The air flow exiting from an axisymmetric noz- 
zle of a wind tunnel passes through a vane swirler with 10 adjustable 
flat blades. The test facility and time-mean performance of the swirler 
are described at length elsewhere [i3]> For the present study the subse- 
quent large chamber confinement was removed and the free jet flow alone 
was studied, with swirl vene angles of 45 and 70 degrees being used. 
Figure 9 gives the specific measurement locations D and E used for mea- 
suring the 45 and 70 degree swirl situations, respectively. These loca- 
tions are in the high shear region of the flow close to the swirler exit. 
They were chosen since it was expected that all six components of the 
stress tensor would be significant, thereby providing a good test of the 
technique for their measurement. The hot-wire was also placed well away 
from the 'Jdge of the recirculation zone, so as to avoid any instantane- 
ous flow reversal on the wire. Rotations of the probe axes have been 
performed conforming to Cases 4 and 5 only, for both of the swirl strengths. 
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The standard six-orientation technique in the configuration of Case 
I produces the properties of both swirling jets as given in Tables XXXIX 
and XL. As can be seen, ali components of the time-mean flow and the 
Reynolds stress tensor are evaluated. The two sets of results are not 
quite at the same flowfield location because of the change in size and 
shape of the recirculation zone as swirl strength increases. However, 
the increase in the turbulent properties of t?vi flow are clearly evident. 

Rotation about the z-axis by -90 degrees to Case A probe configura- 
tion causes a deterioration in the accuracy of the results obtained from 
the technique, as inspection of Tables XL I and XL I I for **5 and 70 degrees, 
respectively, shows. The axial and swirl time-mean velocities are still 
fairly accurate for both flows but the radial velocity has suffered a 
large increase (or decrease) relative to its measurement with Case 1. 

.The normal stresses deduced after the rotation appear to be reasonable 
except for the radial and swirl components in the strong swiri flow of 
Situation E. These two components are greatly overestimated withaigreat 
deal of scatter in the results. It is again fel t that these poor results 
are because of the technique's inability to measure accurately flow 
properties when the dominant flow direction is in the direction of the 
probe holder. That is, when a large velocity is approximately normal to 
the wire in each of the six measuring orientations, insensitivity re- 
sults as discussed earlier. Correspondingly, the shear stresses also 
show a reduction in accuracy, with all three components either over- or 
'underpredicting the Case 1 values. 

The results of the compound angle of Case 5 are presented in Tables 
XL I ! I and XLIV, respectively, for the 1*5 and 70 degree swirl cases:. The 
axial time-mean velocity is seen to be good when compared to the standard 
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case, but the other two components show a reduction in accuracy in these 
highly turbulent flowfields. The inaccuracy of the radial velocity has 
been discussed earlier. The three components of the turbulence inten- 
sity appear to be fairly good with reasonable values deduced compared to 
the standard Case 1 values. Again, however, the radial and tangential 
components are less than desirable for the strongly swirling flowfield. 
The shear stresses for the 45 degree swirl situation D are considered 
good except for the u'w 1 term. This component is subject to great in- 
accuracies for only slight errors in these input data, as described ear- 
lier in section 6.5. for the strongly swirling flow (Situation E) the 
measured stresses are not so good when compared with the standard Case 1 
measurements. As these tables indicate, the time-mean velocities and 
normal stresses are subject to error which is magnified in the determina- 
tion of the shear stresses. 

6.9 Assessment of the Technique 

General results of the present and previous studies are now assess- 
ed in connection with the applicability, accuracy, and directional sensi- 
tivity of the six-orientation single normal hot-wire technique. 

Previously, In his measurements of strongly swirling vortex flows, 
King [48] compared his time-mean velocity and normal stress measurements 
with corresponding measurements obtained using a Laser Doppler Veloci- 
meter. He found excellent agreement indicating the validity of the 
method. He was not able to compare shear stress measurements in his 
swirl flow; however, he was unable to use his LDV for this purpose, in 
fact, despite the existence of advanced multicolor LDV systems and their 
use for shear stress measurement, no one has yet reported them in highly 
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swi r 1 i ng f low si tua t i ons , certainly not one over a range of swirl strengths 
as reported in this thesis. 

The present measurements have been compared with previously avail- 
able data whenever possible. In the nonswirling confined jet case, re- 
sults for time-n.-an velocities u and v, normal stresses u‘ and v' 

rms rms 

and shear stresses u'v 1 compare very favorably (see Reference [41], Fig- 
ures 7 and 8) with those of Chaturvedi [173. He used a cross-wire probe 
for the shear stress measurements. So also did McKi 1 lop [12] for non- 
swirling confined flow in the same facility as the present research. Re- 
sults, with and without exit nozzles, are in good agreement for the above 
quantities (see also Reference [12], Figures 21 through 28 ). The larg- 
est errors of approximately 15 percent occurred with the shear stresses 
in regions of high gradients. 

In the swirling confined jet case, comparison with Janjua and 
McLaughlin [75] for a moderate swirl strength in an identical facility 
has been possible. They made triple-wire hot-wire measurements in a 
flow with an inlet swirl vane angle = 38 degrees, using analog Lu- 
digital signal conversion and computer data reduction. For this purpose 
it was necessary to know in advance the local time-mean velocity vector 
direction; the data of Yoon and Li Hey [10] were used for this purpose. 
Their measurements [75] of time-mean velocity compare very well, within 
5 percent, with those of Reference [10] and hence of Reference [41]. Their 
measurements [75] of the three normal Reynolds stresses and the three 
shear Reynolds stresses are compared to x/D = 0.5, 1,0, and 1.5 with the 
six-orientation single-wire measurements of Reference [4l]. There is ex- 
cellent agreement, with a maximum discrepancy of 8 percent (see Refer- 
ence [75], Figures 10 through 18), indicating again the validity of the 
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present measurement technique. It appears to be an extremely viable, 
cost-effective technique for turbulent flows of unknown dominant direc- 
tion, requiring only accessibility to the measurement location. Probe 
interference appears not to be a major problem. 

For the present study. Figures 26 through 30 summarize measured 
values for the five situations A through E, respectively. Each figure 
presents facility coordinate time-mean velocity, normal and shear stress 
values obtained with each of the five probe holder versus facility con- 
figuration possibilities of Cases 1 through 5- A remarkable observation 
is that, in general, the configuration is of little importance,- results 
appear quite constant across the five cases. 

On the other hand, present production run results have used the 
Case 1 configuration exclusively. In those results, generally, ensemble 
average of deduced results from each of the six possible combinations of 
three adjacent wire orientations has been used. This was because of 
lack of local flow directional knowledge--! f this knowledge is available 
it is expected [48] that the combination with minimum cooling velocity 
in the central of the three wire orientations used will produce more 
accurate estimates of deduced flow quantities. From results of the pre- 
sent study for Situations A through E, Table XLV confirms this, especial 
ly for time-mean values. In any case, the appropriate choice of wire 
orientation for minimum cooling velocity is not known a priori. The 
values given in this table could only be determined after the measure- 
ment. However, for turbulence quantities and in the 45 and 70 degree 
situations, more confidence may be placed in the average of all possible 
wire combinations. This smoothing has been used exclusively in the pre- 
sent study. 


71 


CHAPTER VI I 


CLOSURE 

7*1 Conclusions 

The present research is concerned with the time-mean and turbulence 
properties of a low speed, nonreacting, swirling flowfield. The measure- 
ments have been performed in a model of an axisymmetric gas turbine com- 
bustor with varying degrees of swirl and geometry. The six-orientation, 
single hot-wire technique has been used to produce measurements of the 
time-mean velocities and complete Reynolds stress tensor in all of the 
flowfields of interest. The turbulent eddy dissipation rate and hence 
the dissipation length scale have been deduced for selected flowfields 
via time differentiation of the fluctuating hot-wire signal. 

The existence, size, and shape of both the corner recirculation 
zone and the central toroidal recirculation zone are affected by the de- 
gree of swirl. Measurements show that increasing the swirl strength 
from zero to medium swirl produces shortened corner regions and the gen- 
eration of a central bubble extending to approximately 1.5 chamber dia- 
meters. downstream of the jet exit. - In addition, a precessing vortex 
core is observed on the jet centerline which stretches to the exit of 
the test section. Increasing the swirl strength enlarges the central 
zone and vortex core with negligible effect on the corner region. 

The introduction of swirl to the confined, expanding jet flowfield 
greatly increases the turbulence properties of the flow. High levels of 
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normal and shear stresses are found on the two shear layers in the en- 
trance region of the test section. Dissipation rates are also found to 
be large in the high shear regions of the flow. These rates reduce and 
spread over most of the flowfield as the exit to the test section is ap- 
proached. 

Placement of a strong contraction nozzle at l/D = 2 with an area re 
duction of 4 causes a significant effect on the time-mean swirling flow- 
field. Central recirculation zones are shortened and axial velocities 
along the entire jet axis become positive. The core regions become nar- 
row with strong swirl velocities and gradients. Turbulence levels and 
shear stresses are found to increase along the jet centerline near the 
exit of the contraction nozzle but decrease slightly at the jet exit. 

No recirculation zones are found to be present in the nonexpanding 
flowfield with uniform velocity profiles across the entire test section. 
The swirling flowfield? exhibit strong tangential velocities compared 
with the corresponding expanding flow. A core region is still evident 
with almost solid body rotation along the axis. Turbulence levels are 
found to be low for this flow except near the jet centerline where 
strong velocity gradients are present. 

The ccuracy and directional sensitivity of the six-orientation hot 
wire technique has been performed in axisymmetr ic and swirling free jets 
Variation of input parameters and their effect on the output data has 
shown that the least accurate output quantities are the shear stresses, 
in particular the 0-component. 

The directional sensitivity analysis has shown that the technique 
adequately measures the properties of a flowfield independent of the 
dominant flow di rect ion except when the flow is predominately in the 
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direction of the probe holder, with the six-orientations of the probe 
creating insignificant changes in hot-wire response. The analysis shows 
that this component of time-mean velocity is inadequately deduced. Only 
reconfiguration of the probe can overcome this problem a posteriori. 

7.2 Recommendations for Further Work 

Fundamental research work should be continued in swirling flows in 
several areas. First, the introduction of a round nonswirling jet into 
the existing flowfield should be studied to enable a greater understand- 
ing of the performance of a gas turbine combustor. This would involve 
measurements of the time-mean and turbulence properties of the resulting 
flowfield. Second, the turbulent dissipation rate would have to be de- 
termined in such a flowfield before successful modeling of combustors 
could be achieved. This may involve improving the existing technique to 
account for the high frequency turbulence components and the averaging 
over the hot-wire length. 
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O . 02500 0.96335 0.94240 0.96926 0.89265 0.61281 0.66329 

O.OOCOO 0.97234 0.94179 0.95616 0.91932 0.61995 0.67986 



TABLE I (Continued) 



rms 


TABLE I (Continued) 



rms 



i ABLE I (Continued) 



. 00020 0 00332 0.01791 0.008'lfl 


TABLE II 

TIME-MEAN AND TURBULENCE DATA FOR SWIRLING 
FLOW 6 = 38°. D/d = 7 
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- 0.07267 0.11461 - 0.17983 - 0.15610 - 0.16447 - 0.16890 


TABLE it (Continued) 
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TIME-MEAN AND TURBULENCE DATA FOR SWIRLING 
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TABLE V (Continued) 




o 

a — 

*s> 

X 


O 

© 


© 

© 

© 

O 


© 

CO 

r-* 

O 


in 

r- 

co 


© 

o 

© 


T7 — 

co n 

to 

in o> 


in r> 

© © 

0) CM 
CM — 


tD ^ — 

© © 

— n — 

CD 10 CD 

O O o 


n 

co 


o 


o 
(0 
«* 

5 o in 

o — — 


© 

o 

in 


cn 

n 

a 


© 

CM 

o 

© 


to 

to 


CO 

to 



© 

© 

© 

*7 

© 

r*> 

o 


o 

© 


Q 

© 

© 

© 

© 

© 

© 


O 

o 

© 

CM 

— 

© 

© 

»- 


© 

© 

CM 

— 

— 

o 

— 

— 




o 

o 

o 

o 

o 

o 

o 

o 

CM 


CM 

© 

© 

r- 


© 

© 

CM 

CM 

© 

n 

© 

cr 

o 


© 

© 

© 


© 

CM 

5 

o 

© 

© 

© 

O 

O 

CM 

CM 

— 

O 

o 

o 





o 

o 

o 

o 

o 

o 

o 

o 

© 

o 

© 


is 

CM 

© 

© 

— 

CM 

© 


v 

© 

CM 

© 

© 

© 

© 

© 

CM 

© 

© 

© 

© 


© 

© 

CM 

CM 

o 

© 

o 

— 

o 

o 




o 

d 

P 

o 

d 

o 

o 

d 

o 

CM 

CM 

o 

© 

o 

© 

p 

© 

— 

© 

5 

© 

© 

© 

5 


© 

© 

m 

© 

o 

*7 

© 

n 

© 

© 

© 

© 


© 

© 

CM 

— 


— 

— 

— 

*- 



o 

o 

d 

o 

d 

o 

o 

O 


0} to 

0; - 

K o 

in tj 


O 

to 

to 

© 


o 

to 

o 


to — 

r- io 

o © 

— © 


CD 

O 0) 

r- ip 

CM — 


U> 

CM 

V 

<0 


CM 

(0 


CM 

© 

© 

© 

o 


© 

© 


© 

•e- 

o 


c- © © 

r- © to 

© — r** 

0-0 

CM CM CM 


CM 

© 


© 


tn h- O 

© — to 

in © - 

© © © 


o 

n 

© 

© 


© 

© 

© 


© 

CM 

© 

© 


© © r- © o 

O © rj a> n in 

CM © CM O CM © 

^ CM — 



© 

© 

* 

© 

n ai oi 


© 

© 


© 

CM 

© 

© 

© < 

D r~ 

© 

© 

© 

O 

© 

CM 

o 

N 

v © n 

© 

© 

© 

* 

© 

© 

© 


© © ® 

© 

o 

© 

o 

CM 

CM 

CM 

CM 

CM CM © 

* 


© 


o 

O 

d 

o 

odd 

o 

o 

o 


8 

o 

8 

8 

§ c 

3 s 

p 


o 


© 

6 

© 

o 

© c 

r> U 

? © 

8 

o 

© 

o 

© 


CM 

o 

h* 

© 

© c 

? N 

© 

© 


o 


*3' 

© 

© 

© © © 

© 

© 












ot 

d 

d 

d 

o 

o o o 

o 

o 

O 


ffi CM ® O O © © 

© r* — ® ^ n © 

© © cm © tj © S 

w N m n f) •• (ff 


8 

8 

o 

8 

o 



O 

© 

© 

© 

§ 

u 

© 

r- 

O 

8 

in 

O 

O 

© 

O 

o 

o 

o 

o 

O 

b 

d 

o 


118 

w) 


rms 


TABLE V (Continued) 



TABLE V (Continued) 



n 

© n 

o 

X 

cM 

CD 

a i 

D C 

0 CP 

CM 

IS 

rj 

IP 



u> 

CD X 

•— 

n 

X 

IP 

I 

r> c 

0 f- 

O 

in 

O 

CM 


IP 

in 

X 1 

V 

u> 

D 

r- 

o 

D ! 

P u 

■> IP 

X 

in 

O 

<7> 


• 

Q 

9 ' 

PM 

o 

O 

o 

M7 

O i 

0 c 

? o 

o 

g 


g 

g 

CM 

o 

O ' 

D 

o 

O 

o 

o 

& 1 

5 c 

5 o 

o 

5 

o 

6 

5 


o 

o i 

D 

o 

o 

o 

o 

o i 

D O O 

o 

o 

d 

o 

o 


IP 

X 1 

X 

o 

IP 

r- 

Ml 

O 1 

P U 

w 

n 

IP 

r** 

X 

x 


CM 

X X 

o 

r> 

CD 

o 

GD 1 

7) C 

M — 

•• 


O 

MI 

in 

O 

n 

r> n 

IP 

cm 

D 

r* 

CM | 

P C 

u r> 


IP 

CM 

m 



o 

Q < 


g 

g 

•- 

g 

CM | 


? Q 

— - 

g 


O 

Q 

CM 

o 

O i 

D 

s 

o 

o 

s 

O ( 

5 C 

5 5 

o 

5 

o 

o 

6 


o 

o i 

D 

o 

o 

d 

d 

o i 

D O O 

d 

o 

d 

o 

o 


n 

CO 

P 


IP 

MI 

CM 


CD (P 

IP 

CD 


r** 

m 


r> 

CP 

L> 


MI 

X 

O 

mi 

iP O 

M7 

X 

n 

(p 

o 

X 

CD 

MI 

cm 

o 

•• 

n 

CD 

MI 

tn cm 

<*— 


V 

x 



g 

5 

L> 

o 

g 

g 

Q 

o 

O O 

g 

g 

g 

o 

Q 


6 

O 

5 

o 

o 

o 

O 

O 1 

D O O 

5 

o 

6 

o 

6 


o 

o < 

o 

d 

d 

d 

6 

o i 

6 6 o 

o 

d 

o 

6 

o 




OI 

— 

r- 

0) 

O 

it; 

o 

— 

IP 


CD 

MI 

o 

X 

CM 

IP 


IP 

IP 

CM 

MI 

CD 

g 

CM 

IP 

CO 

X 


CP 

X 

MI 


o 

X 

o 

CM 

0) 

CM 

CM 

f- 

5 

C 

CD 

o> 

r- 


0) 

X 

X 

IP 

X 

X 


CM 


CM 

CM 


— - 

•- 

g 

g 

CM 


Q 

g 

Q 

o 

Q 

Q 

o - 

o 

o 

o 

o 

O 

o 

o 

o 

5 

o 

o 

5 

5 

o 

o 

5 

5 

K 

o 

o 

d 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

d 

d 

o 


X 

o 

•- 


X 

X 

X 

CM 

CM 

X 

X 

«p- 

X 


C 


CM 

CM 


X 


r- 

CM 

X 

n 

MI 

X 

CM 

X 

X 

n 

X 


— - 

o 

MJ 

X 


x 

X 

CM 

03 

X 

CM 

o 

MI 

n 


X 

X 

X 

Ml 

X 

T- 


g 

•- 


g 

g 

o 

g 

g 

Q 


g 

g 

CM 

O 

o 

O 

o 

o 

o 

o 

o 

o 

o 

5 

o 

o 

o 

o 

o 

o 

o 

d 

o 

o 

o 

o 

d 

d 

o 

d 

o 

d 

6 

d 

d 

o 

o 



n 


CM 

— 

MI 

n 

X 

X 

CM 

r- 

o 


a 

CM 

X 

CJ 

X 


n 

X 

X 

X 

X 

X 

X 

X 

CM 

X 

X 

X 

o 

CM 

X 

X 

X 

O 

n 

•— 

X 

X 

X 

MI 

n 


MI 

*p— 

X 

X 

o 

X 

CM 

CM 



■*- 



g 

g 

n 

X 

X 

r- 

**• 

•— 



o 

o 

g 

CM 

6 

o 

o 

o 

o 

o 

q 

o 

o 

O 

o 

o 

o 

o 

q 

o 

q 

o 


o 

o 

d 

o 

c 

o 

o 

o 

o 

o 

d 

o 

o 

o 

O 

o 

o 



o 

g 

g 

O 

C Q 

g 


o 

Q 

o 

g 

<5 5 

5 


X 

g 

X 

5 

x g 

X 


CM 

o 

f- 

X 

CM o 

r- 

o 

V 

MI 

n 

r> 

n n 

CM 

X 

6 

o 

6 

o 

d O 

o 



O 1 

o < 

3 § 

o 

o 

o 

o 

X 1 

D X 

o 

X 

CM ( 

D 

X 

CM 


~ O 

o 

q 

d < 

6 6 

o 

o 


120 


0.00990 0.00403 0.00203 



TABLE V (Continued) 


© O O t* h- .. .. .. w 

ajpjr-o«i» 6 ioovSiHVSic~S 

• Oooooopoo5oog-pf«r> 
•MOOOOOOOOOOOOOOOOO 


— g m tn o o ■to 

id O m o v cv in 

unO'O'-*- 
o g o o o o o 


n in n a> 

<(»#•■ 
— in O) 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

n 

cm 

© 

0 

_ 

© 

O 

© 

CM 

n 

© 

© 

r> 

© 

O 

K 

0 

<7 

in 

rr 

O 


n 

CM 

^7 

© 

O 

© 


© 

in 


n 

*7 


*7 

n 



•*- 

g 

O 

0 


W— 

f- 


V— 

0 

O 

Q 

6 

p 

O 

p 

0 

0 

g 

p 

a 




CM 

O 

0 

5 

O 

5 

O 

0 

O 

0 

O 

0 

0 

0 

0 

p 

0 

0 

d 

d 

0 

d 

d 

0 

d 

p 

0 

d 

0 

0 

d 

d 

0 

*T 

« 

2 

© 

10 


© 

0 

© 



© 

© 


0 ) 

© 

O 

in 

5 

© 

© 

n 

© 

© 

© 

0 ) 

© 

0 


© 

0 

© 

n 

CM 

CM 

•- 

CM 

CM 


CM 

O 

g 


n 


© 

O 

© 


0 

p 

O 

O 

0 

0 

p 

P 

g 

g 

p 

Q 

g 

0 

O 

0 

O 

O 

0 

0 

O 

5 

5 

0 

<5 

0 

0 

0 

5 

0 

5 

d 

0 

d 

d 

d 

0 

0 

d 

d 

d 

0 

0 

d 

d 

d 

d 

© 

© 


n 

© 

•C 

© 


h* 







© 

© 



CM 

© 

© 


n 

in 


© 

O 


© 

r> 


0 

O 

© 

© 

a > 

f- 

© 

c- 



g 

© 

f- 

r** 

© 

© 

O 



g 

g 

g 

O 

g 

g 

g 


Q 

O 

g 

O 

O 

0 

0 

0 

0 

0 

O 

0 

5 

5 

0 

O 

5 

5 

5 

O 

P 

d 

d 

0 

d 

d 

d 

d 

0 

d 

d 

d 

d 

d 

6 

d 

d 

© 


T 7 

© 

« 

© 

© 

r _ 

O 

in 

© 

© 

O 

© 

r> 




© 

© 


O 



^7 


in 

© 

© 

^7 

© 

0 

fx 

© 

Q 



© 

© 

*7 

’ey 

© 


0 

0 

f-* 

0 

© 


Q 


g 

g 

g 

g 

Q 


g 

g 



g 



0 

0 

O 

5 

5 

0 

0 

P 

O 

5 

0 

0 

0 

0 

0 

0 

d 

0 

d 

d 

d 

d 

d 

0 

d 

0 

d 

d 

d 

d 

d 

d 

CM 

r- 

© 

0 

cr 

© 

© 

© 


in 



r> 

CM 

CM 

O 

© 

n 

•- 

0 ) 

© 


© 

*— 

© 

CM 

© 

© 

n 

© 

O 

© 

© 

t*. 

*7 

V 

CM 


© 


O 

•— 

© 


© 

® 

O 

© 

O 

© 

O 

g 

•— 

g 

CM 

*— 


in 

in 

w- 



O 


0 

5 

0 

5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

P 

0 

0 

d 

d 

d 

d 

d 

d 

d 

d 

d 

6 

d 

d 

d 

d 

d 

d 


88888888 S 88800 B 8 

inginotnginoininoftoinotn 

wOf^inwot^incMr-inwof-incN 

OQ^ nrinncMCMni^-^*-*— 0000 

a o 0 0 oooood good dodo 


121 


TABLE V (Continued) 




v'w'/u 


TIME-WEAN AND TURBULENCE DATA FOR NONSWI RUNG 
FLOW $ = 0°, D/d = 2, WITH CONTRACTION 


l 


a 


I 




O) 


in 

o 

to 


b 

0) 

» 

ID 

n 

CM 

CM 

57 

© 

© 

57 




o 

in 

© 

b 

CM 

ID 

b 

57 

b 

0) 

CM 

© 

b 

P) 

57 

b 

© 



o 

in 


© 


5J 

to 

D 

to 

in 

CM 

57 

o 

© 

© 

© 

© 

57 




CO 


5j 

o 

5? 

O 

to 

b 

© 

c 

© 

57 

b 

n 

57 

to 

© 



cm 

o 



CM 

CM 

D 

© 

V 

0 

© 

© 

b 

b 

m 

a 

© 

© 




o 

o 

o 

o 

o 

o 

d 

d 

o 

d 

d 

o 

o 

o 

o 

o 

d 




o 

© 

0) 

01 

b 

57 

CM 

in 

57 

© 

57 

b 

o 


b 

© 

to 






rj 

b 

© 

CM 

0) 

o 

CM 

CM 

to 

D 

to 

b 

CM 

o 

to 



in 


o 



© 

in 

57 

o 

to 

b 

to 

© 

to 

© 

© 

to 

to 





CD 

CO 

*- 

© 

in 

in 

b 

03 

b 

in 

o 

to 

57 

© 

to 

to 





o 


CM 

CM 

© 


in 

ID 

b 

© 

0) 

0> 

0) 

© 

m 

© 




o 

d 

o 

o 

o 

> o 

o 

o 

d 

i o 

o 

o 

o 

d 

o 

6 

o 




© 

in 


0) 

© © 

to 

ID 

_ 

to 

b 


© 

© 

© 

K 

© 




to 

in 

CM 

in 

n © 

in 

to 

b 

to 

to 

CM 

in 

© 

b 

b- 

b* 



o 


CM 

© 

CM 

CM ID 

in 

b 

c 

i ~ 

b 

57 

CM 

CM 

CM 

CM 

© 




0) 

r- 



to m 

b 

o 

c 

► © 

5T 

© 

© 

© 

© 

© 

© 

c 

o 

o - 

o 

o 


CM 

cm n 

57 

ID 

b 

b 

© 

© 

© 

m 

© 

© 

© 

D 

_l 

X 

o 

o 

o 

d 

c 

> o 

o 

d 

c 

> d 

o 

o 

o 

6 

o 

o 

o 

3 

h 




















< 






















o 

in 

CM 

ID 

r- 

ID 

57 

51 

► -r- 

© 

© 

* 

— 

to 

to 


1 fU 

LU 


b* 

CO 

b 

CO 

in id 

0) 

o 

ID ID 

© 

in 

© 

to 

b 

h- 

© 



in 

0) 

a 

b 

b 

0) ID 

© 

o 

a 

) to 

n 

O 

o 

O 

CM 

© 

© 


M 



in 

in 

in 

cm in 

b 

© 

ID 03 

o 




CM 

O 



Nl 

d 

o 

o 

o 

o 

CM 

<7 

ID 

a 

) © 

to 

to 

© 

© 

© 

© 

© 


o 
' z 


6 

» 

6 

i 

6 

i 

6 

i 

d o 

o 

d 

c 

> o 

6 

6 

o 

o 

o 

o . 

o 




o 

r> 

w 

in 

© 

in 

57 

ft © 

CM 

to 

o 

to 

© 

5J 

b- 




n 

Q 

© 

ID 

CM CD 

0) 

© 

a 

> © 

© 

O 

in 

CM 

b 

© 

O 



o 


O 

to 


b 

© 


57 

ID b 

b 

© 

© 


CM 

© 

57 




X 

Q 

Q 

o 

o o 


o 

O - 



CM 

© 

to 

57 

© 



o 

5 

5 

5 

o 

o o 

o 

01 

a> a 

© 

© 

© 

© 

© 

© 

© 




o 

d 

6 

i 

0 

1 

0 o 

1 i 

d 

o 

o o 

o 

o 

d 

d 

o 

e 

o 




n 

o 

o 

o 

r 

3 S3 

o 

Q 


? o 

o 

o 

o 

o 

o 

o 

9 




fi 

o 

5 

5 

c 

3 O 

o 

o 

c 

5 o 

O 

o 

O 

o 

y 

o 

o 




Q 

in 

o 

in 

c 

5 in 

o 

in 

t 

? in 

o 

© 

y 

© 

5 

© 

o 




Q 

r- 

in 

cm 

3 

S b- 

in 

CM 


5 b 

to 

CM 

o 

b 

© 

CM 

6 



o 

tr 

© 

n 

© 

© CM 

CM 

CM 

CM — 

■*" 

5“ 


o 

o 

o 

O 



8 

a 

o 

o 

6 

6 

d o 

d 

o 

o d 

d 

o 

o 

o 

d 

d 

o 



123 



TABLE VI (Continued) 






Of 


«- 

o 

© 

CP 


o 

® 

© 

Of 

o 

Of 

© 

o 



tn 

o 

Of 

0 ) 

in 

o 

n 

Of 

f- 

0 ) 

o- 

0 * 

r> 

*“ 

IP 

© 

03 

o 

o 


9 

9 

ID 

© 

© 

Of 

Of 

® 

© 

Of 

9 

0 * 


9 

Of 



Of 

Of 

0 * 

® 

o 

© 

o 

r* 

© 

in 

o 

Of 

0 ) 

9 

9 

tn 

Of 

•*- 





Of 


Of 

*“ 


Of 

Of 

Of 


Of 

Of 

Of 


6 

o 

d 

o 

o 

o 

o 

d 

O 

o 

o 

o 

o 

d 

o 

o 

o 



CD 


© 

in 

9 


® 

ID 

0 ) 

O 



o 

9 

© 



T- 

© 

o 

ID 

CD 

h- 

0 - 

9 

Of 

© 

9 

Of 

r- 

o 

9 

9 

h* 

tn 

© 

© 

o 

O 



r* 

r- 

® 

in 

9 

Of 

tn 

© 

© 

IP 

© 


Of 

9 

tn 

CD 

ID 


r- 

CD 

9 

a 

0 ) 

O 

® 

© 

© 


03 

— 






Of 

— 


Of 



Of 

«•» 

•** 


Of 

Of 


j® 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

d 


o 

Of 

(71 

9 

CO 

© 

Of 


Of 

® 

*9 

in 

© 

9 


Of 

o 




w- 

o 

CD 

rx 

® 

■9 

o 

in 

9- 

© 

O 

O- 

r~- 


9 

o 

of 



© 

CO 

in 

o 


in 

0 * 

© 

© 

© 

r- 

h* 

© 

9 


Of 

9 

© 

9 

o» 

o- 


Of 

tn 

in 

03 

© 

© 

© 

© 

© 

© 







"" 

Of 

Of 

of 

Of 

Of 








o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

d 

o 

o 

o 

o 

o 

o 



CD 

o 

Of 

© 

03 

Of 

of 

o 

03 

© 



© 

__ 

0) 

© 



03 

© 

T— 

in 

r- 

t- 

0 - 

in 

CD 

© 

© 

D 

© 


o 

© 

tn 

tn 

o 

CD 

cn 

© 

© 



ip 

0 ) 

O 


0 * 

f- 

© 

n 

0 ) 



tn 

in 

CD 


9 

Of 

© 

of 

Of 

O 

CD 



r- 

Of 

© 

6 

o 

o 

o 

o 

*■* 


Of 

Of 

Of 

Of 

Of 


Of 

Of 


Of 



d 

i 

o 

» 

o 

t 

0 

1 

d 

. d 

d 

d 

o 

o 

o 

o 

d 

d 

o 

o . 

o 


o 

0 


CD 


0 ) 

CO 

O) 

o 

© 

h* 

© 

o 

o* 

o 

03 

© 


o 

s 


Of 



ID 

r- 

03 

9 

Of 

03 


9 

© 

9 

Of 

o 


O 

o 

CO 

Of 

o 

® 

© 

Of 

<n 

0 * 

o 

Of 


9 

© 

0 - 


o 

Q 

9 

© 

© 

9 

03 

Of 

O 

Of 

0 - 

© 


0 ) 

O 

O) 

03 

O 

6 

5 

o 

o 

o 

o 

o 




o 

o 




o 


! 

o 

o 

d 

t 

o 

d 

i 

o 

d 

o 

o 

o 

o 

o 

o 

d 

d 

d 

o 


o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 


o 

o 

Q 

Q 


o 

o 

Q 

o 

o 

o 

o 

o 

o 

o 

o 

o 

Q 

o 

g 

o 

Q 


o 

in 

5 

in 

o 

in 

o 

tn 

o 

in 

o 

IP 

o 

in 

O 

IP 

O 


o 


tn 

Of 

o 

o- 

in 

Of 

o 

r~ 

in 

Of 

u 

o- 

© 

Of 

o 

Q 

9 

© 

© 

03 

03 

Of 

Of 

N 

Of 

*“ 


*" 

*" 

o 

c 

o 

O 

CL 

6 

o 

6 

o 

o 

o 

d 

d 

o 

o 

o 

o 

c 

o 

o 

d 

o 


124 



(b) v/u 


TABLE VI (Continued) 


o 

« 


m 


o in f- co tP 2 

<x o> in "a J? 

in « to w ip 2 


v in 
«o o> 
n t» 
o 


in 

CM 

M 

oi 


— in 

r- n 

ps p> 

— r- 


so ® ^ 2 2 S 

tp m o in ® 

*- p> ••- r- in w 

h in ^ - J 5 

— — — — 00 


O ' o 


o ^ ti in o> O jp 

pj 0) k o n in }P 

® ^ ® S ffi m o 

«- «J PI Cl ID ® 0 

^ ^ n 


in t* rt ~ 

co © in in 

® cm ® 

*- *- 0 ) in 

cm n — 


r- 

in 

cm 


ip 

n 

n 

ip 

n 


n id f- cm 

ft o « u> 

0 0)00 

^ o o o 

0000 


o o 


o 

o - 

**v 

X 


co ■«•* p> 

Pi — « 

r- m w 

0) «- « 

o - - 


w 

o 

in 

in 


ai 

o 

in 

in 


n 

01 

o 

CO 


1 - O cm 

r- cn cn 

r-* ip cm 

CO O f* 


n 

o 

<n 

O 


in r- 

n in 

f- o 

o 

o o 


o 

o 

o> 

D 

o 


o 

o 

« 

n 

o 


0 ) ® ® 

in ^ o 

0 ) 01 <T) 

ft ft ft 

OOO 


*- ai in oi *■ S 

tj r- ® r- ^ co 

to 0) - r- £ 2 

o ® O ^ 

O O O O - - 


CM 

o 


cm 

n ® 

r> <? '- 

to co 

— *- o 


co 

© 

n 

O 


- o 


CD 

CM 

o 


«M 

o 


v n ip 

— r- ip 

SO IP IP 


- o 

in <0 

in in 


in 

E 

u 


O 

6 


o 

6 


o 

o 


o 

o 


o 

d 


3 

"o 


o 

o 

o 


8 


01 - O IP 

t» m 01 2 

M CM 01 

CM SM « T* 

OOOO 


o in co 

— m n 

in in ip 

OOO 


co 


— o 

0) 01 

rj in 


w to in oi to 

r n o It ' 

in v sj mj tj 


0 

0 

0 

0 


Q 

0 

Q 

g 

5 

0 

0 

5 

in 

0 

ui 

0 

c* 

in 

PS 

n 

n 

Cl 

pi 

b 

b 

d 

d 


O O 

o o 


o 

o 


o 

o 


o 

o 


o 

o 


o 

o 


o 

o 


8 

in 


cm 

o 


in 

cm 


O 

Q 

0 

g 

O 

8 8 

2 

O 

in 

O 

0 

in 

O 

0 

in 

0 0 

0 m 

0 

0 

Of 

O 

r** 

in 

CM 

0 t; 

in 

M 

CM 


** 


— 0 

O 

b 

b 

b 

b 

b 

b b 

b 


in 

« 

O 




125 


£ 


TABLE V! (Continued) 


o 

n 


c 

D 



0) 

cn 

D 

D 

ct 

CD 

CO 

r- 

in 

10 

CD 

CT 

05 

in 

o 

*7 

— 

n 

cn 

CD 10 

0) 

0) 

h* 

CN 

f- 

ct 

O) 

in 

O) CN 

® 

0) 

m 

10 

o 

© 

«7 


CO 

h* 

10 

in 


o> 

CJ 

c 

CT 

in 


© 

in 

to 

u> 


03 

o* 

CD 

CD 


o 

to 

o 


r- 

in 

P- 

tn 

o 

o 

o 

0 

o 

o 

O 

o 

o 


o 


c 

o 

o 

o 

o 

o 

o 

d 

0 

o 

d 

o 

d 

o 

o 

o 

o 

c 

o 

o 

d 

o 


IP 

O 

V 

1C 

0* 

CN 

CD 

CN 

in 

L 

® 

CN -<7 

in 

03 

Nf 


rt 



CN 

in 

W 

in 

n 

in 

CD 

CT 

r* 

03 

to 

CN 


n 

oi 

h* 

o 

c 

ID 

o 

<*- 

01 

h* 

in 

o 

cn in 

CT 

0) 

t7 

in 

CP 

in 


C E 

© 

0) 

c- 

in 

01 

in 

N 

0) 10 

tp 

CT 

CN 

o 

o 

o 

0 

c 

o 

o 

o 

o 

o 

o 

o 

c 

o 

o 

o 

O 

o 

d 

o 

D 

c 

d 

o 

o 

o 

o 

d 

o 

c 

o 

o 

d 

d 

0) 


o 

n 



to 

n 

10 

03 

^7 

CT 

0) CN 

n 

0) 


in 

<N 

n 

D 

* 

in 

V 

w~ 

ct 

h* 

CO 

01 

CT CN 

ID 



O) 

*7 

CD 

D 


CN 

CD 





CT f" 

h* 

f* 

r- 




i- 

h 

03 

CD 

in 

10 

03 

n 

CT 

r 

CN 

CN 

CN 

CN 

o 

o 

o 

o 

c 

o 

o 

o 

o 

o 

o 

u 

c 

o 

o 

o 

o 

o 

d 

o 

o 

o o 

o 

o 

d 

o 

6 

6 

c 

o 

o 

o 

o 


CN 

CN 

9 

— If) 

O' 

to 


o 

0) 

O 


to 

o 

CT 

_ 

f- 

o> 

00 

10 

in 03 

<N 

n 

0> 

G 

in 

CT 

c 

in 

to 

CT 

to 

0) 

CT 

ct 

CT 


o 



01 

o 

<7 


c 

r- 

® 

r> 

CT 

CN 

CN 

n 

ct 


in 

10 

in 

o 

CN 

CN 


c 

1 o 


o 


o 

o 

o 

o 

o o 

o 

o 

o 

O 

o 

o 

c 

o 

o 

o 

o 

o 

d 

o 

o 

d d 

o 

d 

o 

o 

d 

o 

c 

o 

d 

d 

o 

o 


in 

in 


O 

o 

CD 


r* 

o 

01 

o> n 

to 

CN 

P^ 

5 

?5 

to 

01 



CO 

in 

n 

m 

N 

in 

c 

\ 01 

f- 

CT 

03 

Q 



w 

O 0) 

o 

r* 

in 


r? 

CN 

r- o 

O) 



s 

o 

cn 

04 

CN — 

CN 


CN 


CN 

CN 


■*- 

CN 


o 

o 

o 

o 

o o 

o 

o 

o 

o 

o 

o 

o o 

o 

o 

o 

o 

d 

o 

o 

o o 

o 

d 

o 

d 

d 

o 

d o 

o 

d 

o 


o 

o 

O 

G O 

o 

O 

o 

G 

o 

o 

o o 

o 

O 

p 

5 

5 

6 

5 

C 

? o 

ft 

ft 

Q 

o 

o 

G 

c 

> o 

u 

o 

O 

Q 

in 

o 

in 

f 

> in 

o 

in 

O 

in 

o 

in 

c 

5 tn 

o 

tn 

Q 

Q 

r- 

in 

CN 

f 

5 r* 

in 

CN 

O 


in 

CN 

t. 

5 r* 

in 

CN 

Q 


o 

ct 

n 

n cv 

CN 

CN 

CN 

■** 

•* 


- o 

O 

O 


o 

o 

o 

o 

o d 

o 

o 

d 

o 

o 

o 

o d 

o 

O 

o 


i 


> 

‘‘ 5 * 


t 


(26 


rms 


TABLE VI (Continued) 




o 


tn 


O 





CM 


10 


* 

© 

o 


CO 

© 


© 

CM 

OJ 


© 



o 


IP 

© 


© 

© 

© 

•- 

r> 



CM 

rf 

CM 


o 









o 

o 

o 

o 

o 

o 

o 

o 

o 



o 

o 


IP 

CM 

!*• 

p* 



CD 



D 

IP 

© 

© 

$ 

te 

o 

n 

O 

o 

CM 

P 3 

O 

CM 

O 

r> 

P* 

CM 

© 

CM 

o 

o 

d 

o 

o 

o 

o 

o 

d 




o 

© 



CM 

o 




n 

o 

P* 

© 

© 

© 



o 

00 

p- 

© 

P) 

© 

© 


p* 

00 

© 

0) 

CM 

CM 

o 

© 

o 

o 

o 

o 

o 




o 

o 

0 

o 

o 

d 

o 

o 

o 

o 



IP 


5 T 

© 

© 

CM 



CM 


IP 

CM 


r* 


IP 

n n 

r- 

to 

© 

P* 

r* 

0) 

© 

O 

© 

© 

© 

o 

c 

o 

o 

o 

q 


O 

O 

o 

o o 

o 

d 

o 

o 

o 

o 

o 

c 

> h* 

IP 

IP 

o 

o 




? 

> © 

CM 

© 

(0 

in 

© 


5 

r 



LO 

01 

ip 

P* 

© 



<0 

r> 

© 

r> 

o 

o 

§ 

O o 

o 

o 

o 

o 

o 

o 

o 

o o 

o 

o 

d 

o 

d 

d 


Q 

o 

o 

O o 

O 

Q 

8 8 

5 

5 

o 

G Q 

o 

Q 

O o 

5 

© 

o 

© Q 

© 

o 

£ 8 

X 

!»► 

© 

CM O 

r- 

© 

CM G 

n 

« 

© 

© © 

CM 

CM 

CM CM 

6 

6 

G 

6 o 

6 

o 

o d 


TT 

p> 

h 

p> 

© 

*3 

© 

© 

© 

© 

© 



© 

CM 


© 

© 

O 


p- 

© 


© 

p> 

o 

o 

© 

© 

© 

55 

© 




O 

o 

O 

G 

o 

G 

o 

o 

d 

o 

o 

o 

d 


© 

o 

o 

*» 

p- 

© 

© 

© 

CM 

w 

© 

CM 

© 

o 

o 

© 

© 

r~ 

© 

© 


© 

0) 

© 


© 

© 

© 

© 


m 

© 

— 

O 

o 

o 

o 

o 

o 

o 

o 

o 

d 

o 

o 

o 

o 

d 



CM 


© 

© 

p- 

p> 

© 

© 

© 

CM 

CM 

© 

© 

© 

r- 

o 

© 

o 

© 

*- 

© 

© 

© 

r- 

M 3 

r> 

CM 

CM 




O 

O 

O 

o 

o 

o 

o 

O 

o 

o 

o 

d 

o 

o 

o 

o 


© 

o 

© 

© 

© 

© 

© 

© 

n 

© 

© 

p> 

CM 

© 


© 

© 

p* 

CM 

■r* 

O 

© 

© 

© 

CM 





Q 

o 

g 

O 

o 

o 

o 

o 

O 

o 

o 

d 

o 

d 

o 

d 

o 

d 

o 


© 

© 

r> 

© 


© 


CM 

n 

© 

© 

© 

G 


CM 

© 

to 

© 

© 

© 

© 

p- 

© 

P* 

o 

o 

o 

o 

Q 

o 

o 

o 

5 

5 

o 

o 

O 

o 

o 

O 

o 

d 

o 

o 

o 

o 

o 

o 


o 

0 

o 

Q 

o 

o 

o 

o 

5 

Q 

o 

5 

o 

o 

o 

o 

© 

Q 

© 

o 

u> 

o 

© 

0 

f- 

© 

CM 

c 

t~ 

in 

CM 

o 





o 

O 

o 

o 

o 

o 

o 

o 

6 

o 

o 

d 



2 

« 


i 

I 

i 

: 

i 

!■ 


127 



rms 


TABLE VI (Continued) 



’v’/u 


TIME-MEAN AND TURBULENCE DATA FOR SWIRLING 



129 


TABLE VII (Continued 



0.00000 0.24016 -0.44969 -0.6B424 -0.49251 -O.S736S 


TABLE VI I (Continued) 



(c) w/u 


TABLE VI! (Continued) 


o - 


I tn 


n 

ro 

01 

ro 

o 


57 

in 

oo 

57 


57 

00 

CM 


O 

o 


8 

O 

o 

v 


o 

ip 

ro 

•w- 

o 

ro 

o 

CM 

o 

IP 

* 

in 

IP 

in 

r- 



*J 


c 

o 

o 

O 

o 

O 

d 

o 

6 

o 

o 

ro 

ip 

in 



ro 

in 

0) 

0) 


CM 

ro 

57 

5J 

h* 

h- 

T*~ 

r- 

r- 


O 

o 

O 

O 

O 

o 

o 

6 

o 

o 

CM 

r* 

0) 

CM 

<71 


01 

ip 


in 

*7 

in 

in 

r- 

CM 

Cl 

ro 

CM 


CM 

■** 

— 

** 


* — 

d 

o 

o 

o 

o 

CM 

0) 

r* 

ro 

CM 

o 

ro 

O 

in 

o 

*7 

ro 

^7 

ro 

01 

ro 

ro 

ro 

r~ 

in 

CM 

CM 


*■“ 

**■ 

o 

O 

o 

o 

o 

ro 



ro 

ro 

01 

0) 

O 

o> 

01 

n 

ro 

ro 

*" 


n 

CM 


*— 




*- 

** 

T" 

6 

o 

o 

d 

o 

o 

o 

o 

c 

g 

o 

o 

o 

c 

o 

in 

o 

in 

c 

in 

r- 

in 

CM 

c 

r- 

n 

ro 

ro 

r 

CM 

o 

o 

o 

C 

d 


ffi 

o 

CM 


'J 

cm 

ro 

in 

O 


r- 

CO 

o 


eo 

CO 

O) 


r~ 


ip 

C- 


5 f 

f- 

CX) 

CD 

O 


o 

IP 


Cl 

CM 


M 


Cl 

O 


O) — 

co *- 

cm in 

Oi O 

° r 

d o 


r- oi ci 

r- - cj 

IP d '«fl 

CM CM O 


0 

n 

01 

o 


Ip 

CM 

o 

IP 


0 

in 

01 


<£> 

tP 

IP 


n 

r- 

«p 

oo 


CM 

n 


® *7 < 7 ) 

O cd r- 

O O) IP 

n v in 


room 

in co c* 

C) f- CM 

room 

*- CM CM 


ro in 

ci 

co ro 

— o 

CM CM 


ro 


m 

CM 


<7 

CM 


in 

B 

L. 


IP 


in 

CD 

CM 


f- 

in 

0) 


O 

m 


ro 

r- 


8 

m 

CM 

CM 


0) 

ro 

ro 

CD 

*7 

ro 

in 

ro 

r- 


ip 

IP 

ro 

o 

01 

01 

in 

o> 

01 

m 

cm 

ro 

57 

Ol 



57 

ro 

ro 

o 

o 




— 

*■' 


CM 

CM 

CM 

o 

o 

d 

o 

o 

o 

O 

o 

r*- 

o 


CM 

ro 

n 

CO 

in 

T7 

r- 

«— 

o 

CM 

in 

in 

CM 

r- 

o 


CM 

01 

ro 

r- 

o 

t'- 

*57 

CM 

O! 

•* 

CM 

ro 

ro 



ro 



** 

*~ 


o 

d 

d 

d 

6 

o 

o 

d 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

in 

o 

in 

o 

in 

o 

in 

o 

t- 

in 

CM 

o 


in 

CM 

CM 



*" 

*• 

O 

o 

o 

o 

d 

o 

d 

o 

o 

o 

d 


*7 

CM 

in 


(P 

O 

n 


132 


TABLE VII (Continued) 






€N 

01 

Ol 

CO 

r> 


in 

in 

ID 

in 

ID 

® 

r> 

Ui 


^7 


r> 


in 

m 

<N 

o 


O 

n 


O 

® 

® 

'y 

ffl 

® 

O) 

o 

Cl 


■*- 

V 



*7 

O) 

ID 

in 

in 

CN 

n 

0) 

<7 

CC 

•*- 



i** 

CO 

to 

01 

o 

0) 

CT) 

r- 

O) 

o 


o 

r- 

o 


n 

CN 

o 

o 

o 

o 

c 


o 

o 

o 

o 



*" 


CN 

CN 

IN 


o 

o 

o 

d 

o 

o 

o 

d 

o 

o 

o 

o 

o 

o 

o 

O 

o 


r- 

r> 

CD 

OJ 

in 

r> 

D 

CO 

r* 


CO 

CN 

CN 

01 

® 

h* 

*7 


IX 

CN 


n 

o 

<7 

in 

c r. 

in 

ID 

CD 

o 


ID 

in 

r- 

r> 

in 


o 

w— 

cn 

in 

n 

ID 


CN 

ID 

o 

co 


o 


r> 

n 


m 

(7) 

0> 

o 


0) 

01 

o 


o 

CN 


r* 

® 




— 

o 

o 

o 



O 

o 



■*" 

*- 




CN 

CN 

CN 


o 

o 

o 

d 

o 

d 

o 

d 

d 

d 

o 

o 

o 

d 

o 

o 

o 


CD 

r- 

in 

f- 


o 


n 


CD 

in 


cN 


CD 

ID 

® 


IT 

«7 

ci 



CD 


ID 

r~ 

01 

r> 



® 

CN 


in 

o 

N 

o 

CO 

•*- 

o 

ID 

to 

in 

0) 

ID 

m 

in 

u> 

n 

CN 

n 

•* 


r 

-o 

n 


in 

n 

CN 

CN 

o 

^7 

in 

ID 

n 

r~ 

® 

o 

o 

o - 



*- 


•- 






*- 


CN 

cv 

CN 

r> 

r> 



■ 

♦ 

• 


• 

• 

» 




» 


► 


• 

• 

X 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

d 

o 

d 

o 

o 


O 

in 

o 

o 


o 

n 

in 

o 

n 

O) 

ID 

o 


CN 

0) 

« 


0) 

CO 




ID 

n 

■»— 



o 

ID 

o 

o 

in 

in 

in 

1 in 

o 

c* 

*- 

f- 

fcO 

c- 

*! 

ID 

CN 

•- 

n 

r- 

® 

'I- 

n 

in 

■c 


in 

ro 

CO 


(N 

CN 

n 

r> 



ID 

ID 

® 



01 

o 

d 



— 

"* 


— 


*“ 



** 



** 


*• 

CN 


o 

o 

d 

d 

d 

d 

o 

o 

d 

o 

o 

o 

o 

o 

o 

o 

d 


o 

o> 

n 


CO 


CO 

C- 

r- 

ID 

® 

ID 

ID 

o 

« 

t" 

CN 


c 

o 

10 

0) 

■»— 

in 

n 


r- 

n 


r- 

r- 

*7 

® 

01 

n 

o 

c 


01 

n 

r- 



n 

in 

0) 

CN 

n 

O) 

o 

® 

w 

r> 


c 

n 

rj 

n 

n 

CN 

o 

ID 

in 

*7 


•c 

o 

CN 

o 

o 


d 

o 

— 

** 

~~ 

<r * 



CN 








o 

•• 


c 

o 

d 

d 

o 

d 

d 

o 

o 

d 

d 

o 

o 

d 

d 

6 

o 


r 

o 

o. 

o 

Q 

o 

o 

o 


o 

0 

o 


o 

o 

Q 

o 


2 

o 

o 

o 

o 

o 

o 

o 

Q 

o 

Q 

o 

Q 

o 

o 

O 

P 


2 

in 

o 

m 

O 

m 

o 

in 

Q 

in 

o 

in 

o 

in 

o 

in 

o 


c 

r- 

in 

CN 

O 

r- 

in 

cn 

o 

i- 

in 

CN 

o 

N 

in 

CN 

o 

o 


n 

n 

D 

n 

in 

fN 

CN 

CN 

*- 




O 

o 

O 

o 

CC 

c 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 



> 

o> 





133 


rms 


TABLE VII (Continued) 




0.20938 0.2265B 0.24990 0.19947 0.24578 


TABLE V!l (Continued) 






0 

U> 

c 

o 

eg 


CD 

0) 

IP 

Ip . 

in 

CD 

a 

*7 

o 

0 


cr 

r* 

C 

r- 


0) 

t7 

CD 

DT 

CD 


CD 

D 

■— 

o 

0 

a 

o 

o 

o 

o 

o 


o 



eg 

CD 

CN 

in 

*7 



IP 

CD 


Q 

Q 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 


eg 

rg 

« 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

5 

o 

5 

o 

o 

o 

o 


o 

o 

d 

d 

o 

o 

o 

o 

o 

o 

o 

d 

d 

o 

o 

o 

d 


* 

<7 

f** 

CO 

m 

CD 

CD 

o 

r* 


a. 

CN 

o 

CD 

o 

IP 

< 

0 


CD 

a 

o 

o 



W 


in 

IP 

o 

o 

0 

CN 

CN 

IP 

0 [ 

ip 

•- 

•» 

CN 

eg 

id 

eg 

CD 

eg 

ip 

D 

r*» 

1- 

01 

o 

CD 


0 

^ 1 

* 

o 

o 

O 

Q 

o 

o 

o 

o 

o 

o 

o 

o 




CN j 

m 

5 

o 

o 

O 

o 

o 

o 

o 

o 

o 

S 

o 

5 

o 

o 

o 

° 


o 

d 

o 

o 

o 

o 

o 

o 

d 

d 

o 

o 

o 

o 

o 

o 

o 


If) 

CD 


IP 


<9 

0) 


X7 


CD 

0) 

0 

o 

o 

r- 

0 1 


l- 


CN 

CD 

CN 


o 

eg 

r* 


o 

0 

CD 

o 

0 

o 


o 

0) 

o 

CO 

0 

o 

Cl 

ID 

CD 

CD . 

CD 

CN 

D 

<7 

0 


in 

0 ; 


CJ 


O 

o 


o 

y 

o 

o 

o 




eg 

CN 

O 

o pg n 

o - 

o 

o 

o 

5 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

3 

X 

d 

o 

o 

d 

o 

o 

o 

d 

o 

o 

o 

d 

o 

o 

d 

o 

d 

> 



















> i 



















3 


*7 

a> 

CN 

0) 

*7 

n 

CD 

CP 

•- 

*7 

IP 

0 

r- 

eg 

eg 

w~ 

o 


cd 

in 

ID 

IP 

CD 

CD 

eg 

eg 

eg 

CD 

IP 

eg 

CN 

in 

CD 

o 

e* 1 

0 

»- 

o 

CD 

IP 

o 

0) 

0 

o 

r- 

O 

0 


CD 

0 

0 

0 

Q 

® I 

» 

cn 

ex 


Q 

"T- 

g 

O 


o 

CN 

Q 


O 

o 

O 

O 

CD 

o 

o 

o 

o 

o 

O 

5 

o 

o 

5 

o 

5 

o 

o 

o 

5 

5 

5 


o 

o 

d 

o 

o 

d 

d 

o 

o 

o 

d 

o 

o 

d 

d 

o 

o 



CD 

IN 

CD 

o 

o 

CD 

•*7 

0 


_ 


o 

^7 

^7 

CN 

1 


Q 

ID 



o 


■*- 

CD 

0 

«» 

CD 

0 

o 

NT 

o 

n 

IP 

o 

o 


CD 

^7 

u 

CN 

*7 

■»- 

CD 


IP 

v- 

g 

CD 

r* 


CD 

* 

o 

o 

o 

o 

o 

O 

o 

eg 

o 



r 

CJ 

o 


CD 

*— 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 


o 

d 

d 

o 

o 

o 

d 

d 

o 

o 

6 

6 

d 

d 

d 

o 

O 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 


o 

o 

o 

g 

o 



Q 

o 

o 

o 

o 

o 

o 

o 

o 

y 

o 

D 

CJ 

o 

o 

5 

o 

Q 


O 

IP 

o 

IP 

o 

IP 

IP 

5 

in 

o 

IP 

o 

IP 

o 

IP 

o 


o 

r- 

IP 

CN 

5 

f* 

IP 

CN 

o 

f- 

IP 

CN 

o 

r- 

in 

eg 

o 

a 

*7 

CD 

n 

CD 

ID 

CN 

CN 

eg 

CN 


w- 

«»- 


O 

o 

O 

o 


• 


* 





• 

• 

• 

• 






a 

o 

o' 

o 

d 

o 

d 

d 

o 

o 

o 

o 

d 

o 

o 

o 

o 

o 


135 


TABLE VII (Continued) 


i 



O 


in 

P) 

© 


CD 

© 

oi 

© 


o 


© 

CN 

© 



cm n 

n 

n 

n 


CD 

CN 

CN 

CN 

in 

© 

u 

© 


© 

o 


o 

o o 

O 

n 

o 

O 

O 


CN 

CN 

CN 

© 

© 


CN 

n: 

© 




o 

o 

ft 

o 

o 

o 

p 

o 

o 

o 

Q 

O 

•— 

© 



CN 

8 o 

5 

o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 



6 o 

o 

o 

o 

o 

d 

d 

d 

o 

o 

d 

o 

o 

* 

o 

o 

, 



n 


CO 

in 

CN 

CD 

0) 

© 

CD 

© 

o 


© 

© 

•7 





in 

in 

o 

in 

P" 


CN 

© 

© 


© 

© 

CN 

P- 







CN 

CN 


CN 

tN 

© 

r> 

© 

CN 



© 



a o 

o 

Q 

o 

O 

o 

o 

O 

o 

o 

o 

o 


T* 


2 


— 

o 8 

o 

5 

o 

o 

o 

o 

8 

o 

o 

o 

o 

o 

o 

o 

o 



d o 

d 

r 

o 

d 

d 

o 

o 

o 

d 

d 

d 

d 

o 

d 

o 




CN 

© 

P- 

n 

m 

o 

r** 

© 

© 

© 


© 

r- 

© 

© 




o 

IT 

r* 

in 

r** 

ID 


o 

cN 

© 

© 

© 

O 


© 




© 

© 

in 

in 

rv 

© 

© 

IT 

in 

© 

© 

© 

© 

© 

© 




o 

o 

o 

O 

p 

o 

o 

O 

o 

o 

U 

o 

*■* 


~ CN U 


8 8 

8 

8 

8 

8 

8 

o 

o 

o 

o 

o 

o 

G 

o 

o 

o 



o o 

o 

o 

o 

o 

d 

o 

6 

d 

d 

o 

d 

o 

o 

o 

o 



















s 


















3 





o 

<r 

CD 

o 

CN 

O 

© 

N 

© 

© 

CN 

CN 

© 




CN 

in 

<71 


0) 

CD 

© 

in 

© 

O 

© 

© 

© 


o 




in 

t- 

ID 

o 

CD 

in 

r- 

h* 

© 

r- 

«*» 

P- 

© 

** 

© 

s 



CN 

o 

Q 


g 

o 

o 

o 

u 

O 

■*— 

o 





o 

o o 

o 

o 

8 

o 

5 

o 

o 

O 

o 

o 

c 

o 

o 

o 

o 



d o 

o 

o 

d 

o 

d 

d 

o 

o 

o 

6 

d 

o 

d 

o 

d 



O © 

cn 

o 

NJ 

CD 

n 

in 

CD 

r- 

CN 


o 

r- 

CN 


© 

i 





CN 

CN 

© 

CD 

© 

CN 


© 


CN 


XT 

o 

< 

o 


ID 

o 

in 

© 

CN 

CN 


in 


CN 


© 

© 

© 

IN 

1 


fi o 

o 

o 

o 

0 

O 

CN 

w 

o 

© 


o 

G 


o 

o 


o 

8 8 

o 

5 

o 

8 

o 

o 

O 

5 

o 

o 

o 

u 

O 

o 

O 

1 


o o 

o 

d 

o 

d 

o 

d 

d 

d 

o 

o 

o 

d 

o 

o 

o 

i 


o o 

o 

o 

n 

o 

o 

o 


o 

o 

o 

9 

o 

o 

o 

9 




5 

o 

Q 

o 

9 

o 

ft 

o 

G 

o 

O 

o 

o 

o 

o 

> 



o 

U) 

Q 

in 

o 

in 

o 

in 

O 

IT 

g 

© 

o 

© 

o 

[ 



© 

CN 

o 

r* 

in 

CN 

o 

r*- 

in 

CN 

o 

t- 

© 

CN 

o 


a 

^ r? 

n 

n 

n> 

CN 

CN 

CN 

CN 





o 

o 

c 

o 

i 

a 

d d 

o 

o 

o 

o 

o 

d 

o 

o 

o 

d 

o 

o 

o 

o 

o 

i 


136 


E' 


TABLE VII (Continued) 



CD 


w 


o 

w 

h* 

eo 


r- 

o 

ID 

r> 

« 

K 

ID 

01 

n 

Q 

o 

Q 

O 

o 

CD 

o 

o 

O 

O 

o 

o 

o 

d 

O 

d 

6 

d 


0) 

o 

CD 

m 

*7 


in 

»- 



1- 

h- 

n 

8 

t- 

o 

o 

0) 

§ 

O 

01 

8 

0) 

cm 

O 

o 

6 

d 

o 

o 

o 


CO 

CP 

<0 

CM 

01 


r- 


0) 

CD 


n 

CM 

n 

o 

CD 

n 

o 






cm 

o 

O 

O 

o 

o 

O 

o 

o 

o 

o 

o 

o 



m 

n 

r> 

-- g 


n 

CD 

in 

O 

® 

o 

o 

O 0) 

o 8 

CO 

8 

oi o 

8 o 


o 

o 

o 

o 

o 

o 

in 

CM 

CM 

in 

in 

CN 

CD 

CM 

CO 


® 

CO 

n 

0) 

P5 


in 


T* 

Q 

O 

o 

O 

o 

o 

O 

o 

S 

o 

o 

o 

o 

o 

o 

d 

o 


o 

5 e 

£ 

o 

© e 

g 

o g 

© 

o 

o 5 

in 

in g 

in 

o 

tn o 

CM o 

r- 

in 

CM 6 


•* . 

G 

o 

o o 

o 

o d 

o 

o 

o d 



0.02215 0.01146 0.02074 0.02447 0.06375 


TIME-MEAN AND TURBULENCE DATA FOR SWIRLING 
FLOW = 70°, D/d - 2 , WITH CONTRACTION 


vr 

© 

IP 

rv 


r* 


rv 



© 


IP 

rv 


© 

n 

© 

IP 

(P 

r* 

CN 

IP 

•c 

© 

rv 

CD 

© 

© 

O 


•c 




ip 

o 

f- 

rv 

r- 

(D 

© 

r- 

in 

© 

U) 

rv 

© 


<7 

cn 

« 

LC 

in 


o 

^7 

u> 

r* 

O 


0) 


CD 

© 

p) 

o 

© 

n 

p> 

*7 

in 

in 

© 

in 

tp 

in 

in 

© 

© 

d 

o 

o 

o 

o 

o 

d 

o 

o 

d 

d 

o 

o 

o 

o 

6 


CN — © 

— © r- 

© © © 

O © 

- O o 


v r- 03 

© Oh* 

Gi O LP 

i- O cn 

cm n n 


6 

d 

o 

d 

o 

o 

o 

o 

o 

d 

o 

o 

o 

o 

o 

d 

o 

05 

h- 

© 

cv 

’J 


in 

© 


© 

c 

© 

© 

rv 

© 

O 

r- 

© 

CM 

in 

© 


© 

0) 

© 

© 

in 

o 

r* 

in 

© 

O 

© 

0) 

h- 

O) 

rv 

CN 

r- 

© 

rv 

© 

cv 



O 

rv 

o 

o 

r- 

r- 

N7 

rv 

*- 

rv 


in 

■*- 

O 

■*7 

© 

© 

r- 

o 

rv 

T“ 

*” 

■*“ 

■*“ 





CN 

© 

© 

© 

© 

© 

o 

o 

o 

o 

d 

d 

o 

d 

o 

o 

o 

o 

d 

o 

o 

o 


W CM 

ip cd r- 

oo o cv 

« -c in 

© © © 


o 

0 

1 

0 

1 

0 

1 

o < 

6 

i 

0 

1 

6 

d 

o 

o 

d 

o 


*cr 

© 


© — 

© 

© 

© 

© 

© 

o 

© 

0) 

O 


r- 

— 

r* 

© 

© 

1- 

© 

r** 

© 

o 

r- 

© 

© 

v < 

N 

© 

© 

© 

© 


ID 

© 

r- 

© 

© 

CM 

r- i 

IN 

r- 

rv 

© 

© 

© 

CD 

o 

O 

O 

O 

■*" 

© ( 


r- 

© 

o 

© 

© 



0 

1 

6 

i 

0 

1 

o 

*- i 

!N 

*" 

o 

0 

1 

o 

< 

o 

« 

6 

t 

d 

o 

o 

o 

o 

o 


o 

o 

o 

o 

o 

n 

o 

o 

o 

5 

S 

o 

u 

O 

o 

o 

o 

o 

o 

o 

© 

o 

© 

o 

IP 

0 

m 

o 

© 

o 

© 

o 

© 

r- 

© 

CN 

o 

r* 

n 

CM 

Q 


© 

CN 

o 


© 

© 

r> 

n 

CJ 

N 

CM 

CM 

— 




o 

d 

c 

o 

o 

O c 

D 

6 

o 

o 

o 

o 

d 

d 


138 


0.00000 0.42242 0.3273G 0.33767 0.45736 0.35764 


TABLE VIII (Continued) 



r- 

in 

in 

IP 

r** 


o 

0) 

*7 

h- 

w 

CV 

CO 

tP 

f"* 

CD 

CV 


s 

01 

ip 

t- 

in 

o 

■»— 

«- 

a 

Cl 

*3 

CO 

**• 

^7 

CD 

IP 

Ip 

o 

cv 

in 

CD 

m 

cv 

CO 

IP 

CD 

W 


CV 


CV 

IP 

CD 

CD 

■*? 


to 


IP 

<N 

in 

tn 

CV 

01 

D 

CD 

CD 

CD 

*7 

CD 

CD 

U> 

IP 

cv 

co 

CO 

CO 


CO 

CO 

V 

*3 


*7 

o 

in 

IP 

tn 

in 

U> 



o 

d 

o 

t 

0 

1 

d 

o 

d 

* 

0 

1 

d 

i 

o 

* 

o 

» 

o 

< 

0 

1 

0 

1 

0 

1 

o 

* 

o 


IP 

IP 

CO 

tn 

tn 

in 

o 


CD 

cv 

cv 

CD 

r- 

CO 

*7 

CD 

CV 



cv 

in 

cv 


tn 

in 


■*- 


in 

«0 

cv 

c- 

** 

o 

cv 

ID 

10 


0) 

Cl 

CO 

**■ 

■*7 

r- 

•— 

o 

tp 

tn 

CD 

cv 

W 


0) 


o 


CD 

in 

cv 


r- 

tv 

CO 

<7 

r- 

<0 

*- 

tp 

o 

CV 

<7 


CO 

CV 

CV 

CO 

to 

CO 

CO 

<7 

■«7 

in 

<? 

in 

tn 

in 

in 

to 

IP 


o 

6 

i 

d 

i 

d 

i 

d 

1 

o 

* 

0 

1 

0 

1 

o 

t 

d 

i 

0 

1 

6 

i 

d 

i 

0 

1 

d 

i 

0 

1 

0 

1 


o 

<7 

cv 

o 

IP 

_ 

CO 

CO 

CD 

rv 

cv 

cv 

cv 

in 

in 

ip 



o 

0) 

CV 


to 

o 

CD 

cv 

CD 

cv 

r*- 

o 

rv 

D 

to 

r*- 

in 

O 

o 

CD 

in 

CD 

<7 

(P 

0) 

IP 

h* 

CO 

r- 

<7 

CD 

in 


co 

in 


o 

in 

f* 

0) 

co 


IP 

0) 

r** 

IP 


in 

XX 

CV 

in 

cv 

tp 

o - 

o 

CV 

«N 

CV 

cv 

to 

co 

<0 

*7 

CO 


<7 


n 

to 

CO 

CO 

X 

o 

0 

1 

o 

d 

0 

1 

o 

6 

i 

o 

d 

1 

d 

i 

0 

1 

d 

i 

o 

» 

o 

• 

0 

1 

d 

t 

d 

i 


cv 

0) 

n 

*7 

CO 

0) 

o 

0) 

K 

r- 


oo 

tp 

ip 

r- 

*7 

CO 


ip 

0) 

CO 

in 

o 

CD 

cv 

IP 

CO 

cv 


CO 

to 

tn 

o 

o 

IP 

ID 

to 

0) 

o 

r- 

in 

CO 

o 


CO 


o 

01 

cv 

ip 

in 

0) 

GO 

> 

<p 

IP 

CO 

o 

r- 

f* 

h- 

h- 

-c 

CO 

CV 

CO 

tp 

»- 

cv 

in 

r- 

o 

n 

to 

<0 

V 

*- 

CO 



CO 

^7 

*7 


*7 


CO 

T7 

CO 


6 

o 

o 

d 

o 

o 

t 

d 

d 

0 

1 

0 

1 

0 

1 

0 

1 

6 

t 

0 

1 

O 

» 

o 

j 

d 

i 


n 

U> 


0) 

IP 


(P 

CV 

in 

r- 




!*• 

to 


CD 



in 

IP 

r*' 

^7 

Q 

CO 

in 

«- 


cv 

*7 

in 

CO 


in 

o 

o 

Q 

CO 

IP 

CO 

•— 

o 

in 

ID 

<0 


cv 

in 

£0 

CO 

o 

ip 

0) 


Q 

cv 

IP 

vt 

CV 

CO 




r- 

r- 

tp 


o 

IP 

*7 

tn 

6 

o 

<0 

CO 

<7 

*7 

in 

<P 



u> 

in 

tn 

in 

in 

in 

IP 

r- 


o 

0 

1 

d 

i 

0 

1 

0 

1 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 


o 

O 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 


o 

o 

0 


o 

6 

c 

o 

o 

5 

o 


o 

Q 

5 

o 

o 

Q 

o 

5 

O 


in 

O 

in 

o 

in 

o 

in 

o 

in 

5 

in 

o 

in 

o 

in 

o 

o 


cv 

O 

r- 

in 

cv 

o 

r- 

in 

cv 

o 


in 

cv 

o 


in 

o 

o 

-o- 


10 

CO 

CO 

CO 

cv 

cv 

cv 

cv 


**" 

<T~ 


o 

o 

5 

a 

o 

o 

d 

o 

o 

d 

o 

d 

d 

o 

d 

o 

6 

d 

d 

o 

o 



139 


(b) v/u 


TABLE VIII (Continued) 





1 40 


TABLE VI I I (Continued) 



rms 


TABLE Vi I I (Continued) 



142 


rms 


TABLE VIII (Continued) 




i 

! 

i 

I 

! 

I 




o 

in 


© 



© 



in 


CM 

w 

© 

CM 

© 


© 

O ) 

tn 

CM 

f- 


© 

CV 

<<3 

© 

© 

o 

© 

r- 

in 

o 

© 

o 

r- 

o 

ip 

tv 

o 

© 

© 



tn 

r- 


CV 

tn 

© 


•— 


r* 

r* 

ID 


© 

r- 

0) 

© 

CV 

10 


© 

r- 



*- 

10 

CM 

o 

o 

o 

o 

o 

o 

o 

o 


** 

CM 

CV 

CM 

© 

<3 

*3 

© 


o 

o 

d 

o 

d 

o 

d 

o 

d 

o 

d 

o 

o 

o 

d 

o 

o 




IP 

CD 

© 

V 

© 

t . 



© 


o 


CV 

© 

© 



o 

O 

© 

0) 


h» 

© 

© 

© 

© 

© 

CV 

»* 

© 

tO 


© 

o 

in 

U> 


© 


cv 


r* 


© 

(0 

o 

tn 

h* 

<3 


• 

a 


!*• 

© 


0) 

© 

© 

© 

O 


vj 


V 

*3 

o 

© 


o 

o 

o 

o 

o 

o 

o 

o 

o 




CV 

cv 

cv 

© 

O 


6 

d 

o 

o 

o 

o 

o 

o 

o 

o 

d 

o 

d 

o 

o 

o 

o 



in 

o 

0) 

U> 

o 

© 


® 


VT 


© 

© 

© 

CV 

© 


Q 

VT 

CD 

© 

— - 

5 


r- 

<0 


V 

^3 

cv 

© 

© 

»r— 

© 

o 

o 

CD 


© 

10 

CV 

o 

r- 

© 

© 

in 

in 

© 

o 

r- 

ffi 

© 


o 

© 


© 

© 


© 

© 

© 

CV 

n 

cv 

m 

© 

in 

© 

cv 

o «- 

N 

o 






** 

**■ 



*■* 

** 


cv 

« 

CV 

CM 

X 

o 

6 

o 

o 

o 

o 

o 

d 

o 

o 

o 

6 

o 

o 

d 

d 

o 


in 


01 

o 

tv 

© 

© 

f- 

o 



© 

© 

© 

r- 

*3 



w- 


o 


tv 

<3 

© 

© 


'C 

© 

o 

ID 

CV 



1" 

© 

ffi 

o 

o 

O) 

10 

tv 

CV 

*- 

© 

© 

in 

© 

© 

■<3 

© 

© 


• 

© 

© 

© 

o 

CO 

© 

© 

© 

© 

© 

•*■ 

© 

© 

CD 


r- 

c 

o 


CM 

tv 

CM 

o 


o 

o 

*• 


CM 


CV 

CV 

© 

cv 

cv 


d 

o 

d 

o 

d 

d 

o 

o 

o 

d 

d 

6 

d 

o 

o 

o 

o 



2 



o 

CO 

(0 

© 

t*- 

CM 


I** 

*0 

o 

o 

© 

v- 

© 

© 

© 

© 


Q 

CD 

r- 

VJ 

© 


© 

© 

to 

h* 


© 

© 

© 

© 


© 

o 

Q 

CO 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

*3 


Q 

© 

© 

© 

o 

tv 

r- 

© 

© 

<7 

© 

© 

o 

© 

© 

© 

© 

o 

o 




CM 

© 


© 

© 

V 


© 

© 

© 

© 

© 



o 

o 

o 

d 

o 

o 

o 

d 

d 

o 

o 

d 

o 

o 

o 

d 

6 


D 

N 

a 


o ° 

o 

Q 

o o 

o 

o 

in o 

© 

o 

w o 


© 

■*r 

© 

© 

o o 

o 

d 


8 8 8 8 8 


© 


© 

o 

© 

© 

5 

t- 

© 

© 

o 

© 

© 

© 

© 

o 

o 

o 

d 

o 



o o o o o o o 


143 



rms 


TABLE VIII (Continued) 



0.01424 0.01140 


TABLE VIII (Continued) 


m 


o 

« 


o 
o - 

X 


in 


O' 

ID 

o 

CD 

o 



ID 

CD 

CN 

CN 



CN 

in 

CN 

in 



0) 

r- 


in 

•— 

D 


ID 

re 

(D 

CN 

r- 


0> 

o 



o 

o 



in 

f*> 

in 

CD 


in 

f- 

r* 

in 

N 

0) 

o 

o 

o 

o 

o 

o 

§ 

8 

8 

o 

o 

o 

o 

8 

8 

o 

o 

o 

o 

re 

O 

re 

O 

O 1 

o 

d 

o 

d 

o 

d 

d 

o 

o 

o 

o 

o 

o 

o 

o 

o 

6 

!? 

0) 

*0 

CD 

r- 

in 

CD 



CN 

o 

CD 


re 

ID 


r* 

! 


r> 

0) 

O 

m 

O 

0) 

f- 

ID 

ID 

in 

in 

0) 

in 

CN 

CN 

o 



o 


O 

9* 

CN 

•» 


n 

CN 

re 

CD 

in 

O 

N 


D 

o 

o 

Q 

O 

O 

o 

O 

o 

o 

o 

Q 

o 

9* 


O 

CN 

o 

5 

o 

5 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

d 

o 

d 

o 

d 

o 

6 

o 

o 

o 

o 

6 

d 

o 

d 

o 

n 

CD 

© 




in 

O 

o 

o 



ID 

CD 

ID 

2 i 

o 




r*- 

0) 

n 



CO 

© 

CD 

in 

O 

in 

O 

0) 

6 



n 

CJ 


^7 

(N 

9- 

w 

CN 

CN 

NJ 


N 

<7 

0) 

o 

o 

8 

o 

o 

o 

o 

o 

o 

o 

o 

8 

8 

c 

o 

o 

o 

o 

o 

8 

8 

o 

o 

CN 

o 

.01 

2 

u 

o 

o 

c 

d 

o 

d 

o 

o 

o 

o 

6 

o 

o 

o 

o 

o 

o 

o 

N 


















_5 



LD 

o 

r- 


in 

t- 


CD 

ID 

0) 

CN 

— 

CN 

0) 


3 

in 



O) 

re 

D 

ID 

TT 

re 

r- 

* 

CD 

CN 

o 

re 

o 

— 



CN 

o 

CM 

O 

O 

CN 

CN 


o 

in 

o> 

r- 

r- 

a> 


c 

o 

o 

O 

o 

o 

Q 

o 

O 

O 

o 

O 

O 

•* 

CN 



o 

c 

o 

5 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

O 

O 

o 

XI 

o 

o d 

o 

o 

o 

d 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 



in 

ID 

CN 


CN 

CN 

o 

0) 

IP 

h* 

CD 

CN 


*- 

ft 

in 

in 

CD 

n 

f* 

n 

re 

o 

CD 

U> 

ID 

o 

re 

n 

in 

X 

c 

re 

CN 

n 

LD 

CD 

in 

o 

f» 

in 

0) 

re 

CN 

© 

ID 


o 

f 

o 

o 

O 

o 

0) 

in 

03 

N 

CN 


CN 

CN 



CD 

o 

c 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

O 

O 

d 

o o 

e 

d 

o 

o 

d 

d 

o 

o 

o 

d 

o 

o 

o 

o 

o 


> o 

o 

o 

o 

0 

o 

o 

Q 

o 

p 

o 

o 

o 

o 

{ 

o 

f 1 

> o 

o 

o 


5 

o 

o 

o 

o 

O 

o 

o 

o 

g 

Q 


f 

} in 

o 

in 

ft 

in 

o 

in 

o 

in 

o 

in 

o 

in 

o 

O l 

CN 

> 

5 r- 

IP 

CM 

o 


in 

-N 

o 

r- 

U) 

CN 

o 

r- 

in 

Q 


n n 

n 

re 

n 

CN 

CM 

IN 

CN 





O 

O 

. i 

d 

o d 

o 

o 

o 

d 

o 

o 

d 

o 

o 

o 

o 

o 

d 

o 


145 


h 


TABLE VI I I (Continued) 



X3 

f- 

00 

o 

CM 

in 

o 

ro 

r- 

ID 

ffl 

in 

CD 

f- 

01 

in 

— 



CM 

o 

0) 

® 

CD 

n 

© 

o 

r** 

in 

o 


n 

o 

n 

ID 

o 

n 

CM 

«— ■ 


r> 

CM 

in 


o 

r 

cn 

IP 

r- 

ID 

to 

h- 

in 


o 

O 

o 

Q 

O 

Q 

o 

O 


»- 

— * 

CM 

CM 

CO 

n 

CM 

© 

CM 

5 

o 

o 

O 

O 

3 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 


d 

o 

o 

o 

d 

o 

o 

o 

o 

O 

o 

o 

o 

d 

6 

o 

o 


in 

m 

n 

D 

to 

CD 

o 

o 

in 

in 

CD 

ro 

CD 

C5 

ID 

ID 

in 


to 

CM 

CD 

CO 


0) 


r> 

in 

if 

CM 

© 

•- 

o 


ID 

r* 

in 

n 

CM 

UP 

•— 



CO 

D 



Ol 

m 


tn 

in 

CD 

0) 


o 

o 

O 

o 

o 

o 

CD 

o 

p 

c 

P 

o 

CM 

•- 

— 

»- 

o 

— 

o 

o 

o 

o 

o 

5 

O 

5 

3 

3 

5 

o 

o 

o 

o 

o 

o 


o 

o 

o 

o 

d 

d 

o 

o 

d 

o 

o 

d 

o 

o 

o 

o 

o 


o 

o 

o 

n 

r* 

CD 

h- 

in 

IP 

Q 


n 

in 

o 

IP 

D 

r- . 


o 


ID 


CM 

CM 

rr 

CM 

in 


v 

o 


in 

o 

in 

w 

o 

o 

to 

O) 

to 

ra 

It 

K 


in 

If 


(D 

CD 

r* 

CM 

r> 

to 


o 

o 

o 

o 

ii> 

o 

o 

p 

o 

C 

o 

O 

o 

— - 

— 

— 


o - 

N 

X 

o 

d 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

5 

o 

5 

o’ 

o 

o 

o 

o 

o 

d 

o 

o 

o 

o 

o 

o 

o 

d 

o 

o 

O N 3° 

o j> 


















h 


P*- 

ID 

in 


O) 

n 

cn 

ID 

— 




^3 

to 

CO 

to 

00 !> 


•-* 

t 


in 

in 

o 

to 

00 

n 

<J) o 


ID 

to 

t- 

r> 

in 

in 

*— 


X— 

o 

CM 

CD 

•• 


o 


o 

CM 

P- 

to 


w 


■ 



•— 


o 

o 

O 

o 

o 

c 


p 

«- 

•- 

CM 

U) 

— - '"'* v 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

5 

o 

o 

o 

o 

o 


o 

o 

o 

o 

d 

o 

3 

o 

d 

o 

d 

o 

d 

o 

d 

o 

o 


o 


CM 

0) 

o 


r> 

o 

_ 

o 

_ 

CM 

CD 

© 

ID 

_ 

in 


o 

CM 

o 

o 

<3 

in 

o> 

ID 

0) 

r- 

CO 

in 


<3 

CO 

o 

ro 

o 

o 

CD 

o 


n 

CM 

r> 

ID 

to 

ID 


in 




CO 



o 

O 


w- 

•— 

— - 

in 


m 

ID C 

PM 

CM 

CM 

PM 

n 

•— 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 


o 

d 

d 

6 

o 

d 

o 

o 

o 

c 

o 

d 

o 

o 

o 

o 

o 


o 

o 

o 

o 

o 

o 

o 

o 

p 


p 


o 


p 

o 

2 


o 

o 

O 

o 

5 

o 

o 

y 

5 

c 

5 

Q 

o 

Q 

5 

P 

3 


in 

CD 

in 

CD 

in 

o 

in 

5 

in 

c 

in 

o 

in 

Q 

in 

o 

o 


cm 

o 


in 

CM 

CD 


in 

CM 

c 

r* 

in 

CM 

o 

r- 

in 

o 

o 

* 

* 

r> 

CO 

CO 

CO 

CM 

CM 

CM 

CN 


•** 

— * 


o 

o 

o 

or 

o 

6 

o 

d 

o 

o 

d 

o 

o 

c 

o 

O 

d 

d 

d 

o 

o 


146 


r 

i 


TIME-MEAN ANU ‘i URBULENCE DATA FOR NONSWIRL- 




147 





CM 


n 

CV 

^7 

r- 

ID 

o 



a> 

n 


CD 


CD 

D 

CV 

V 


in 

CD 

IT) 

O) 

CV 

tv 

(V 

in 

CP 

cv 



o 

o 


o 


o 

O 


O 


cv 

CM 

PV 

CM 

CV 

cv 

CV 

CM 

CM 

cv 



6 

o 

6 

o 

o 

d 

6 

6 

o 



W 


0> 


o 

P) 

® 

in 

TJ 



cn 

1- 

o 


D 

P) 

CM 

*7 

CD 


o 


v> 

r- 

ID 

0) 

CD 

05 

CD 

r- 



CM 

t- 

o 

o 

O 


o 

o 

0) 


rt 

CM 


CV 

p; 

CV 

CV 

PV 

tv 

— 



6 

6 

d 

c 

o 

o 

o 

o 

o 



_ 

U> 


in 

o 

CD 

in 

o 

ID 



C 

CO 

CV 

n 

ID 

0) 

CD 

o 

CV 


in 

CD 

in 

to 

CV 

in 

CV 



CV 



a 

o 

*- 

— 

■*- 

o 

P> 

cv 

— 

“O 


— 

CN 

CM 

rV 

cv 

cv 

rv 

CV 

cv 

0) 

o 


o 

d 

6 

d 

o 

d 

d 

o 

o 

c 






















c 


f~ 

CO 

o 

o 

CD 

ID 

in 

n 

® 

o 


D 

«? 

0) 

IT 

o 

CV 

r» 

rv 

PJ 

(_> 

c 

r- 

CO 

in 

n 

ID 

n 

o> 

n 

n 


O 

0) 

o 


O 

— - 

»- 

o 

»- 


o - 

CM 


CV 

cv 

CV 

cv 

cv 

CV 

cv 

X 

s 

X 

6 

o 

o 

O 

o 

o 

d 

d 

o 












LU 













00 

® 

— 

a 

in 

0> 

ro 


® 

<c 


i «> 

co 

o 

*— 

co 

CD 

p* 

w 

CV 

in 

1 ^ 

m 

<3 

in 

o 

a> 

cv 

^7 

cv 

J— 


IT 

o 

n 


O 

CV 


CV 

in 


6 

! « 

tv 

CV 

cv 

rV 

cv 

cv 

cv 

cv 



o 

d 

o 

o 

o 

o 

6 

o 

o 



a ) 

CV 

o 

*7 

VJ 

05 

05 

0) 

_ 



n 

n 

w~ 

in 

o 

r~ 

*7 

n 

in 


o 

r- 

o 

r~ 

CO 


in 

0) 




» 

CO 

O) 


CD 

CD 

c 

cc 

o 



o 


■*“ 


**" 

*■* 

CV 

— 

CV 




o 

o 

o 

o 

6 

o 

o 

o 

o 



o 

Q 


o 

o 

o 

o 

o 

Q 



o 

5 

o 

o 

O 

O 

o 

O 

o 



o 

o 

5 

o 

Q 

O 

o 

o 

o 



o 

in 

o 

in 

O 

in 

o 

in 

o 


o 

V 

r> 

d 

tv 

CV 



O 

o 


\ 



• 

• 

. 

r 





a 

o 

o 

o 

o 

o 

O 

d 

o 

o 



{ 

i - 


i 


* 


148 


i 






1 50 




r* 

«- 

IP 

n 

ip 

0) 

o 

in 

o 

u> 

<n 

r* 

cn 


o 

o 

O) 

co 

ip 

ID 

in 

CO 

CD 

IP 

CD 

in 

ip 

o 

o 

o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

O 

O 

o 

o 

o 

o 

6 

6 

6 

6 

o 

o 

o 

d 

o 


ip 


*- 

CN 

*7 

co 

W 

o 

o 

CT> 


o 

IN 


t?s 

o 


T- 

in 

O 

IP 

ip 

(P 

IP 

ip 

ID 

IP 

(P 

o 

o 

o 

g 

g 

o 

o 

O 

8 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

d 

d 

o 

o 

o 

o 




Oi 

CO 

<7 


CN 


V 

CN 

in 




CN 

co 


o 

o 

CN 

in 

0) 

o 



in 

IP 

10 

IP 

ID 

to 

O 

in 

in 

IP 



o 

o 

o 

o 

o 

g 

O 

o 

o 


~o 


o 

o 

o 

U 

o 

5 

o 

o 

o 


CD 

3 


d 

o 

o 

o 

o 

d 

d 

d 

d 


C 












■M 

C 



r- 

CD 

O' 

in 

CO 

tn 

CO 

CO 

o 

3 

o 


in 

in 

CO 

CN 

n 

CO 


o 

in 


o 

o 

ip 

ip 

IP 

IP 

10 

IP 

IP 

IP 

ip 




o 

o 

g 

o 

o 

o 

g 

o 

o 

£ 


o — 

o 

o 

8 

o 

o 

5 

8 

o 

o 

- u 

X 

X 

d 

d 

o 

o 

o 

d 

o 

d 

o 

3 













ui 

-J 


ip 

CO 

(0 

n 

CN 

n 

CN 

CN 

in 

"aT 

< 


CD 

® 

0) 

a> 

<T 

in 

CD 

IP 

o 


in 

10 

in 

in 

in 

IP 

ip 

in 

m 

ip 


h 

o 

o 

g 

g 

o 

o 

o 

g 

o 



o 

o 

o 

8 

8 

o 

o 

o 

5 

o 




o 

o 

d 

o 

o 

o 

d 

d 

d 






in 

CO 

0) 


0) 

c-* 

CO 




CO 

CO 

r** 


TJ 



CN 




o 

in 

in 

in 

in 

in 

in 

in 

in 

in 




o 

o 

o 

o 

g 

o 

o 

O 

O 



o 

o 

o 

o 

o 

5 

a 

o 

o 

o 




o 

d 

d 

o 

o 

£> 

d 

d 

o 




o 

O 

o 

o 

o 

o 

o 

o 

g 




o 

O 

o 

Q 

o 

o 

o 

o 

O 




o 

o 


8 

8 

o 

o 

o 

o 




o 

in 

8 

in 

o 

in 

o 

in 

o 



D 


CO 

co 

CN 

CN 


* r ’ 

O 

Q 



a 

o 

o 

6 

o 

d 

d 

o 

o 

d 



I 


151 


TABLE IX (Continued) 



152 


TABLE 




153 







! 

in 

o 

to 

n 

in 


o 

CD 

CO 



n 

CO 


O) 

0) 

to 

in 

CM 

in 


1 ^ 

CM 


CO 

0) 


in 

o 

in 

CO 



in 

10 


CM 

CD 

in 

CM 


CM 


cs 

o 


n 

O 

CM 

GO 

in 

to 

O 






— 


6 

o 

d 

o 



in 


CO 

in 

10 

10 

O 


0; 



IT 

*7 


ti 

<0 

o 

CO 

o 

n 


o 

CO 

m 

in 


to 

Q 

o 

CM 

in 



10 

CO 

CM 

10 

in 


CM 

in 

tj 


CM 

o 


«0 


CO 

0) 

10 

* 

o 


1 





■*“ 

d 

o 

d 

d 


1 


0) 

in 

CO 


in 

_ 

m 

o 


j 


r* 

»• 

CD 

in 

CM 

CO 

CP 

to 


in j 

® 

o 

•— 

CO 

CD 

to 

01 

CD 




o 

O 


01 

t- 

CO 

o 

0) 


T5 



CM 

n 


CO 

0) 

to 

o 

CM 

0) 

3 


— 

- 




6 

o 

d 

d 

c 






















*-» 











C 


in 

CO 

in 

■O 

0) 

CM 

o 

o 

to 

o 




(0 

in 

CO 

O} 


CO 

10 

CJ 

o 

n 

r- 

r*- 

in 

t- 

CM 

m 

co 

CM 



U) 


o 

*- 

o 

to 

00 

r- 



o - 


CM 

V 

in 

to 

Ol 

ID 

in 

o 

X 

X 

T-' 

*■» 

j. 


-■ 

o 

d 

d 

o 

UJ 











_J 











<r 




<0 

o 

01 

tn 

r* 

© 

(0 

j± 


01 

Q 

o 

'O’ 


<? 

CO 


o 


in 

Of 

o 

P) 

a> 



CO 

in 

CD 



o 

CM 

CO 

■* 

o 

in 


CO 

O 


6 


CM 

CO 

in 


o 


in 

CM 




•“ 




— • 

d 

o 

d 



© 

CO 

CM 

o 

eo 

iD 

»•> 

to 

*T 




01 

CP 

Q 

-j 

10 

10 

10 

© 


o 

eo 

r- 

CO 

O 

<n O 


01 




CM 

cr 

10 

O 

cn o 

^7 


CD 


6 

10 

CO 

01 

r- 

CM 

o> 

CD 

i> 





*■ 




d 

o 

6 

o 



o 

o 

s 

o 

£ 

g 

o 

g 

g 



8 

8 

c 

8 

G 

o 

8 

o 

p 



O 

in 

c 

in 

c 

in 

o 

in 

o 


■e 

c 

CO 

r 

M 

O 



O 

o 


X 











a. 

o 

6 

c 

» 0 

c 

l o 

d 

o 

o 



j 


I 

I 


155 



* 




157 



158 


I » 




c 

► tx 

— 


*- 

X? 

V 

n 

10 



I s - 

r- 

CN 


n 

w 

r 

i tn 

a) 


in 

n - 

n 

n 

tn 

CN 

o 

i tn 

CN 


• 

c 

' 2 

o 

o 

Q 

o 

c 

> O 

o 


cn 

c 

> o 

o 

o 

O 

o 

c 

l o 

o 



o 

' o 

d 

o 

o 

d 

c 

» d 

d 



in r- 

r- 

tn 

h- 

o 

c 

> x» 

CN 



ip ~ 

in 

r- 

O) 

tP 

cn <p 

r> 


o 

cn n 

N 

tn 

xy 

n 

CN CN 

X? 



o 

' 2 

O 

o 

o 

o 

r 

► 2 

Q 


cn 

c 

1 o 

O 

o 

o 

o 

c 

► o 

5 



d 

i d 

o 

o 

6 

o 

c 

i o 

o 



ID — 

in 

0) 

o 

tn 

r< 


r> 



c 

► n 

o 

XT 

in 

<** 

r> — 

cn 


tn 

n 

i in 

r> 

XT 

in 

n 

cn n 

XT 



c 

1 o 

o 

o 

o 

o 

r 

> o 

O 

•V 

*" 

c 

► O 

o 

o 

o 

o 

c 

1 o 

o 

0) 


c 

i o 

o 

o 

o 

o 

c 

i o 

d 

c 










2 

o 



n — 

ID 

O) 

0> 


u> tp 

3 

CD .N 



o> x? 

n 


n 

XT 

c 

> r- 

IP - 

o 

o 

0) CT) 

t£» 



D 

CM CN 

tn > 

u 


c 

> o 

o 

CN 

*■* 

O 

c 

> O 

CN _ 

"W 

D <- 
\ 

c 

> o 

o 

o 

o 

5 

c 

i o 

Is 

X 

X 

c 

i d 

o 

o 

o 

o 

c 

i o 

o 

LU 










'oi 

CD 


r- 


IP 


r- 

GO 

01 o 

tn 

< 


c 

0) 

'**■ 


tn 

f" 

in co 

h- 

h- 

in 

— ■ 

r- 

r- 

o 

cr> 

XT 

CN CN 

ID 


» 


2 

O 


o 

O 

c 

I O 

O 


o 

c 

5 

o 

o 

5 

o 

c 

► O 

5 



c 

o 

o 

d 

o 

d 

c 

> 6 

d 



n n 

in 


0) 

CO 

IP CN 

r- 



in n 

IP 

CM 

xy 

IP 

r* 

CD 

CN 


o 

r> 


N 

in 


XT 



cn 



c 

o 


CM 


w— 


CN 

»— 


O 

c 

o 

o 

o 

o 

o 

c 

o 

o 



c 

o 

d 

o 

d 

o 

c 

o 

d 



c 

o 

Q 

8 


3 

g 

£ 

g 

§ ’’ 



£ 

o 

tn 

8 

tn 

o 

in 

£ 

in 

O 


o 

XJ 

r> 

n 

<N 

CN 

•* 


o 

o 


\ 


• 

* 

♦ 







a. 

c 

o 

o 

o 

o 

o 

c 

• o 

d 





159 




o 

in 

n 

o 


in 

<P 


_ 



in 

f- 

in 

— 

0) 

in 

ip 

IP 

in 

O 




*“ 



n 

cm 

tM 

CM 




» i 

o 

P 

g 

o 

g 

O 

o 

o 

Q 



CM ! 

o 

o 

o 

o 

o 

o 

o 

5 

5 




6 

o 

o 

o 

o 

o 

o 

d 

o 





<J> 


o 

IP 

0> 

CM 

C3 

IP 




o 

n 

n 


CO 

0) 

!*• 

tJ 

CT) 



O 

CM 

— 

in 

CM 

n 

CM 


CM 

CM 

O 




s 

o 

o 

O 

g 

g 

g 

o 



N 

o 

5 

V 

o 

o 

o 

o 

o 

5 




o 

o 

6 

o 

o 

o 

o 

d 

o 





_ 

CD 


o 

r- 

o 


o 



tn 

IT 

W 

r> 

oi 

CM 

n 

m 

r> 



in 

Q 



* 

o 

o 

o 

g 

Q 

O 

O 

Q 


X 

*- 

o 

o 

o 

5 

o 

o 

o 

o 

5 


0) 

3 


o 

o 

o 

o 

o 

o 

d 

o 

o 


c 






















2 

o 

u 











3 

c 


o 

n 

n 

CM 

IP 

w- 

in 

o 

CD | 

s 

0 


5 

CM 

00 

CM 

*- 


ip 

« 

CM 

“ 

u 

o 



m 

in 

IP 

O 


CM 

U> 



* 

o 

d 

p 

g 

9 

g 

g 

g 

o 



D — 

o 

O 

o 

5 

O 

5 

o 

o 

p 1 

3 

X 

"■h 

X 

o 

o 

o 

o 

o 

o 

d 

6 

o 


LU 











> — . 

_l 











-C 















0) 

IP 

o 

® 

CO 

o 

in 


CM 



in 

r- 

r- 

cm 

IP 

in 

IP 

ip 




i n 

in 


V 

in 

CM 


CM 

in 

O 




g 

p 

o 

g 

g 

o 

n 

o 



o 

o 

o 

o 

5 

o 

o 

o 

o 

5 




o 

o 

o 

o 

o 

d 

o 

d 

o 




o 

IP 

r- 


O 

in 


0> 




o 

r*- 

«T 


n 

IP 


n 

CM 

CM 



IT 

■c 

in 

CD 

03 

n 

CM 

CD 

n 




O 

O 

O 

o 

o 



CM 

CM 



o 

o 

o 

o 

o 

o 

o 

o 

o 

o 




o 

o 

d 

d 

d 

o 

o 

o 

o 




o 

o 

o 

o 

o 

o 

o 

p 

P 




o 

o 

o 

g 


Q 

© 

ci 

5 




o 

o 

o 

o 

5 

o 

Q 

o 

o 




1 o 

in 

o 

in 

o 

in 

O 

in 

o 



o 


rj. 

n 

CM 

cm 

«*■ 


o 

o 






♦ 

» 

* 






a 

1 d 

o 

o 

o 

o 

o 

d 

d 

d 



160 




161 


TABLE XI 


s 




o 

O 1C 

in 

10 

1" 


n 



n 

r- O 

n 

*7 

n 

<0 

O) 


re 

in 

in — 


r- 

in 

CN 

© 


• 

re 

h* in 


re 

0) 

CN 

© 



in 

re 

0 

10 

in 

in 

10 

D 



~~ 


F* 

o 

o 

o 

o 

o 



a 

0) m 

CM 

*•< 

o 

r> 

CN 

C3 


TJ 

CN | 

*- 

in 

0) 

CN 

r** 


z 

o 

CO 

cm n 

10 

o 

re 

CN 



• 


in in 

re 

CD 

<T 

0) 

© 

■J 

«M 

TJ 

« 0) 

10 

*7 

in 

in 

re 

CC 


- 

- 

6 

o 

d 

o 

o 

d 

3 










to 










CC 


CN 

O r» 

01 


10 

0) 

CN 

cv 


n 

in «> 

® 


CN 

m 

in 


IT 

in 

ip 

n 

0) 


CN 

CD 



*3- 

n r- 

CM 

10 

re 

co 

in 



in 

re o> 

ID 


m 

re 

re 

s " 


— 

- 

3 

d 

o 

d 


d 

O -a 










s 










UJ o 


CN 

to in 

CN 

CN 



IP 



to 

in ip 

0) 

CD 


© 

in 

*• 

O 


re in 

D 

(0 

n 

CD 

01 

UJ o 

« 


O tfl 

re 

re 

® 

r* 


_J o 

O 

n 

CM CD 

in 

re 

re 


re 


s 


• 


» 

, 




CD 

X 


- < 

3 

o 

o 

o 

o 

o 

oc -e- 










3 










(- 3 










_ O 


e> 

in mj 


CN 

o 

re 


o _J 


mt 

o n 

O 

CN 

h* 


re 

3S U. 

in 

in 

in oi 

® 

r-* 

n 

CN 

r- 

< 

* 

01 

*7 in 

n 

o 

O 

CN 

re 


o 

n 

CM CV. 

10 

in 

in 

10 

^7 

z 

< 


- 

- 6 

o 

o 

o 

d 

d 

UJ 










s 










UJ 


10 

r- co 

MI 

D 


G 

CO 

Z 



D 10 

CN 

0) 

CO 

CM 

GO 


o 

01 

r- 0) 



re 

re 

CD 


* 

10 

co in 

0) 


re 

CD 

re 


o 

0) 

10 <i 


MI 


M3 

CN 



o 

d d 

o 

o 

o 

d 

o 



o 

o o 

o 


o 


o 



o 

o 

5 

o 

6 

O 

A 

o 



o 

o c 

> 

o 

Q 

o 

Q 

o 



o 

in c 

> 

in 

Q 

in 

o 

in 


D 


re re 

CN 

CN 



o 


s 

♦ 


» 

• 

• 

« 

• 



a 

o 

o o 

o 

o 

o 

o 

o 


r*- 

n 

r* 

0 > 

cn 

6 


r- 

«o 

in 

n 

6 


o 



o 


162 

I . . 





163 







<o 

© 

© 

© 

© 

i** 

•- 

o 

o 


© 

CO 

© 



O 

w- 

CM 


in 

© 

CM 

to 

CM 

CM 

n 

© 

© 

© 


w 

o 

O 

V 

© 

© 


© 


cm 



CM 








d 

o 

o 

o 

6 

o 

d 

6 

o 


« 

o 

o 

© 

© 

© 

V 

© 



CV 

© 


© 

CM 

CM 

CM 

05 


C 

O 

CO 

© 


© 



o 

© 

. * 

in 

© 



CM 

r** 

© 

© 

CD 


** 





*- 





c 

o 

d 

d 

o 

6 

d 

d 

d 



© 

o 



to 

© 

o 




© 

© 

© 

© 

© 

© 

© 


© 

in 

© 


05 

© 

© 

© 

r* 

© 

* 

CO 

O 

h* 

© 

© 




© 



CM 




*• 



© 


o 

d 

o 

o 

o 

o 

o 

d 

d 


0) 

CM 

r- 

© 




© 



© 

*•" 

© 

© 


*— 

o 



o 

© 

CM 


© 

© 

CM 

© 

© 

© 

* 

« 


0) 

f*- 

co 

© 

to 

© 


D *“ 

\ 



*“■ 





*- 

© 

X 

o 

d 

o 

o 

o 

o 

o 

o 

o 


OJ 

o 

© 

© 

© 

o 

© 

© 

© 



© 

© 




© 

© 


in 

CO 



© 

o 

CM 

© 



* 


© 

CM 

o 

© 

© 

© 

r» 


O 



CM 

CM 





© 


o 

o 

o 

o 

d 

d 

o 

o 

o 


CO 

© 

h- 

N 

CM 

o 

© 

© 

© 


o 

CM 

© 

© 

© 

© 

rr 

f- 

© 

o 

tn 

CM 

<0 


© 

o 

CM 

© 

© 



© 


© 

O 

CM 

© 



o 

CO 

n 

CO 

CM 

CO 

CM 

CM 

© 

© 


d 

o 

o 

o 

o 

o 

o 

d 

d 






166 



TABLE XI (Continued) 




tn 

cm 


H 

f- 

CM 

0) 

in 

a) in 


0> 


n 

ID 

o 

id r> 

CD 


w 

ID 


CM 

r\ 

in <D 

« 

o 

O 

8 

8 

8 

8 o 

o 

o 

o 

o 

o 

o 

o o 

o 


O 

n 


in 


o 

CO 

o 

w 

n 

O 

o 

w 

s 

O 

* 

o 

TT 

8 

10 

cm 

n 

8 

in 

tn 

<T 

o 

o 

m 

ID 

O 

a 

ID 

o 

o 

O 

o 

o 

o 

o 

d 

o 

o 

d 

d 

o 

o 

r, 

CO 

n 

0) 


r> 

o 

cm 

c- 

a 

« 

o 

O 

o 

<*> 

o 

® 

8 

CM 

CM 

o 

o 

n 

O 

o 

r* 

in 

8 

r* 

CO 

o 

CM 

o 

d 

o 

d 

d 

o 

o 

o 

o 

o 


O 
O - 
X 


in 

d 


10 

■» 

O 

t- 

CM 

ID 

CM 

0) 

ID 

® 

n 

Q 

o 

w 

CM 

o 

o 

**♦ 

6 

<0 

tn 

ID 

Q 

tj 

tn 

tn 

o 

o 

o 

O 

o 

O 

o 

o 

O 

o 

o 

o 

o 

o 

d 

o 

d 

o 


O) 

tn 

CM 

in 

at cm 

CD 

CO 

0) 

CM 

00 

in *0* 

0) 

CO 

o 

CM 

0) 

•• CM 





O 

"Z 2 

2 

O 

o 

o 

5 

o o 

o 

o 

o 

o 

o 

d o 

d 


oi O 
3 


> 

’D 


0> 


O 

G 


"C 

n 

O 

o 


c 

in 

0) ID 

c 

Oi 

CM 

IP 

co 

co n 

CM 


(X> 



*“ 

ID 


0) 


CM 

CM 

«r- 


CM 

O 

O 

o o 

o 

c 

o 

d 

o 

o o 

6 

o 

o 


GO 

in 

« 

o 


o 

a 


8 

8 


o 


s 

tn 

§ 

o 

in 

8 

CM 

o 

in 

o 

C) 

CM 

** 

o 

o 

d 

d 

d 


O Q 
O O 
o 5 

8 8 



l 





168 





169 


V « 



EDDY DISSIPATION RATE FOR NONSWIRLING 
FLOW 4 » 0°, 0/d = 2 



cT/D 


TURBULENT DISSIPATION LENGTH SCALE FOR NON- 



EDDY DISSIPATION RATE FOR SWIRLING 



TURBULENT DISSIPATION LENGTH SCALE FOR 



ooooo o 



TURBULENT DISSIPATION LENGTH SCALE FOR 
SWIRLING FLOW <J> = 70°, D/d = 2 





p- 

o 

n 

0) 

P 

a ) 


CD 


in 

p 

o 

0) 

f- 

<? 

'X 

«r- 


1 


CN 

V 

<7 

id 

<D 

t— 

m 

CO 

in 

in 

CN 

p 

in 

CO 

r 

n 


in 

p* 

P* 

in 

rj 

n 

CN 

re 

co 

■r- 

0) 

CN 

w- 

o 

CD 

CN 




• 

to 

ID 

CD 


id 

<N 

re 

01 

<r 

P) 

in 

CD 

CN 

n 

in 




cn 

o 

o 

o 





o 

o 

o 

o 

O 


*— 


*— 




d 

o 

d 

o 

o 

o 

d 

o 

d 

d 

o 

o 

o 

o 

o 

o 

o 



co 

0) 

0) 

I- 



P> 

CD 

(X) 


CO 

CN 

P 

CN 

CO 

0) 

m 



01 

n 

in 

CN 


p 

p* 

in 

0) 

CD 

ID 

in 

0) 

O 

<7 

If 

^7 


o 

cn 

CN 

f- 

CM 

o 

CN 

o 

p* 

*7 


in 

p 

CO 

CD 


c 

CD 



r- 

'T— 

•v- 

n 

0) 

9 

n 

CN 

<N 

CN 

CN 

n 

*r 

in 

in 

*7 

ro 


CN 

o 

o 

O 

o 

o 

o 

p 

p 

o 

o 

o 

o 

o 

o 

o 

o 

o 



o 

o 

o 

o 

d 

o 

o 

o 

o 

o 

d 

d 

o 

o 

d 

o 

d 



*7 

P* 

Ol 

p . 

01 

cn 

<0 

o 

o 

in 

P 

Cl 

CD 

P> 

P) 

o 

in 



6 

<D 

O 

00 

id 

in 


o 

CP 

t- 

CD 

CN 

CD 

PI 

01 

IX 

CN 


in 

5 

C0 

O 

o 

^7 

p- 

*0 

p 

CN 

o 

01 

* 


»- 

r- 


«•— 




O 

*r* 

•— 


CN 


CN 

t* 

CN 

CN 

D 

P) 

P) 

pi 

P 

O) 


**■ 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

c 

o 

o 

o 

o 

p 



o 

o 

o 

d 

o 

o 

o 

d 

o 

d 

o 

o 

o 

o 

d 

o 

d 



CO 


o 

in 

in 

P- 

CN 

o 

CN 

o 

P 

CD 


CD 

^7 

in 

a 



0) 

<7 

P) 

0) 

pi 


re 

<7 

P 

in 

O 

o 

CN 


in 

p- 



o 

o 

o 


n 

p* 

PI 

id 

00 

o 

p 

p 

n 

o 

O) 

p* 

0) 

CD 



CN 

P) 

n 

CN 

cn 

n 

0) 

re 

<7 

<? 


in 

CD 

CD 

in 


CN 

o 

a - 

o 

o 

o 

o 

o 

p 

p 

o 

o 

o 

o 

o 

o 

O 

o 

C 

o 







* 














X 

d 

o 

o 

d 

o 

o 

o 

o 

o 

o 

o 

d 

o 

o 

o 

o 

o 



P- 

r* 

0) 


IP 


C 

CD 

o 

IT 

30 


CN 

IP 

CN 


CD 



o 

CN 

o 

P> 

CO 

to 

CN 

in 


Pi 

01 

P 

n 

CO 


CN 

CO 


in 

CN 


CN 

01 

0) 

01 

*y 


CO 

CD 


*7 

0) 

O 

IP 

■y 




CN 

CM 

CN 




CN 

CN 

CN 

CN 

C) 

P) 

p> 

^7 

P) 

P 

CN 


6 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

a 

o 

o 

o 

o 

o 

o 



o 

o 

o 

o 

d 

6 

6 

d 

d 

d 

o 

d 

o 

o 

o 

a 

o 



o 

*7 

01 

o 

o 

P- 

tn 

CN 

n 

re 

CD 

PI 

Pi 1 

in 

P) 

•7 

CN 



0) 

m 

P* 

CO 

CO 

o 

CN 

P 

01 


CN 

CD 

CD 

in 

in 

cx 

a) 


o 

o 

O 

o 

y 

p 


o 

g 

o 

Q 

O 

Q 

Q 

o 

O 

r 

Q 



q 

o 

o 

5 

o 

o 

p 

o 

o 

o 

p 

o 

o 

o 

o 

c 

o 



d 

o 

o 

o 

d 

o 

o 

o 

o 

o 

o 

d 

d 

d 

o 

o o 



o 

c 

o 

o 

o 

o 

O 

o 

o 

o 

g 

o 

n 

Q 

n 

c 

n 



o 

u 

O 

o 

p 

o 

5 

o 

o 

o 

ft 

C5 

g 

ft 

o 

ft 

Q 



y 

in 

o 

in 

p 

in 

o 

in 

o 

in 

o 

in 

o 

in 

o 

IT 

O 



o 

p- 

in 

CN 

o 

t- 

in 

CN 

u 

p 

in 

CN 

u 

p- 

tn 

P 

Q 


Q 


co 

n 

re 

n 

<N 

CN 

CN 

CN 





o 

o 

c 

o 


\ 

• 







• 




» 

■ 






Oc 

o 

6 

6 

d 

6 

o 

o 

o 

d 

o 

o 

o 

o 

o 

o 

c 

o 



176 



TABLE XVI I I 


<&u><0ipu>h»r'Or> 

«§§ 8 §§§§§§ 

666666666 


ICLPLDLOI^CDCOOO 

ooodododd 


to ri tn in in.in in in is 
inPOOOpOQpc 

* Q pooogogo 

O'- Ooooooooc 

x 666666666 


tD IP IP tp 


r- rN» ® 


888888088 

000000000 

odo dddddd 


CV 

cv 

CV 

cv 

CV 

CV 

CV 

0 

CV 

CV 

O 

u 

0 

Q 

p 

0 

O 

Q 

0 

0 

O 

0 

9 

Q 

0 

Q 

5 

0 

0 

O 

0 

0 

0 

0 

O 

O 

0 

d 

d 

d 

d 

d 

d 

d 

0 


O 

0 

O 

0 

O P 

p 

O 

5 

O 

0 

O P 

0 

O 

0 

O 

Q 

O P 

0 

ip 

0 

in 

5 

in 0 

in 

r> 

n 

cv 

CN 

■»- 

O 

6 

d 

d 

d 

d d 

0 




cc 

o 


X 

X 


m 


«! ■ 
O TJ 
N 
Q 
X 

i— ► 
O o 

Z o 

UJ 

-I II 

z -e- 

0 

h- § 
< —i 
CL U. 

t/> cs 

01 z 

Q 3 

ad 
I- — 

2 3 

UJ V) 

-I z 
=5 O 

S 38 


o 

1 o 
o 

I ^ 

o 

© 

o 

CN 

© 

Ifi 

© 

O 

cn 

© 

© 

CN 

© 

o 

© 

f"* 

tn 

in 

© 

© 

© 

cn 

o 

0 

8 

o 

8 

8 

8 

8 

8 


o 

o 

d 

o 

c 

o 

d 

o 

d 

O 

o 

u 

O 

Q 

© 

N 

CN 

© 

O 


© 

© 

© 

n 

o 

© 

in 

© 

LO 

© 

h- 

CN 

o 

O 

5 

o 

8 

8 

8 

8 

8 


o 

o 

d 

c 

o 

o 

d 

d 

o 

in j 

*£ 

in 

■ji 

00 

O 

f'c 

cn 

© 

© 

o 

<y 

n 

© 

o 

C~ 

CM 

© 

CN 

© 

in 

n 

m 

© 

C5 

w 

o ^ 

N 

O 

O 

8 

8 

§ 

O 

8 

8 

8 

X 

O 

o 

d 

o 

d 

o 

o 

o 

o 

in 

<7 

to 

© 

o 

© 

© 

K 

f"v 

!-~ 

® 

8 

0) 

in 

© 

© 

© 

© 

© 

o 

© 

© 

© 

8 

d 


8 

8 

8 

8 

8 

o 

o 

8 


o 

o 

o’ 

o’ 

d 

d 

o 

o 

d 

O 

CJ 

CN 

© 

CN 

o 

>T 

r*. 

t 

n 

© 

n 

fN 

© 

© 

r> 

© 

« 

o 

© 

O 

M 

5 

o 

O 

8 

8 

O 

O 

o 


d 

o 

o 

d 

o 

o 

o 

o 

o 


Q 



! 

i 

i 


178 


TABLE XX 




180 


TABLE XXI I 


r 



181 



TABLE XXIV 


< 

Z O 
— CO 

(/> 

Ui U_ 
— O 

— Ui 
H — i 
Z 13 

< z 

g< 

Ui 
UJ Z 
o < 
z > 

Ui 

-J —i 
3 o: 
00 _ 

£ ^ 
3 CO 
K* 

31 

Z I— 
o — 

ce 3 
ui o 
fr— — t 
Ui u 
21 

CD 


CO 

o 

a. 

o ^ 

-J LA 
UI CSj 

LL O 

3 

c II 

-J 

Ll O 

Ui u 

> 


2 

< 

a. 


3 3 

Cl CO 
Z 

— o 

Ui 
Ll. Z 

o — 
u 
u- Z 
C-» o 
ui o 
Ll. 

U. 


H II 

z 

ui o 
co \ 
ui x 
a: ' 

CL 
UJ 

cn 


5 

=> 


& 
- w 
* 


-E 


i * 


i> 


55 


§ 

3 


♦ 


luT |uT 


» 

<* 

w 

in 

I 


O 

'•> 


to 

+ 


I^J 


O - -H 


ur 


m — • »-« 


a 

o 


09 

o 

«4 

I 


©* 

c 

*■4 

# 


o 

N 


CD 

ei 

o 

i 


o 

o 


o 

o 

♦ 


^ 


N — « 
° ° 
-« o 
» ♦ 


(O 


o 

N 

I 


rf' rf- 


183 


♦0.21 +0.01 +0.05 -1.63 +0.13 -0.79 



TABLE XXV 


MEASUREMENTS FOR SITUATION A, CASE 1 


Combination 

Used 

u /u 0 

Measured 

v/u 0 

w/«o 

612 

0.978 

0.203 

0.009 

123 

0.961 

0.199 

NR 

234 

0.970 

0.126 

0.020 

345 

0.968 

NR* 

0.013 

456 

0.976 

0.011 

0.008 

561 

0.967 

0.197 

0.017 

Mean x 

0.970 

0.147 

0.013 

S.D. o 

0.006 

0.074 

0.005 

o/x 

0.0C6 

0.503 

0.353 


* Not Resolved 
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TABLE XXVI 


MEASUREMENTS FOR SITUATION A, CASE 2 


Combination 

Used 

u/u 

0 

Measured 

V/U 0 

w/u o 

612 

0.978 

0.031 

0.035 

123 

0.972 

0.036 

0.040 

234 

0.990 

0.056 

0.015 

345 

0.982 

0.043 

U .022 

456 

0.988 

0.024 

0.049 

561 

0.986 

0.022 

0.041 

Mean x 

0.983 

0.021 

0.033 

S.D. o 

0.006 

0.012 

0.012 

o/x 

0.006 

0.553 

0.353 
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TABLE XXVII 


MEASUREMENTS FOR SITUATION A, CASE 3 

n 


Combination 

Used 

u/u 

0 

Measured 

v/ u o 


612 

0.954 

-0.030 

-0.008 

123 

0.954 

-0.031 

-0.008 

234 

0.967 

-0.039 

-0.051 

345 

0.948 

-0.103 

0.080 

456 

0.971 

0.001 

-0.003 

561 

0.951 

-0.030 

0.006 

Mean x 

0.958 

-0.039 

0.001 

S.D. o 

0.008 

0.031 

0.039 

a/x 

0.009 

0.805 

38.986 
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TABLE XXVIII 


MEASUREMENTS FOR SITUATION A, CASE 5 


Combi nation 
Used 

u/u 0 

Measured 

v /u 0 

w/u 

0 

612 

0.976 

0.203 

0.010 

123 

0.960 

0.200 

NR 

234 

0.971 

0.126 

0.019 

345 

0.968 

NR 

0.013 

456 

0.976 

0.010 

0.008 

561 

0.968 

0.198 

0.017 

Mean x 

0.969 

0.147 

0.013 

S.D. o 

0.005 

0.074 

0.004 

o/x 

0.006 

0.503 

0.317 
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COMPARISON OF AVERAGED RESULTS FROM SIX POSSIBLE WIRE 
COMBINATIONS WITH OPTIMAL COMBINATION 







Figure 2. The Flowfield Being Investigated 








8. 0096 ♦ 3.9943U* - 0.0B2HJ 




PITCH, YAW FACTORS 



HOT-WIRE VOLTAGE E(VOLTS) 

Figure 6. Pitch and Yaw Factors as a Function of 
Hot-Wire Mean Effective Voltage 
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Figure 9. Free Swirling Jet Measurement Locations 
for Situations D and E 




Figure 10., Free Jet Calibration Apparatus 
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Figure 16 . Time-Mean and Turbulent Flowfield <j> 
Degrees, D/d = 2 


















Figure 17. (Continued) 
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Figure 19. Time-Mean and 
Wi th Strong 
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Figure 22. (Continued) 
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